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ABSTRACT 

The  electronic  interactions  responsible  for  the  observed  electronically 
induced  phase  transformations  in  several  important  classes  of  materials  were 
investigated  theoretically  in  close  collaboration  with  experimental  programs 
at  other  institutions.  These  include  destruction  of  superconductivity  and 
onset  of  a  normal  magnetic  state  at  low  temperatures  in  ternary  rare-earth 
superconducting  compounds  and  the  possible  occurrence  of  a  new  mixed  phase, 
the  origin  of  the  phase  transition  to  the  superconducting  phase  in  the  Chevrel 
compounds,  the  possibility  of  excitonic  superconductivity  in  CuCf,  the  origin 
of  phonon  anomalies  in  Pd  and  Pt  metals,  electronically  driven  phonon  insta¬ 
bilities  and  the  formation  of  charge  density  waves  in  layered  transition 
metal  dichalcogenides 
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I .  INTRODUCTION 

This  Final  Technical  Report  under  AFOSR  Grant  No.  76-2948  covers  a  highly 
prolific  period  of  solid  research  accomplishments  marked  especially  oy  their 
forefront  and  innovative  character.  New  theoretical/computational  methods  were 
developed  and  applied  to  the  study  of  electronically  driven  phase  transitions 
in  a  wide  variety  of  materials.  To  report  in  detail  on  all  these  accomplishments 
would  require  a  massive  tome.  Since  four  interim  annual  technical  reports  have 
already  been  submitted,  this  Report  will  summarize  progress  and  accomplishments 
achieved  on  several  projects  in  this  last  year. 

The  high  level  of  productivity  is  seen  from  the  appended  total  list  of 
publications  for  the  entire  period  -  numbering  some  50  publications  in  all. 

This  is  but  one  measure  of  the  accomplishments  and  recognition  achieved.  It 
does  not  include  the  numerous  invited  (and  contributed)  papers  presented  at 
national  and  international  conferences,  invited  seminars  and  colloquia  presented 
at  universities  and  national  labs,  international  conferences  and  seminars  which 
the  PI  helped  to  organize  (often  as  chairman  or  co-chairman) ,  service  on 
national  and  international  scientific  committees,  etc. 

II.  RESEARCH  OBJECTIVES 

The  aims  and  purposes  of  this  research  were  to  (i)  develop  sophisticated 
theoretical  and  computational  methods  which  would  (ii)  permit  the  accurate 
theoretical  determination  of  the  elec  r  -  ic  structure  of  complex  systems  and 
the  origin  of  their  electronically  driven  phase  transitions.  It  is  widely 
recognized  that  the  many  important  and  increasingly  complex  properties  of  such 
materials  require  the  most  developed  form  of  the  quantum  theory  of  matter. 

Unique  information  derived  from  such  studies  is  increasingly  attracting  the 
interest  of  material  scientists  seeking  solutions  and/or  understanding  of 
phenomena  (including  some  in  the  more  classical  metallurgical  domain)  thought 
previously  to  be  too  complex  for  first  principles  theoretical  treatment. 
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III.  METHODS  AND  APPROACH 

The  newly  discovered  exotic  phenomena  in  complex  materials  such  as  the 

rare-earth  ternary  superconductors  have  challenged  the  traditional  understanding 

« 

acquired  from  approximate  theoretical  treatments  within  the  quantum  theory  of 
matter.  This  challenge  has  been  met  by  recent  developments  of  the  energy  band 
method  which  have  made  it  an  increasingly  powerful,  sophisticated,  and  successful 
tool  of  the  solid  state  physicist.  The  success  of  the  energy  band  approach  is 
based  on  the  work  of  Hohenberg,  Kohn  and  Sham  in  developing  the  local  density 
functional  formalism  for  investigating  various  ground  state  properties  of  solids. 
Aside  from  providing  a  rigorous  basis  and  justification  of  the  single  particle 
energy  band  description,  it  has  led  to  accurate  tractable  computational  schemes 
for  describing  them  from  first  principles. 

Because  of  the  importance  of  charge  transfer  between  different  atom  species, 
meaningful  results  are  possible  only  if  carried  out  self-consistently.  For 
magnetic  systems,  local  spin  density  functional  theory  carried  out  as  self- 
consistent  spin-polarized  energy  band  calculations  (different  potentials  for 
different  spins)  are  obviously  necessary.  Both  self-consistency  (which  requires 
iterative  procedures  to  obtain  convergence)  and  spin  polarization  make  these 
calculations  enormously  more  complicated  and  time  consuming  and  so  have  been 
carried  out,  using  traditional  methods,  only  on  simple  systems. 

Our  work  has  been  at  the  forefront  of  these  developments  in  the  local  density 
theory-energy  band  approach.  Several  novel  schemes,  such  as  the  self-consistent 
linear  combination  of  atomic  orbitals  (LCAO)  method  and  the  self-consistent 
linearized  augmented  plane  wave  (AFW)  and  the  self-consistent  spin-polarized 
linear  muffin  tin  orbital  (LMTO)  methods  are  now  among  the  basic  methods  we  have 
developed  and  are  using  in  our  fundamental  theoretical  studies  of  both  bulk 
solids  and  surfaces.  The  examples  described  later  under  progress  indicate  some 
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of  the  richness  and  variety  of  properties  and  phenomena  which  are  being  and 
can  be  studied  today.  Comparisons  with  experiment  provide  today  such  a  demanding 
test  of  the  predictions  of  one-elecLron  theory  that  our  work  has  the  additional 
virtue  in  permitting,  by  the  comparison  with  experiment,  'accurate  determinations 
of  the  relative  magnitude  and  importance  of  many-body  effects  in  real  solids. 

IV.  SCIENTIFIC  SIGNIFICANCE  OR  IMPACT  ON  TECHNOLOGY 

The  important  phenomena  of  phase  transitions,  superconductivity  and  magnetism 
in  the  complex  systems  being  invescigated  are  highly  exciting  from  the  scientific 
point  of  view.  In  addition,  these  materials  are  also  important  technologically 
in  that  they  have  potential  device  applications,  and  so  the  theoretical  knowledge 
obtained  from  the  underlying  physical  phenomena  may  have  important  technologic-'1 
consequences.  To  give  one  example,  the  Chevrel  phase  ternary  molybdenum  chalcogenides 
are  examples  of  compounds  with  remarkable  physical  properties.  Not  only  do  they 
present  the  highest  superconducting  transition  temperatures  found  for  any  ternary 
system  so  far  (~14K) ,  but  they  also  have  critical  magnetic  fields  (~600  kG)  which 
are  by  far  the  highest  ones  ever  observed.  This  latter  property  has  created 
the  main  interest  in  these  materials  because  they  could  be  suitable  for  certain 
technological  applications  such  as  the  production  of  very  high  magnetic  fields. 

Some  of  their  properties  are  related  to  their  unusual  phonon  spectrum  or  to  the 
extreme  pressure  sensitivity  of  their  transition  temperatures  and  others  to  the 
occurrence  of  lattice  instabilities  at  low  temperatures  which  are  often  believed 
to  be  closely  related  to  superconductivity  in  general.  Hence  the  additional 
excitement  caused  by  the  observation  that  the  addition  of  magnetic  impurities 
like  Eu  causes  a  large  increase  in  the  critical  field.  Most  recently,  major 
interest  has  focussed  on  the  magnetic  behavior  at  low  temperatures  and  its 
influence  on  superconductivity  including  (i)  the  possible  coexistence  of  super¬ 
conductivity  and  magnetic  order  at  low  temperatures  and  (ii)  the  observation  of 
reentrant  magnetism  from  the  superconducting  state  back  to  the  normal  state  at 
low  temperatures.  Our  own  work  has  given  the  first  understanding  of  the  complex 
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phenomena  associated  with  the  origin  of  magnetism  and/or  superconductivity 
in  these  materials  and  appears  to  explain  the  very  high  upper  critical  fields 
observed  in  these  systems. 

( 

V.  RECENT  PROGRESS  ACHIEVED 

A.  High  Field  Superconductivity  of  Chevrel  Phase  Compounds 

The  properties  of  the  Chevrel  phase  compounds,  particularly  SnMo  S  and 

'68 

EuMo  S  ,  are  of  great  scientific  and  technological  interest  because  of  their 

O  O 

unusual  high  magnetic  field  properties  including  the  enhancement  of  the  upper 

critical  field  H  observed  when  Eu  is  added  to  SriMo  S  .  For  example,  up  to  a 
c  L  oo 

Eu  concentration  of  0.5,  T^  is  hardly  changed  with  the  depression  occurring 
abruptly  only  at  high  concentration.  This  behavior  is  contrary  to  observations 
on  all  other  materials  with  the  addition  of  a  local  magnetic  moment  and  the 
theory  of  Abrikosov  and  Gor'kov. 

We  have  carried  out  self-consistent  LPTCO  energy  band  calculations  for 

these  compounds  including  all  electrons  in  all  15  atoms  per  unit  cell.  As  a 

result  of  these  unique  studies,  we  have  demonstrated  that  supercond”ctivity  is 

due  to  the  high  Mo  d-band  density  of  states  at  the  Fermi  energy  resulting  from 

the  unusual  large  charge  transfer  of  Mo  electrons  to  the  chalcogens.  There 

is  also  a  large  charge  transfer  from  the  metal  site  to  the  cluster  giving  Eu 

essentially  no  occupied  conduction  bands  and  for  Eu  a  typically  divalent  isomer 

shift  in  very  good  agreement  with  the  experiments  of  Dunlap,  e_t  al^.  Furthermore, 

we  find,  also  in  close  agreement  with  their  spin  lattice  relaxation  rate 

measurements,  that  the  conduction  electron  density  of  states,  and  hence  [jN(E  )] 

r  , 


is  reduced  by  an  order  of  magnitude  from  its  metallic  state  value.  This  results 
in  a  very  weak  coupling  of  the  4f  electrons  to  the  conduction  electrons  and  hence 
only  a  very  weak  Ruderman-Kittel  interaction  which  leaves  the  dominant  magnetic 
interaction  to  be  the  weak  dipole  interaction  surmised  by tRedi  and  Anderson. 

B.  Isos tructural  Phase  Transition  in  fee  La  and  High  Superconducting 
Transition  Temperature 

Fee  La  shows,  at  low  temperature  under  pressure,  an  apparently  isostructural 
phase  transition  at  ~53  kbars.  The  importance  of  understanding  this  transition 
is  related  to  the  importance  of  the  observed  high  superconducting  transition 
temperature,  T  ,  and  its  unusually  large  pressure  derivative,  dT^/dp.  This 
latter  observation  has  raised  wide  speculation  that  a  new  pairing  mechanism 
involving  f  electrons  was  operative. 

We  have  been  carrying  out  extensive  studies,  as  proposed  last  year,  of  the 
electronic  structure  and  properties  of  fee  La  at  ambient  and  high  pressures  by 
means  of  a  self-consistent  relativistic  LAPW  study  of  the  band  structure  of 

fee  La  at  three  values  of  the  lattice  constant.  An  important  feature  of  this 
approach,  especially  in  view  of  the  care  which  is  necessary  when  treating  f 
states,  is  that  few  restrictions  (i.e.  constraints)  are  built  in;  for  example, 
the  charge  density,  and  the  resulting  self-consistent  potential,  are  a1 lowed 
a  general  (i.e.,  non-muff in-tin)  form,  and  important  relativistic  effects  are 


included  either  exactly  (mass-velocity  and  Darwin  terms)  or  non-perturbative ly 
(spin-orbit  interaction). 


We  find  that  the  f  bands  lie  ~2-2.5  eV  above  the  Fermi  level  and  are 

~1  eV  wide,  resulting  in  a  very  small  (0.05  electrons)  localized  f  occupation. 

Under  pressure  the  f  bands  rise  and  broaden  appreciably,  resulting  in  only  a 

slight  increase  in  f  occupation.  The  rigid-muff in- tin  approximation  for  the 

electron-phonon  interaction,  X,  overestimates  the  superconducting  transition 

temperature  T  by  407.,  but  we  find  that  the  drastic  increase  in  T  under 
c  c 

pressure  can  be  attributed  primarily  to  changes  in  the  electronic  stiffness. 
Structural  transitions  which  occur  at  25  kbar  and  53  kbar  may  be  related  to 
changes  in  Fermi  surface  topology  which  we  find  to  occur  approximately  at  these 
pressures.  Indeed,  the  theoretically  observed  behavior  of  the  bands  under 
pressure  strongly  suggests  a  Lif shitz-Dagens  structural  instability  due  to  a 
strongly  d-f  hybridized  band  crossing  E  in  both  pressure  ranges. 

r 

C.  Origin  and  Nature  of  the  a-v  Phase  Transition  in  Cerium 

As  discussed  in  our  proposal  last  year,  cerium  metal  has  a  number  of 
unusual  but  not  yet  well  understood  properties  which  are  thought  to  arise  from 
the  nature  and  character  of  its  4f  electrons.  These  unique  properties,  which 
include  a  phase  diagram  with  at  least  five  solid  phases,  the  occurrence  of 
antiferromagnetism  at  low  pressures  and  superconductivity  in  the  high  pressure 
phases,  have  stimulated  much  experimental  and  theoretical  work.  Theoretical 
efforts  have  been  mainly  concerned  with  model  investigations  of  the  iso- 
structutai  o-y  transition.  This  work  has  emphasized  the  local  character  of  the 
4f  states  while  describing  the  band  structure  by  relatively  crude  models. 

We  have  undertaken  relativistic,  self-consistent  linearized  augmented 
plane  wave  (LAPW)  band  structure  calculations  for  fee  cerium  at  various  values 
of  the  lattice  constant  equal  to  or  greater  than  that  of  o-cerium  in  order  to 
study  the  isos tructural  ar-v  phase  transition  which  occurs  under  pressure 
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In  contradiction  to  the  4f  electron  promotional  model  of  the  transition,  the 

results  yield  essentially  one  4f  electron  to  be  occupied  in  each  phase  but  with 

the  4f  wave  function  somewhat  less  localized,  and  therefore  more  bandlike,  in 

the  "collapsed"  a  phase.  A  singly  occupied  4f  state  is  shown  to  be  consistent 

with  the  available  experimental  data.  The  weaker  bonding  of  the  4f  electrons, 

compared  to  that  of  the  6s-5d  valence  electrons,  accounts  for  n-ceriun  appearing 

to  have  3. 5-3. 7  bonding  electrons.  Our  calculation  of  the  superconducting 

transition  temperature,  T^,  suggests  that  a  small  spin  fluctuation  contribution 

detrimental  to  superconductivity  is  necessary  to  account  for  the  very  lew  value 

of  T  in  a-Ce. 
c 

D .  Surface  Phase  Transitions  in  '.'(OOP 

The  unusual  phase  transition  of  the  (001)  surface  of  tungsten  metal  has 

been  investigated  using  results  of  our  calf -consistent  calculation  for  the 

tungsten  (001)  surface,  obtained  by  means  of  our  recently  developed  self- 

consistent  linearized-augmented-plane-wave  (LAPU’)  method  for  thin  films. 

Above  T  ~  300K  the  ideal  (lxl)  surfaces  are  observed,  and  the  reconstructed 
c  ~ 

c(2x2)  and  c(2.2x2.2)  superlattices  are  observed  for  W  and  Mo,  respectively, 

below  T  . 
c 

Our  work  has  centered  on  the  generally  accepted  view,  which  we  have 
developed,  that  instabilities  in  the  electronic  ground  state  with  respect  to 
a  spatially  inhomogeneous  perturbation  are  associated  with  sharp  structure 
in  the  electronic  generalized  suscepti biiicy  function.  In  bv1v  e,rs terns,  such 
structure  has  been  related  to  the  existence  of  parallel  pieces  of  Fermi  surfaces 
(nesting  features)  and  is  responsible  for  the  occurrence  of  span  density  waves 
in  chromium  metal  and  spiral  spin  structures  in  the  rare-earth  metals.  Using 
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our  self-consistent  results  for  W(100)  we  have  calculated  the  layer-projected 
generalized  s’- “  '-eptibility  function  to  detect  such  incipient  instabilities.  Ke 
find  th*.  .lectronic  surface  states  or  surface  resonance  states  near  the  Fermi 
energy  appear  to  drive  the  phase  transition  through  a  surface  phonon  softening 
and  a  gapping  of  the  two-dimensional  Fermi  surface.  This  provides,  for  the 
first  time,  theoretical  evidence  chat  the  phase  transition  may  indeed  be  elec¬ 
tronically  driven  —  in  agreement  with  the  earlier  proposal  of  P.W.  Anderson 
(Nobel  Laureate)  and  Tosatti. 
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Relative  isomer  shifts  of  3 d,  Ad.  and  Sd  ions  are  obtained  by  calculation  of  atomic  electron  charge  densities, 
using  numerical  Hartree-Fock  and  multiconfi duration  Dirac-Fock  theory,  with  exact  exchanee.  to  describe  the 
ionic  systems.  The  relativistic  calculations,  for  both  point  and  finite  nucleus  models,  are  compared  with  results 
obtained  from  the  nonrelattvistic  calculations  augmented  ov  one-e.ectron  reiat.  •  >:ic  corrections.  The  rela¬ 
tivistic  charge  densities  and  experimental  isomer  shifts  are  employed  to  deduce  estimates  of  changes  in 
mean  square  nuclear  radu. 


L  INTRODUCTION 

With  the  discovery  of  the  MQssbauer  effect,  it 
was  widely  recognized  that  the  measurement  of 
isomer  shifts  had  great  potential  for  becoming  an 
important  tool  for  the  study  of  both  nuclear  and 
electronic  properties  of  atoms  in  various  en¬ 
vironments.1  The  isomer  shift  provides  a  mea¬ 
sure  of  the  product  of  the  difference  of  mean- 
square  radii  of  two  nuclear  isomers  and  the  dif¬ 
ference  in  electron  density  in  the  region  of  the 
nucleus,  between  two  isomeric  atoms  in  different 
chemical  environments.  Whereas  these  quantities 
are  impossible  to  measure  by  themselves,  iso¬ 
mer-shift  measurements  coupled  with  independent 
information  about  either  the  isomeric  radii  or  the 
electronic  densities  provides  knowledge  of  the 
other.  Thus,  in  principle,  a  measurement  of  the 
isomer  shift  between  two  tons  in  different  en¬ 
vironments,  and  an  accurate  calculation  of  the 
electronic  charge  densities  would  allow  one  to 
deduce  the  change  in  nuclear  isomeric  radii.  Un¬ 
fortunately,  the  enormous  difficulties  involved  in 
determining  theoretically  either  electron  or  nu¬ 
clear  densities  has  left  the  promised  potential 
largely  unfulfilled. 

In  the  face  of  the  greater  difficulties  inherent 
in  the  state  of  nuclear  theory,  some  progress  has 
been  made  in  this  field  by  assuming,  as  a  first 
crude  approximation,  that  the  electronic  charge 
difference  between  several  ions  in  different  chemi¬ 
cal  environments  may  be  ascribed  to  a  change  in 
tonicity.5**  Thus,  Hartree-Fock  calculations  have 
been  performed,  on  various  charge  states  of  an 
atom’  and  the  resultant  electron  charge  densities 


at  the  nucleus  are  then  used  to  estimate  relative 
isomer  shifts.  Implicit  in  this  procedure  is  the 
assumption  that  the  difference  in  nuclear  radii  is 
the  same  for  the  various  ionic  states.  Using  these 
calibrations,  comparison  with  experimental  iso¬ 
mer  shifts  then  provides  estimates  of  other  effects 
due  to  the  chemical  environment,  such  as  ligand 
bonding  and  electron  correlations. 

The  isomer  shift  SAB  between  a  source  B  and  an 
absorber  A  is  given  by  the  well-known  expres¬ 
sion*'9’7 

SAB-  (27c/3£y)Ze2[|<l'A(0)|2  -  j<ts(0)|*]fi(r*> , 

(1) 

where  £„  is  the  y-ray  energy,  Z  is  the  nuclear 
charge,  t’i d- A (0); 2  and  ej.s(0).z  are  the  electron 
charge  densities  at  the  nuclei  of  the  source  and 
the  absorber,  respectively,  and  i.  r2)  is  the  change 
in  mean-square  nuclear-charge  radius  between 
isomers.  For  a  given  pair  of  isomers,  <5vrJ>  is 
constant,  and  i c 4 (0 )[ 2  -  j^fl(0)|J  measures  the  change 
in  electronic  charge  densities  as  the  ions  are 
placed  in  different  compounds.  Implicit  in  this 
formulation  are  the  assumptions  that  the  elec¬ 
tron  density  is  constant  over  the  nuclear  volume 
and  that  the  change  in  electronic  state  {aa  to  „s) 
does  not  affect  the  nuclear  densities.  Thus,  the 
calculated  values  of  i>AB  for  various  charce  states 
of  the  same  isomers  yield  the  dependence  of  iso¬ 
mer  shift  on  changes  in  atomic  electron  density; 
the  differences  between  measured  and  calculated 
S4a  allow  one  to  deduce  the  additional  effects  of 
the  chemical  environment. 

This  paper  reports  results  of  a  study  of  the  rela- 
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ELECTRONIC  STRUCTURE. 
AND  NEUTRON  MAGNETIC 


The  magnetic  field  induced  magnetization  density  and  neutron  magnetic  form  factor  of  Pd  metal  is  obtained  from  an  ab 
initio  APW  energy  band  study  of  its  electronic  structure  and  properties.  The  magnetization  consists  of  the  spin  density  cal¬ 
culated  for  states  on  the  Fermi  surface  and  a  much  smaller  orbital  contribution.  The  solid  state  wavefunctions  are  found  to 
yield  a  spatial  localization  of  the  spin  density  which  is  greater  than  that  of  the  very  contracted  Hartree-Fock  density  of 
the  free  PdJ  *  ion.  The  theretical  magnetic  form  factor,  which  is  dominated  by  the  contribution  of  the  fifth  band,  is  found 
to  be  in  excellent  agreement  with  the  measurements  of  Cable.  Woiian.  Felcher.  Brun  and  Hoinfeldt. 


1.  Introduction 

The  transition  metais  have  been  studied  more  ex¬ 
tensively  than  have  other  metal  groups  in  the  periodic 
table.  Their  electronic  structure  and  properties  nave 
been  measured  using  a  variety  of  experimental  meth¬ 
ods  which  yield  both  macroscopic  information  (e.g. 
magnetization,  magnetic  susceptibility  and  specific 
heat)  and  microscopic  information  (e.g.  de  Haas-van 
Alphen,  magnetoresistance,  optical  reflectivity,  photo- 
emission)  which  have  been  related  directly  to  the  elec¬ 
tronic  band  structure.  By  contrast,  there  are  few  ex¬ 
periments  which  yield  direct  information  about  the 
nature  of  wavefunctions  in  these  metals.  The  determi¬ 
nation  of  the  neutron  magnetic  form  factor  in  magnet¬ 
ically  ordered  systems  which  gives  the  Fourier  trans¬ 
form  of  the  magnetization  density  is  perhaps  the  one 

t  Supported  by  the  National  Science  Foundation,  the  US  Air 
Force  Office  of  Scientific  Research  and  the  US  Energy  Re¬ 
search  and  Development  Administration  under  Contract 
Nos.  W-7405-eng-82  and  W-3l-109-eng-38. 


outstanding  exception.  While  theory  has  been  success¬ 
ful,  in  most  cases,  in  yielding  an  adequate  account  of 
the  energy  dependent  properties,  it  has  had  distinctly 
poorer  success  in  explaining  observed  form  factors  be¬ 
cause  of  the  greater  theoretical  (and  computational) 
problem  associated  with  treating  magnetically  ordered 
systems.  [In  ferromagneis.  one  must  include  the  ex¬ 
change  interactions  between  ail  Bloch  electrons  (which 
results  in  so-called  spin  polarized  bands)  which  are 
treated  only  approximately  by  spin  polarized  band 
calculations  and  obtain  the  magnetization  density  as 
the  difference  between  spin-up  and  spin-down  band 
magnetization  densities  each  obtained  as  sums  over  all 
occupied  (spin)  states.]  The  magnetic  form  factor 
studies  ( 1  — 3  J  have  found  that  for  the  3d  transition 
metai  ferromagnets  the  spin  densities  are  considerably 
contracted  spatially  and  are  well  represented  by  those 
for  the  Hartree— Fock  ions  in  the  +2  ionization  state 
[4],  By  contrast,  the  spin  densities  obtained  from 
band. structure  calculations  yield  form  factors  which 
are  not  in  as  good  agreement  with  experiment  (5-Sj. 
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The  role  of  matrix  elements  in  the  calculation  of  the  generalized  susceptibilities  y(q)  of  metals  is  discussed 
using  the  results  of  an  augmented-plane-wave  energy-band  calculation  for  the  eigenvalues  and  eigenfunctions 
of  Sc  metal.  The  inclusion  of  the  oscillator  strength  mains  elements  is  found  to  significantly  alter  the 
structure  obtained  for  yiqj  m  the  constam-matnx-elemem  approximation.  The  effect  of  local-field  corrections 
on  the  phonon  anomaly  in  Sc  and  the  observed  magnetic  ordering  of  dilute  rare-earth  Sc  alloys  are  described: 
in  the  latter  case,  a  crude  estimate  of  this  elfect  is  found  to  restore  a  (broad)  peak  in  y(q)  at  q  s  [0.0,0.5(ir/c)|, 
in  good  agreement  with  experiment. 


I.  INTRODUCTION 

The  generalized  susceptibility  function  which 
measures  the  response  of  the  system  to  an  ex- 
ternai  probe  is  central  to  the  understanding  of 
many  physical  phenomena  in  solids.  .Although  the 
calculations  of  these  response  functions  are 
rather  involved,  the  considerable  refinement  in 
the  computational  techniques  which  has  been 
achieved  in  recent  years  makes  accurate  ab  irtitio 
studies  possible.  The  greater  realization  that 
(except  perhaps  for  very  simple  metals)  the  true 
response  function  of  a  real  solid  is  quite  unlike 
the  frequently  used  free-electron  Lmdhard  func¬ 
tion  has  led  to  attempts  to  include  the  real  band- 
structure  effects1'2  into  the  calculations.  In  the 
random-phase  approximation  (RPA),  the  dynamic, 
i.e.,  wave-vector  and  frequency  dependent,  elec¬ 
tron-electron  dielectric  function  is  written 


c (q,  u>)  =  1  +  (4r e2/fty2)  x(q,  w)  , 

(1) 

where  $2  is  the  crystal  volume  and 

x(q.  <*>) =Xi(5.  w) +ix*( q,  <a) 

(2) 

(3) 

Here  £,  is  the  Bloch  energy  of  the  state  k  with 
Bloch  function  &k  and  occupation  number  nk;  k  in¬ 
cludes  both  the  band  index  and  the  wave  vector 
k  which  is  restricted  to  the  first  Bnllouin  zone. 
The  oscillator-strength  matrix  element 

=  ,  (4) 

imposes  the  restriction  k'  =  £  +  q  -i-G,  where  G  is  a 
reciprocal-lattice  vector;  w  denotes  the  frequency 
of  the  external  field,  and  \(q,  w)  is  the  bare  sus¬ 


ceptibility  of  the  solid,  i.  e. ,  the  effects  due  to 
local  fields  are  not  included.  The  zeros  of  e-(q,a') 
=  Re[€(q,  j)J  determine  the  piasmon  dispersion  re¬ 
lation  and  €2(q,  w)  =  Ir.gs iq,  -a)j  is  related  to  many 
properties  including  ii'.e  dynamic  form  factor 
S(q,  j),  the  optical  conductivity,  and  the  optical 
reflectivity.1 

The  static  response  function  x(5)>  i.e.,  the  re¬ 
sponse  function  for  a>~0  is  itself  of  very  great 
interest  because  of  its  possible  role  in  causing 
electronically  driven  instabilities  which  are  mani¬ 
fested  as  spin-density  waves,  charge-density 
waves,  or  in  structural  phase  transformations. 
From_£q.  (3)  we  can  write  down  the  expression 
for  x(qh 

x(q)  =  x(q.  a>  =  o) 


!AWq)l2 


kk'  ^  k  ~~ 


(5) 


It  has  been  argued  that  due  to  the  ‘‘nesting"  fea¬ 
tures,  i.  e. ,  the  existence  of  flat  and  parallel 
areas  of  the  Fermi  surface,  the  denominator  in 
Eq.  (5),  will  tend  to  be  vanishingly  small,  giving 
rise  to  a  large  peak  in  \(q)  at  the  nesting  wave 
vector.  Such  a  large  peak  in  y(q)  coupled  with 
exchange  enhancement  effects  in,  for  instance, 
paramagnetic  chromium,4’5  has  been  thought  to 
be  the  main  reason  lor  driving  the  paramagnetic 
phase  unstable  with  respect  to  an  antiferromag¬ 
netic  phase,  and  has  since  been  used  by  a  number 
of  workers  to  account  for  various  other  instabil¬ 
ities  related  to  Fertm-surface  geometry.  The 
occurrence  of  anomalies  in  phonon  spectra  of 
metals  has  also  been  related  in  a  similar  way  to 
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Molecular-cluster  modeU  are  developed  to  describe  ground-state  properties  of  EuO  and  EuS  in  the  Hartrte- 
Fock-Slater  one-electron  approximation.  Spin-polanz-d  (Eu.Y,)lJ"  complexes  are  examined  using  both  neutral- 
lad  ionic-model  potentials  which  incorporate  a  part  of  the  effect  of  the  crystalline  environment.  Self-consistent 
calculations  are  made  for  (EuO„)'0_.  From  the  charge  and  spin  densities,  the  transferred  hyperfine  held  at  the 
O  site  in  EuO  is  found  to  be  —  8  ±  2  kG  and  a  small  solid-state  bonding  effect  is  predicted  for  the  neutron 
magnetic  form  factor.  The  pressure  dependence  of  "the  charge  density  at  the  Eu  nucleus  in  EuS  is  determined 
as  a  function  of  bond  length  in  the  (EuS,)'0-  cluster  and  used  to  obtain  from  the  experimental  data  an  isomer- 
shift  calibration  constant  a  =  — 0.48a(  mm/sec.  The  one-electron  energy  leveis  of  the  (EuO.)1’*  cluster  are 
found  to  be  in  good  agreement  with  the  augmented-plane-wase  results  of  Cho  when  both  calculations  are 
performed  with  similar  model  potentials.  The  extension  to  selt-consistency  leacs  to  significant  energy -level 
rearrangement  which  indicates  the  importance  of  final-state  relaxation  and  Coulomb  correlation  effects  in  the 
interpretation  of  experimental  spectra. 


L  INTRODUCTION 

The  europium  chalcogenides  Eu.Y  (.Y  =  0.  S,  3ej  and 
Te),  an  important  class  of  magnetic  semiconduc¬ 
tors,  have  magnetic  ordering  properties  which 
vary  from  ferromagnetic  (EuO)  to  antiferromag¬ 
netic  (EuTe),  and  exhibit  a  red  shift  in  their  opti¬ 
cal  spectra  with  increasing  magnetic  order.  In¬ 
sulators  in  the  pure  state,  they  become  semicon¬ 
ducting  when  doped  with  rare-earth  impurities, 
and  show  large  negative  magnetoresistances  near 
their  transition  temperatures.  A  review  of  the 
experimental  properties  and  of  the  theoretical  in¬ 
terpretations  has  been  given  by  Methfessel  and 
Mattis.1  The  general  picture  which  has  emerged 
is  one  of  magnetic  moments  localized  on  Eu2'  ions 
with  a  high- spin  configuration  (4/i)7(4/f)°.  These 
moments  are  couplecLby  an  antiferromagneiic  130° 
superexchange  interaction,  acting  through  the  X:~ 
ligand  ions,  and  by  a  ferromagnetic  exchange 
which  depends  on  the  covalent  mixing  of  Eu  4/  or¬ 
bitals  with  the  unoccupied  5 d  orbitals  on  nearest- 
neighbor  Eu  sites.  Now  most  rare  earths  are  tri- 
valent  yielding,  in  their  respective  monochalco- 
genides,  two  electrons  to  the  ligand  ion,  and  one 
to  the  conduction  band,  so  that  the  compounds  are 
metallic.  Added  as  impurities  to  the  insulating 
Eu  chalcogenides,  rare-earth  ions  populate  the 
conduction  band  with  a  fraction  of  their  third  val¬ 
ence  electrons.  Since  the  conduction  band  is  in¬ 
volved  in  magnetic  ordering  (as  shown  by  the  mag¬ 
neto-optical  red  shift,  and  as  expected  from  the 
nature  of  the  ferromagnetic  interaction),  the  phe¬ 
nomena  associated  with  the  magnetic  phase  tran¬ 
sition  are  modified  and  the  negative  magnetore¬ 
sistance  effects  are  obtained.  More  detailed  eluci¬ 
dation  of  these  effects  has  been  attempted  in  terms 


of  both  impurity  and  energy -band  models,1  but  no 
complete  theory  can  be  said  to  have  been  found. 

The  spin-polarized  one-electron  theory  has  been 
quite  successfully  applied  in  interpreting  magnetic 
properties,  optical  data,  and  photoemission  spec¬ 
tra  of  lighter  materials,  notably  the  2d  metals  and 
their  compounds.  While  the  limitations  of  such 
a  model  are  evident,  for  example,  in  their  in¬ 
ability  to  predict  correctly  multiplet  splittings  ob¬ 
served,  it  is  clearly  worthwhile  to  investigate 
their  applicability  to  open- shell  rare-earth  sys¬ 
tems. 

This  paper  describes  the  results  of  theoretical 
molecular-cluster  calculations  for  EuO  and  EuS 
in  the  Hartree-Fock-Slater  (HFS)  approximation. 
Both  model-potential  and  self-consistent  studies 
are  undertaken  to  better  understand  the  role  of 
"localized”  4 f  electrons  in  the  bonding  and'mag- 
netic  interactions  of  rare-earth  compounds.  Sev¬ 
eral  ground-state  properties  of  these  systems, 
notably  charge  and  spin  densities,  are  used  to  de¬ 
termine  such  experimental  properties  as  neutron 
magnetic  form  factors  and  transferred  hyperfine 
fields.  Results  for  the  pressure  dependence  of 
the  electron  density  at  the  Eu  nucleus  in  EuS  are 
used  in  conjunction  with  the  experimental  Moss- 
bauer  isomer-shift  data  of  Kalvius  et  al }  to  de¬ 
termine  a.  the  isomer-shift  calibration  constant, 
for  151Eu. 

As  is  well  known,  parametrized  band -structure 
models3, 1  have  been  presented  which  are  compati¬ 
ble  with  a  one-electron  interpretation  of  optical, 
x-ray,  photoemission,  and  magnetic  data  for  the 
europium  chalcogenides.  Thus,  the  applicability 
of  the  one-electron  HFS  model  to  the  description 
of  some  of  these  excited-state  properties  is  also 
tested  in  this  work,  with  special  attention  to  4/ 
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Abstract 

This  novel  approach  combines  a  discrete  variational  treatment  of  all  potential  ta  ms  arising  from  the  su¬ 
perposition  of  the  spherical  overlapping  atomic  ^r.arge  densities  with  a  rapid iv  convergent  Fourier  series  rep¬ 
resentation  of  all  multicenter  nonspnencal  potential  terms.  The  basis  set  consists  of  me  exact  numerical  atomic 
valence  orbitals,  augmented  by  charge  transier  states,  virtual  atomic  states,  and  single  ar.aivtic  Slater  orbitals 
for  increased  variational  flexibility.  The  initial  potential  is  a  non-mulfin-iin  overlapping  atomic  potential  in¬ 
cluding  nongradient  local  density  exchange  and  correlation  terms.  Full  self-consistency  is  obtained  by  a  pro¬ 
cedure  that  combines  an  iterative  scheme  within  the  superposition  mode!  wuh  a  self-consistent  optimization 
of  the  Fourier  components  of  the  nonspherieal  cnarge  density  terms.  Urounu-siate  prope,  ties  sucn  as  structure 
factors  and  cohesive  energy  are  computed.  The  results  for  diamond  show  very  good  agreement  with  experiment. 
Comparison  of  the  results  with  the  Hartree-Fock  calculation  is  discussed. 

1.  Introduction 

It  is  now  widely  recognized  that  energy  band  theory  has  become  a  powerful  and  so¬ 
phisticated  method  for  studying  a  wide  spectrum  of  solid-state  properties.  The  prolifer¬ 
ation  of  energy  band  schemes  and  their  application  to  increasingly  diverse  problems  attests 
to  the  current  popularity  of  band  theory .  Increasingly  too.  however,  its  very  success  in 
describing  a  host  of  sophisticated  new  solid-state  experiments  has  led  to  a  “tide  of  rising 
expectations"  which  the  various  existing  computational  energy  band  schemes  (such  as 
APW.  OPW.  KKR.  etc.)  have  been  hard  pressed  to  satisfy.  In  addition  to  challenging  ex¬ 
periments  performed  on  important  materials  having  complex  crystallographic  structures, 
these  new  experiments  have  demanded  not  only  theoretical  descriptions  of  high-resolution 
eigenvalue  phenomena,  but  also  detailed  and  precise  solid-state  wave  functions  with  which 
to  determine  the  expectation  values  of  different  observable  operators.  Such  a  demanding 
test  of  the  predictions  of  one-electron  theory  has  the  additional  virtue  in  permitting,  by 
their  comparison  with  experiment,  accurate  determinations  of  the  relative  magnitude 
and  importance  of  many-bodv  effects  in  real  solids. 

In  this  paper  we  describe  a  new  approach  to  the  fully  self-consistent  solution  of  the 
one-particle  equations  in  a  periodic  solid  within  the  focal  density  functional  formalism 
[IJ.  It  is  specifically  designed  and  developed  to  incorporate  special  features  with  which 
to  overcome  difficulties  encountered  by  other  methods.  Specifically,  as  will  be  shown 
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The  soft— x-ray  L ..  emission  and  absorption  spectra  of  Mg  metal  have  been  determined 
by  means  of  high-resolution  ah  initio  aucmented-plane-wave  calculations.  The  inclusion 
of  transition-matrix  elements  calculated  from  the  augmented-plane-wave  wave  functions 
changes  the  sharp  structure  found  in  the  joint  density  of  states,  yields  an  L -  .  emission 
spectrum  in  good  agreement  with  experiment,  and  indicates  that  band-structure  effects 
play  an  Important  role  in  the  observed  x-ray  threshold  shape. 

Despite  both  intensive  and  extensive  efforts,  threshold  (e.g.,  the  L,j  edges  of  Na,  Mg,  and  Al) 

the  theoretical  understanding  of  x-ray  emission  or  broad  and  rounded  (e.g.,  the  K  edges  of  Li. 

and  absorption  edges  in  metals  is  still  a  subject  Be.  Mg,  and  AO  over  energies  =>  1  eV.  Since 

of  controversy  and  disagreement  mainly  because  these  features  and  especially  the  sharply  peaked 

of  the  differing  roles  assigned  to  band- structure  nature  of  the  edges  do  not  seem  explicable  in 
(one-electron)  and  conduction-electron-core-  a  one-electron  theory  based  on  the  free-electron 

hole,  final-state  interaction  (many-bodv)  effects. 1,2  model,  many-body  models  were  developed  by 
In  the  simple  metals  Li,  Be,  Na,  Mg,  and  A1  the  Mahan3  and  NoziSres  and  de  Dominicis4  in  which 

shapes  of  the  observed  threshold  edges  in  both  the  conduction-electron-core-hole  interactions 

emission  and  absorption  fail  to  exhibit  the  sharp  play  a  central  role.  These  theories  predict 
step-function  edges  expected  to  be  characteristic  threshold  edges  whose  shapes  are  characterized 

of  free-electron  metals.  Instead,  experimentally  by  threshold  exponents  a9  (for  edges)  and  a, 

determined  edges  are  either  sharp  and  peaked  at  (for  K  edges);  for  photoemission  the  asymmetry 
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The  generalized  susceptibility  of  NbC,  determined  from  its  augmented  plane  wave  en¬ 
ergy  band  structure  (with  constant  matrix  elements',  shows  large  maxima,  arising  in 
large  part  from  Fermi  surface  “nesting,”  to  occur  at  precisely  those  q  values  at  which 
soft  inodes  have  been  found  in  the  phonon  acoustic  o  ran  ones  by  Smith.  These  maxima 
can  be  represented  by  a  warped  cube  of  dimension  ~  1 .2 1 2 -  a)  in  momentum  space,  in 
striking  agreement  with  the  soft  mode  surface  proposed  by  A'eber. 


Long  known  for  their  remarkable  properties 
such  as  high  melting  points  ai.d  ultrahardness, 
the  transition- metal  carbides  and  nitrides  are 
considered  an  important  class  of  materials  be¬ 
cause  of  the  occurrence  of  high  superconducting 
transition  temperatures,  Tc.  The  systematic 
neutron-scattering  determinations  of  phonon  dis¬ 
persion  curves  w(q)  in  the  carbides  of  the  transi¬ 
tion-metal  IV  and  V  series  by  Smith  and  Glaser1 
revealed  the  existence  of  soft  phonons  and  led  to 
attempts  to  correlate  these  soft-phonon  modes  to 
the  high  Te  values.  Such  a  direct  correlation  be¬ 
tween  lattice  instabilities  and  high  Tc’s  is  becom¬ 
ing  more  accepted;  an  important  problem  is  the 
theoretical  understanding  of  the  origin  of  phonon¬ 
mode  softening  in  superconductors.  While  rigor¬ 
ous  formulations  of  dynamical  matrix  theory  have 
been  proposed,  such  as  the  work  of  Sinha,5  based 
on  the  augmented  plane-wave  (APW)  method,  no 
ab  initio  calculational  predictions  of  phonon  spec¬ 
tra  have  been  reported.  A  phenomenological 
model  theory  utilizing  the  modified  shell  and  su¬ 
pershell  model  has  been  employed,  particularly 


for  NbC,  by  Weber,  Bilz,  and  Schroder.5  Weber4 
was  able  to  show  that  the  introduction  of  the  su¬ 
pershell  interactions,  which  describe  the  d-d 
metal  interaction,  is  essential  for  obtaining  the 
observed  dips  in  ouiq).  An  important  result  was 
obtained  oy  Mostoller5  when  he  showed  that  the 
general  shapes  of  the  dispersion  curves,  but  not 
the  phonon  anomalies,  could  be  obtained  using  a 
free-electron  dielectric  screening  function  and 
describing  the  electron-ion  interaction  in  a  pseu¬ 
dopotential  representation.  More  recently,  Sin¬ 
ha  and  Harmon’  were  able  to  reproduce  the  longi¬ 
tudinal  branches  of  the  phonon  dispersion  curves 
in  Nb  and  NbC  by  using  a  local  pseudopotential 
to  represent  the  electron-ion  interaction.  In  this 
work,  which  necessarily  introduced  other  approx¬ 
imations,  the  dielectric  screening  matrix  was 
expressed  in  terms  of  several  parameters  which 
were  determined  by  fitting  the  calculated  longi¬ 
tudinal  modes  to  the  experimental  values.  While 
noting  the  possible  importance  of  structure  in  the 
diagonal  d-band  part  of  the  susceptibility  func¬ 
tion,  Sinha  and  Harmon5  assumed  *(q)  to  be  q  in- 
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Large  maxima  in  the  generalised  ■.usceptibilitv  of  x(<7>  of  TaC  calculated  from  is  APW  band  structure  die  found 
to  correlate  very  well  with  observed  suit  phonon  modes.  Such  maxima  are  not  found  for  HfC  where  ononon  anoma¬ 
lies  ate  not  observed. 


Systematic  neutron  scattering  determinations  of 
phonon  dispersion  curves  cc( q)  in  the  caroides  of  the 
transition  metal  IV  and  V  series  by  Smith  and  Glaser 
[  I  ]  revealed  the  existence  of  soft  phonons  and  led  to 
attempts  to  correlate  these  phonon  anomalies  to  their 
observed  high  superconducting  transition  temperatures. 
Tc.  By  contrast,  phonon  anomalies  are  not  observed 
in  low  Tc  (<0.05  R)  compounds  such  as  HfC  and  ZrC. 
Theoretical  efforts  to  understand  phonon  dispersion 
curves  of  the  transition  metal  carbides  have  been  based 
either  qn  phenomenological  models  [2-4|  or  para¬ 
meterized  theories  [S,6|  asub  initio  calculations  of 
the  dynamical  matrix  within  a  microscopic  formation 
are  still  not  available.  We  have  recently  reported  p] 
that  a  detailed  ab  initio  study  of  the  generalized  sus¬ 
ceptibility,  \(q),  of  NbC,  determined  trom  us  aug¬ 
mented  plane  wave  (  APW)  energy  band  structure, 
shows  large  maxima  arising  (in  part  I  from  Fermi  sur¬ 
face  “nesting”,  to  occur  at  precisely  those  q  values  at 
which  minima  (soft  modes)  have  been  found  in  the 
phonon  acoustic  branches  by  Smith  [  1  ]  .  These  maxi¬ 
ma  can  be  represented  by  a  warped  cube  of  dimension 
~1.2(2rr/<z)  in  momentum  space,  in  striking  agreement 

0  Supported  by  An  Force  Office  ot  Scientific  Research.  Na¬ 
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with  the  soft  mode  surface,  proposed  phenomenologi¬ 
cally  by  Weber  [3] . 

In  this  letter,  we  report  results  of  similar  investiga¬ 
tions  of  TaC  w  hich  is  isoelectronic  to  NbC.  has  a  T  = 
10.35  R  ivs  1  l.i  R  for  NbC)  and  also  shows  strong 
phonon  anomalies.  As  in  NbC  we  find  large  maxima 
in  \( t7 j  at  precisely  those  q  vaiues  at  winch  the  soft 
modes  have  been  observed  [  1 ) .  By  contrast,  the  ca!c» 
lated  \(<7 )  for  HfC  i  based  on  a  rigid  band  shift  of  c\ 
in  the  TaC  band  structure)  shows  no  maxima  and 
hence  no  phonon  softening  -  again  in  agreement  with 
experiment. 

— .  It  has  long  been  known  that  the  screening  of  a  oho- 
non  mode  with  wave  vector  q  can  be  strongly  enhanced 
by  a  iarge  scattering  of  the  conduction  electrons 
across  "nesting”  portions  of  the  Fermi  surface,  spanned 
by  the  vector  q.  Thus,  the  generalized  function  \tq ) 
plavs  an  important  role  in  the  lattice  dynamics  of  the 
system,  and.  under  certain  simplifying  assumptions, 
it  can  be  shown  that  if  \  becomes  large,  the  renorma¬ 
lized  phonon  frequencies  (due  to  electron-phonon  in¬ 
teraction)  can  become  small  and  consequently,  sort 
phonon  modes  (or  lattice  instabilities)  can  be  expec¬ 
ted.  In  the  constant  matrix  element  approximation. 
\(q)  is  defined  as: 

x«7)  =  -  S  (/?- )/(££ -£*%)• 
n.  n 
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ABSTRACT 

*  Large  maxima  in  the  generalized  susceptibility,  x  (  q )  of 
NbC  and  TaC  calculated  from  their  APW  band  structures  are 
found" to  correlate  very  well  with  soft  phonon  modes  observed 
hy  Smith.  Such  maxima  are  not  found  for  TrC  and  HfC  where 
phonon  anomalies  are  not  observed. 
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ABSTRACT 

Whereas  a  variety  of  experimental  methods  yield  information  which 
has  been  related  directly  to  the  electronic  band  s true ture ,  there  are 
few  experiments  which  yield  direct  information  about  the  nature  of  wave- 
functions  in  metals.  Perhaps  the  one  outstanding  exception  is  the  deter¬ 
mination  of  the  neutron  magnetic  form  factor  which  gives  the  Fourier 
transform  of  the  magnetization  density.  In  the  past  few  years,  accurate 
neutron  experiments  have  shown  the  inadequacy  of  the  commonly  used  free 
arom/ion  Hartree-Fock  form  factors.  This  paper  describes  recent  advances 
in  high  resolution  energy  band  studies  of  APW  wave f ur.c ’"ions ,  magnetiza¬ 
tion  densities,  and  neutron  form  factors  in  metals.  Illustrative  results 
obtained  in  our  Laboratory  for  several  important  transition  and  rare-earth 
metal  systems  are  presented  and  corrmared  with  experiments  carried  out  at 
ORNL,  ANL,  and  Ames  Laboratory.  Unlike  the  case  of  the  ferromagnetic 
metals,  we  find  very  good  agreement  with  experiment  for  the  paramagnetic 
mecals  studied  to  date. 
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Marked  structure  found  in  an  APW  study  of  both  the  density  of  states  and  .ransition  matrix  elements  contribut¬ 
ing  to  the  emission,  absorption  spectra  of  Na  emphasizes  the  need  for  including  band  structure  effects  in  all  X-ray 

edge  analyses. 


i 


Sodium  metal  has  been  widely  assumed  to  be  the 
prototype  of  a  free-electron  metal,  because  its  Li  3 
emission  spectrum  [lj  is  free-eiectron-iike.  its  Fermi 
surface  is  spherical,  its  effective  mass  is  nearly  equal 
to  the  free-electron  mass,  and  its  optical  properties 
are  well  described  by  a  Drude  model.  Here  we  report 
augmented  plane  wave  (APW)  band  structure  calcula¬ 
tions  which  are  consistent  with  the  above-mentioned 
free  electron  properties  of  Na,  but  which  produce  sig¬ 
nificant  and  qualitative  deviations  from  the  free  elec¬ 
tron  picture  of  sodium's  L-.  3  absorption  edge.  Our 
calculations  reveal  marked  structure  in  both  the  con¬ 
duction  band  density  of  states  and  in  the  core-conduc¬ 
tion  band  dipole  matrix  elements;  they  indicate  that 
at  least  part  of  what  had  once  been  thought  to  be  a 
many-body  threshold  anomalv  at  the  Li  -<  X-ray  ab¬ 
sorption  edge  of  Na  may  be  a  band  structure  effect; 
and  they  call  into  question  ail  X-rav  edge  analyses 
which  do  not  account  for  band  structure  effects. 

In  the  one-electron  approximation  the  intensity 
1(E)  for  the  emission  or  absorption  of  X-rays  of  ener¬ 
gy  £  is  given  by 


I(F)=c  fdkt<^cl7l^k)l25(E  -£*  +  £c) 


f  dsk 


i<cg7iv*>i:. 


(D 
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Here  Ck  and  Ek  jre.  respectively,  the  Bloch  wavefunc- 
non  and  rhe  Bloch  energy  of  an  electron  with  wave 
vector  k  and  band  index  n.  uc  and  E:  are  the  core 
function  and  the  core  energy'  of  the  optical  electron 
(2p  state  for  the  Li  1  spectrum  of  Na).  and  c  is  a  con¬ 
stant  of  proportionality .  The  integral  in  eq.  ( 1 )  is  over 
a  constant  energy  surface  E  =  Ek  -  Ec  in  k- space.  For 
the  emission  spectrum,  all  the  energy  states  Ek  are  oc¬ 
cupied  and  for  the  absorption  spectrum,  the  states  are 
unoccupied. 

Note  from  eq.  ( 1 )  that  in  the  constant  transition 
matrix  element  (TMF.)  approximation,  the  intensity 
depends  on  the  partial  density  of  states  I  DOS)  because 
of  the  dipole  selection  rule  and  not  on  the  total  DOS. 
in  fig.  1  we  show  the  partial  DOS  for  Na  obtained 
from  our  APW  calculation.  The  free-electron  DOS  is 
shown  in  fig.  2  for  comparison.  We  notice  that  the 
APW  total  DOS  and  the  partial  DOS  are  significantly 
different  from  their  counterparts  in  the  free-eieetron 
model,  especially  at  energies  above  the  Fermi  enerey, 
£p .  In  particular,  we  find  structure  in  the  APW  DOS 
which  is  absent  in  the  free-electron  DOS.  The  posi¬ 
tions  of  the  peaks  tr.  our  s  and  d  DOS  correlate  very- 
well  with  those  observed  in  the  Li  3  absorption  spec¬ 
tra. 

From  figs,  i  and  2,  it  :s  clear  that  only  in  the  free- 
electron  approximation  and  neglecting  dipoie  selec¬ 
tion  rules  completely,  can  one  say  that  the  X-ray  in¬ 
tensity  ought  to  be  parabolic.  One  notices  that  the  in¬ 
clusion  of  dipole  selection  rules,  even  for  a  free-elec- 
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APW  band  calculations  ol  the  suit  X-rav  K  absorption  and  emission  'oectra  ol  Li  are  used  to  discuss  the  extent  to 
which  various  conflicting  theories  and  measurements  of  absorption,  emission,  electron-energv  -loss,  and  photoemission 

can  be  reconciled 


The  K  X-ray  emission  and  absorption  edges  of  Lt. 
although  intensely  studied  for  almost  half  a  century, 
remain  the  subject  of  both  experimental  and  theoreti¬ 
cal  controversy.  Experimentally,  a  concensus  has  yet 
to  be  reach  on  (i)  whether  the  high  energy  emission 
threshold  and  the  low  energy  absorption  edge  are  very 
broad  [1-4]  (~0.4  eVj.  broad  [4.  5]  ( —  0.2  eVj,  or 
sharp (~0. 05  eV)  [6] : (ii)  whether  [7]  or  not  [3] 
absorption  and  emission  edges  overlap.  ( ml  whether 
[2,  8]  or  not  (3,  9]  the  edge  breadths  exhibit  more 
temperature  dependence  than  a  Fernn  factor,  (iv) 
whether  the  maximum  of  the  emission  spectrum  is  at 
threshold  or  eV  below  (4,  5]  ( v)  whether  [6]  or 
not  [3j  the  absorption  spectrum  has  a  relative  maxi¬ 
mum  eV  above  threshold,  (vi)  the  importance  ot 
impurities  in  shaping  the  edges  [5] .  and  (vii)  the  quan¬ 
titative  shape  of  X-ray  photoemission  imes  [8,  10). 
Theoretically,  there  are  still  debates  on  the  relative 
importances  of  (i)  band  structure  effects  (11.  12) 

(ii)  Auger  broadening,  [  1 3 j ,  (hi)  phonon  broadening. 
[11.  14]  ,  (iv)  electron-hole  scattering  resonances 
[16,  17]  and  (v)  the  Nczieres  -de  Dominicis  many- 
body  effect  [18].  in  this  paper  we  report  augmented 
plane  wave  (APW)  band  structure  calculations  for  Li. 
and  indicate  the  extent  to  which  they  constrain  the 

•  Supported  by  the  NSF  (Grants  DMR  72-03019.  72-03101 
A02,  73-07661  and  72-03026)  and  the  Air  Force  Office  of 
Scientific  Research  (Grant  76-2948). 


possible  interpretations  of  the  K  edge  spectra  of  Li. 
Since  similar  calculations  of  the  L-,  3  X-ray  spectra  of 
Mg  successfully  described  the  spectral  positions  of 
all  the  observed  structures,  and  produced  spectral 
shapes  in  excellent  semi-quantitauve  agreement  with 
the  data  [19],  we  have  considerable  confidence  in 
the  method's  ability  to  predict  semi-quantitatively 
the  shapes  of  the  K.  emission  and  absorption  spectra 
of  Li. 

The  X-ray  emission  or  absorption  intensity  1(E)  is 
given  in  the  one-electron  approximation  by 

SjV*.-v,.T  :5 t£-E„k+£.)  (1) 

where  Urk  and  Enk  are  the  Bloch  function  and  energy 
respectively  of  a  conduction  electron  in  band  n  with 
wavevectorA:  2C  is  a  Is  Hermcn-Skillman  atomic  core 
function  of  energy  £'c,  and  Enk  is  restricted  to  occu¬ 
pied  (empty)  conduction  stares  foremission  (absorp¬ 
tion).  We  have  evaluated  eq.  ( 1)  using  eigenvalues  and 
wave  functions  found  at  285  non-equivalent  points  A 
in  the  ii  reducible  1 ,48th  portion  of  the  Bnllouin  zone, 
using  the  accurate  tetrahedron  method  [20.  21  ].  This 
method  gives  us  an  advantage  over  earlier  calculations 
which  were  limited  by  histogram  techniques.  The 
warped-muffin-tin  version  ot  the  APW  method  was 
employed,  permitting  the  potential  to  vary  in  the  in¬ 
terstitial  region.  The  potential  was  computed  assuming 
a  lattice  constant  =  0.5971  au,  a  superposition  of 
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The  possible  role  of  electronic  structure  on  observed  phonon  anomalies  in  high-temperature  superconducting 
transition-metal  carbides  is  studied  by  means  of  accurate  ab  initio  calculations  of  the  conduction  electron 
response  function.  From  augmented-plane-wave  determinations  of  the  electronic  band  structure,  density  of 
states,  and  Fermi  surface  of  NbC  and  TaC.  the  vsase-sector-dependent  generalised  susceptibility.  \tqt,  is 
calculated  in  the  constant-matnx-eiement  approximation.  For  noth  NbC  and  TaC.  \ufl  has  strong  maxima  at 
precisely  those  q  values  at  which  soft  modes  were  observed  by  Smith  ano  Glaser  Maxima  in  apt q i  are 
predicted  for  other  directions.  The  locus  of  these  q,.j,  values  can  be  represented  by  a  warped  cube  of 
dimension  —  I2(27r/a)  tn  momentum  space — in  striking  aereement  with  the  sot  (-mode  surface  proposed 
phenomenologically  by  Weber.  In  sharp  contrast,  the  \ i q i  calculated  for  noth  ZrC  and  HfC — for  which  no 
phonon  anomalies  have  been  observed — fall  off  in  all  s>  ametry  directions  awav  from  the  zone  center.  In 
agreement  with  Phillips,  we  thus  interpret' the  phonon  anomalies  in  the  transiuon-meul  carbides  as  due  to  an 
“overscreening"  effect  resulting  from  an  anomalous  increase  of  lhe  response  function  of  the  conduction 
electrons. 


I.  INTRODUCTION 

The  transition -metal  carbides  and  nitrides  of 
the  transition-metal  IV,  V,  and  VI  series  owe 
their  wide  commercial  use  to  their  remarkable 
physical  properties,  essentially  their  great  hard¬ 
ness  and  high  melting  points— properties  usually 
related  to  the  existence  of  covalent  bonds  betw  een 
the  metal  and  ligand  atoms.  Aside  from  their  in¬ 
teresting  physical  properties,  the  transition- 
metal  carbides  (and  nitrides)  under  study  are  good 
metallic  conductors  with  conductivities  comparable 
to  those  of  the  parent  transition  metal.  The  9  and 
10  valence-electron  compounds  are  superconduc¬ 
tors  with  high  superconducting  transition  tempera¬ 
ture  Tc  ( Tc=  11.1  K  for  NbC,  10.35  K  for  TaC,  and 
16.8  K  for  NbN),  while  the  8  valence -electron  com¬ 
pounds  like  ZrC  and  HfC  have  a  low  Te£  0.05  K. 

An  empirical  correlation  between  the  occurrence 
of  high  Tc  and  the  presence  of  phonon  softening  in 
the  phonon  spectra  of  the  transition-metal  carbides 
(TMC)  was  pointed  out  by  Smith  and  GlUser,1  and 
Smith2  after  a  systematic  study  of  the  phonon  dis¬ 
persion  curves  of  the  TMC  by  coherent  inelastic 
neutron  scattering.  Such  a  correlation  between 
phonon  anomalies  and  high  values  of  Tc  has  been 
also  observed  for  a  large  number  of  d-electron 
superconducting  compounds  or  alleys,  e.g.,  the 
well-known  Kb -Mo  system  for  which  Powell  et  al.3 
showed  that  the  low-7"c  alloy  Mb,  ,5Mo1  r5  (Tc 
=>0.04  K)  no  longer  possesses  any  phonon  anomaly. 
The  A-15  superconductors  like  VjSi  and  Nb.Sn  also 
display  a  correlation  between  high  values  of  T c  and 


the  occurrence  of  a  low- -temperature  structural 
phase  transformation'  (from  the  cubic  to  tetragonal 
phase),  or  the  existence  of  anomalous  mode  soften¬ 
ing  for  the  nontransforming  compounds.  However, 
even  though  instability  or  anomalous  phonon  soften¬ 
ing  in  many  systems  appear  to  be  correlated  with 
the  occurrence  of  a  high  Tc,  there  also  exist  some 
counter  examples.’  Quantitative  relations  between 
these  experimental  facts  are  more  difficult  to 
draw  from  first-principles  calculations.  However, 
in  the  existing  model  theories,  such  as  the  strong - 
coupling  theory,5  the  occurrence  of  soft  phonon 
modes  reflected  in  the  variation  of  the  second 
moment  of  the  phonon  frequency,  vai')  (defined  in 
Ref.  6),  does  not  appear  to  be  sufficient  to  explain 
the  variation  of  the  electron-phonon  coupling  con¬ 
stant  \  which  essentially  determines  T c.'"' 

Since  the  two  physical  phenomena,  phonon  mode 
softening  and  high  Tc,  occur  concomittantly  in  the 
TMC,  it  is  very  important  to  understand  the  origin 
of  the  phonon  anomalies  in  these  materials.  In 
this  work,  we  examtne  the  role  played  by  the  re¬ 
sponse  function  of  the  conduction  electrons  in  the 
occurrence  of  soft  phonon  modes,  namely,  we  cal¬ 
culated  from  first  principles  the  wave-number- 
dependent  generalized  susceptibility  function  y(q) 
in  the  constant-matrix-element  approximation 
for  the  9  and  8  valence-electron  TMC:  NbC,  TaC, 
and  ZrC,  HfC.  We  find  an  anomalous  increase  in 
the  susceptibility  of  the  8  valence -electron  TMC 
at  exactly  the  q  values  where  the  anomalies  in  the 
phonon  spectrum  are  observed  experimentally. 

By  studying  the  response  function  along  symmetry 
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The  generalized  susceptibility,  ;<(q) ,  of  both  NbC  and  TaC 
'determined  from  APW  energy  band  calculations  show  large  maxima  to 
occur  at  precisely  those  C(max  values  at  which  soft  phonon  modes 
were  observed  by  Smith.  Maxima  in  <(q)  are  predicted  for  other 
directions.  The  locus  of  these  qmax  values  can  be  represented  by 
a  warped  cube  of  dimension  n.1 . 2(2~/a)  in  momentum  cpace  -  in  striking 
agreement  with  the  soft  mode  surface  proposed  phenomenologically 
by  Weber.  In  sharp  contrast,  the  XlR)  calculated  for  both  ZrC 
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Abstract 

Prediction*  of  the  locai  donuts  formalism  approach  to  concsiso  properties  of  cosalentls  bonded 
solids  (diamond  and  cubic  BN  t  and  promts  pe  ionic  >;> stern  i  L1K1  ere  aescribed  usina  results  of  our 
recently  developed  fulls  self-consistent  numerical  basis  set  cc  su  ns  st  approach.  Comparisons 
with  restricted  Hartree-Focis  results  ana  experiment  lor  cor.esise  energies  anu  equilibrium  lattice 
constants  are  presented.  Some  of  the  principal  bonding  mechanisms  n  these  crystals  are  discussed 
in  terms  of  the  contributions  of  locai  e.xonunac  ana  correlation  to  the  binding  and  the  charge  re¬ 
distribution  relative  to  the  noninteracting  atoms. 


1.  Introduction 

Despite  the  manifold  successes  of  energy  band  theory  in  describing  a  wide 
spectrum  of  solid  state  properties,  the  ab  initio  determination  of  cohesive  energies 
and  equilibrium  lattice  constants  remains  as  a  relatively  unexplored  and  un¬ 
derdeveloped  field  of  study.  Thus,  there  has  developed  considerable  interest  in 
the  application  of  the  local  density  formalism  ILDF)  [1.2]  for  the  study  of  the 
ground  state  electronic  properties  of  molecules  and  solids.  Investigations  of  the 
cohesive  properties  of  small  moiecuies  [3-5],  metais  [6],  and  rare-gas  solids  [7] 
have  elucidated  the  possibilities  of  obtaining  a  reasonably  accurate  description 
of  binding  energies  and  equilibrium  geometries  by  incorporating  exchange  and 
correlation  effects  directly  into  a  one-bodv  potential.  Similar  studies  on  com¬ 
pounds,  ionic  insulators,  and  covalent  semiconductors  arc  beset  with  the  dif¬ 
ficulties  of  having  to  consider  the  full  ( nonmuffin-tin)  crystal  potential  ana  to 
explicitly  account  for  charge  redistribution  and  hy  bridization  processes  by  means 
of  a  fully  self-consistent  treatment.  Such  effects  cannot  be  conveniently  treated 
within  the  standard  Augmented  Plane- Wave  i  \PW)  or  Konn-Korringa-Rostoker 
(KK.R)  techniques  previously  used  for  such  studies. 

We  have  recently  developed  [8.  9]  ihe  fully  self-consistent  numerical  basis 
set  linear  combination  of  atomic  orbitals  I  LCaO)  discrete  variational  method 
(Dvm)  for  treating  ground  and  excited  state  properties  of  solids  in  the  LDF  ap¬ 
proximation.  This  scheme  permits  the  treatment  of  general  (i.e..  analytic  or 
numerical)  basis  functions  and  crystal  potentials,  and  the  determination  of  fully 
self-consistent  solutions  of  the  LDF  one-particle  equations  without  restricting 
the  iterative  path  to  a  superposition  of  spherical  charge  densities  [4]  or  to 
muffin-tin  models  (6,  7],  We  have  demonstrated  a  rapid  convergence  of  the 
self-consistent  (SC)  cycle  when  the  treatment  of  the  full  crystal  charge  density 
is  suitably  apportioned  between  real-space  and  Fourier-transformed  recipro- 
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Generalized  suscer*'  elites.  x°Cfl.  have  been  obtained  for  the  I  T  polymorphs  of  both  TaS,  and  TaS<i  from 
their  Korringi-Kohn-Rostoker  band  structures.  Peaks  are  found  which  correspond  to  the  charge-density-wave 
(CDW)  nesting  vectors  found  by  Wilson,  Di  Salvo,  and  Mahajan  and  appear  to  confirm  the  role  of 
electronically  driven  instabilities  as  the  origin  of  the  observed  CDW  in  ihese  metals. 


L  INTRODUCTION 


The  recent  discovery1,5  of  charge-density  waves 
(CDW)  in  the  layered  dichalcogenides  has  stimu¬ 
lated  great  interest  in  the  subject  of  electronically 
Induced  instabilities  in  metals  first  discussed  in 
the  pioneering  work  of  Overhauser.3  Along  with 
spin-density  waves  and  phase  transformations, 
CDW  are  thought  to  be  manifestations  of  these 
Instabilities  which  result  from  a  divergence  in  the 
generalized  static  electronic  susceptibility/  x(q). 
In  the  approximation  of  linear-response  theory/ 
X°(q)  may  be  written 

X(q)=X°(q)/[l  -/(q)x°(q)],  (1) 


where  /( q)  is  the  electron-electron  interaction  and 
X(q)  (the  bare  susceptibility)  is  defined  as 


X*(q)=  Z.  .  (2) 


Fermi  surface  “nesting,”  i.e.,  the  existence  of 
large  parallel  pieces  of  Fermi  surface,  leads  to  a 
divergence  in  x°(q)  and  possibly  to  a  divergence  in 
X(q)  itself,  depending  on  the  value  of  /( q). 

Electron-diffraction  studies1,5  have  revealed  that 
the  previously  observed’  anomalous  electromag¬ 
netic  properties  of  the  group-V  B  layered  transi¬ 
tion-metal  dichalcogenides  may  be  related  to  the 
formation  at  some  temperature  of  CDW’s  and  their 
accompanying  periodic  lattice  distortions.  This 
anomalous  behavior  has  been  associated  with  a 
nearly  two-dimensional  Fermi  surface  supporting 
a  CDW  formation  which,  in  turn,  introduces  a 
periodic  lattice  distortion  generally  incommensu¬ 
rate  with  the  lattice.1 The  amplitude  of  these  in¬ 
commensurate  waves  grows  from  a  high-tempera- 
ture  (below  550  K  in  IT-TaSe,)  down  as  far  as  the 
first-order  transition  temperature  which  repre¬ 
sents  the  point  of  conversion  (or  "lock-in”)  to  a 
commensurate  geometry.5  In  the  absence  of  any 


firm  theoretical  calculations  of  x(q)  on  these  ma¬ 
terials,  Wilson  et  aL.'  and  Williams  et  aL}  made 
use  of  the  augmented-plane-wave  results  of  Matt- 
heiss7  on  11-TaSj  available  to  them,  generated  a 
Fermi  surface,  and  interpreted  the  electron  dif¬ 
fraction  and  resistivity  data  observed  for  IT-TaSe,. 
This  Fermi  surface  is  nearly  two-dimensional  in 
the  basal  plane  {k,  =0)  with  its  walls  parallel  to 
the  FA  direction. 

This  paper  reports  results  of  detailed  ab  initio 
studies  of  xJ(q)  for  the  IT  polymorphs  of  both  TaS, 
and  TaSe?  derived  from  their  Korringa-Kohn- 
Rostoker  (KKR)  band  structures’  and  full  Brillouin- 
zone  Fermi  surfaces.  As  was  found  earlier/  the 
Fermi  surface  for  both  the  disulfide  and  the  di- 
seienide  are  very  similar,  having  cress  sections 
which  are  approximately  constant  in  planes  that 
are  perpendicular  to  the  z  axis  and  have  large  sec¬ 
tions  which  can  be  nested  by  approximately  the 
same  wave-vector  parallel  to  the  I\Vf  direction. 
However,  the  response  function  of  the  system 
contains  not  only  information  about  Fermi-surface 
(FS)  features  (“nesting”)  but  also  about  band- 
structure  effects  for  states  just  above  and  below  the 
Fermi  energy  (“volume  effects”).  Although  pre¬ 
vious  discussions  have  emphasized  FS  nesting 
features,  we  find  that  both  contribute  substantially 
to  the  important  structure  found  in  x3(q).  Accurate 
calculations  of  x°(q)  for  both  systems  show  major 
peaks  to  occur  at  those  q  values  corresponding  to 
the  CDW  vectors  found  by  Wilson  et  aL.'-  The  theo¬ 
retically  obtained  maxima  in  x°(q)  appear  to  pro¬ 
vide  a  detailed  confirmation  of  the  role  of  elec¬ 
tronically  driven  instabilities  as  the  origin  of  the 
observed  CDW  in  these  metals. 


II.  BAND  STRUCTURE  AND  FERMI  SURFACE 

The  energy-band  structure  of  both  lT-TaS,  and 
IT-TaSe,,  determined  by  the  KKR  method  in  the 
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A  new  approach  to  the  fully  self-consistent  solution  of  the  one-particle  equations  in  a 'periodic  solid  within  the 
Hohenberg-Kohn-Sham  local-density-t'unctional  formalism  is  presented.  The  method  is  based  on  systematic 
extensions  of  non-self-consistent  real-space  techniques  of  Ellis,  Painter,  and  collaborators  and  the  self- 
consistent  reciprocal-space  methodologies  of  Chaney,  Lin.  Lalon,  and  co-workers.  Specifically,  our  approach 
combines  a  discrete  variational  treatment  of  all  potential  terms  iCoulomb,  exchange,  and  correlation)  arising 
from  the  superposition  of  spherical  atomiciike  overlapping  charge  densities,  with  a  rapidly  convergent  three- 
dimensional  Fourier  senes  representation  of  all  the  multicenter  potential  terms  that  are  not  expressible  by  a 
Superposition  model.  The  bas's  set  consists  of  the  exact  numencal  valence  orbitals  obtained  from  a  direct 
solution  of  the  local-density  atomic  one-particle  equations  and  (for  increased  vanational  freedom)  virtual 
numencal  atomic  orbitals,  charge-transfer  (ton-pair)  orbitals,  and  “free''  Slater  one-site  functions.  The  initial 
crystal  potential  consists  of  a  non-muilm-tin  superposition  potential,  including  ncngradient  free-electron 
correlation  terms  calculated  beyond  the  random-phase  approximation.  The  usua:  multicenter  integrations 
encountered  in  the  linear-combination-of-atomic-orbitals  tight-binding  formalism  are  avoided  by  calculating  all 
the  Hamiltonian  and  other  matrix  elements  between  Bloch  states  by  three-oimer.sional  numencal  Dtophanune 
integration.  In  the  first  stage  of  self-consistency,  the  atomic  superposition  potential  and  the  corresponding 
numencal  basis  orbitals  are  modified  simultaneously  and  nun .; ecu r.y  bv  varying  (iteratively)  the  atomic 
occupation  numbers  (on  the  bas  s  of  computed  Bnllouin-rone  averaged  band  populations)  so  as  to  minimize 
the  deviation,  April,  between  'he  band  charge  density  and  the  superposition  charge  density.  This  step 
produces  the  "best"  atomic  corfiguranon  within  the  superposition  model  for  the  crystal  charge  density  and 
tends  to  remove  all  the  sharp  ''localized"  features  in  the  function  Apt  r  i  by  allowing  for  intra-atomic  charge 
redistn  tion  to  take  place.  In  the  second  stage,  the  tnree-dimensional  multicenter  Poisson  equation  associated 
with  April  through  a  Founer  senes  representation  of  April  is  solved  and  solutions  of  the  band  problem  are 
found  using  a  self-consistent  cntenon  on  the  Founer  coefficients  of  April.  The  calculated  observables  include 
the  total  crystal  ground-state  energy,  equtlibnum  lattice  constants,  electronic  pressure,  x-ray  scattenng  factors, 
and  direcnonal  Compton  profile.  The  efficiency  and  reliability  of  the  method  is  illustrated  by  means  of  results 
obtained  for  some  ground-state  properties  of  diamond;  comparisons  are  made  with  the  predictions  of  other 
methods. 


L  INTRODUCTION 

The  current  popularity  of  energy-band  theory 
stems  from  its  successful  application  to  the  study 
of  increasingly  diverse  problems  in  solid-state 
physics.  Recent  new  sophisticated  experiments 
on  both  traditional  materials  and  those  having 
complex  crystallographic  structures  have  de¬ 
manded,  however,  not  only  theoretical  descrip¬ 
tions  of  phenomena  related  to  the  one- electron 
eigenvalue  but  also  detailed  and  precise  wave 
f unctions  with  which  to  determine  the  expectation 
values  of  different  observable  operators.  Such  a 
demanding  test  of  the  predictions  of  one-electron 
theory  has  the  additional  virtue  in  permitting,  by 
their  comparison  with  experiment,  accurate  de¬ 
terminations  of  the  relative  magnitude  and  impor¬ 
tance  of  many- body  effects  in  real  solids.  Thus, 
there  has  developed  considerable  interest  in  ap¬ 
plying  the  Hohenberg-Kohn-Sham1, 11  local-density- 
functional  (LDF)  formalism  to  the  investigation  of 
various  ground-state  properties  of  solids,  despite 
the  usual  difficulties  of  solving  the  associated  one- 
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particle  equation  characterized  by  a  multicenter 
nonspherical  potential.  A  variety  of  well-known 
approximations  have  been  introduced  in  order  to 
reduce  the  complexity  of  the  problem.  For  exam¬ 
ple.  the  various  forms  of  orthogonalized-plane- 
wave  (OPVV)  techniques3,4  suffer  from  serious  con¬ 
vergence  difficulties  when  applied  to  solids  con¬ 
taining  first- row  atoms  with  no  p  states  in  the 
atomic  cores.  Recent  studies  with  extended  sets 
of  orthogonalized  plane  waves3  have  indicated  er¬ 
rors  in  the  eigenvalues  of  the  order  of  up  to  a  few 
volts  due  to  poor  convergence  in  calculations  em¬ 
ploying  several  hundreds  of  OPW’s.  The  muffin- 
tin  approximations  applied  to  augmented-plane- 
wave  (APW)  calculations6, 7  have  been  recently 
criticized  as  introducing  errors  of  up  to  a  rydberg 
in  the  potential  in  diamond8  and  other  covalent 
structures,3  and  overestimating  the  binding  ener¬ 
gy  in  some  covalently  ^bonded  molecules  by  as 
muchas(100-300)lc. 10  Recent  LCAO  (linear  orbit- 
als)-type  calculations11*13  have  overcome  thediffi- 
cultv  of  treating  non- muffin- tin  potentials  and  have 
demonstrated  that  efficient  convergence  with  respect 
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We  use  our  previously  reported  method  for  solving  self-consistently  the  local-density  one-particle  equations  in 
■  numencaJ-basis-set  linear  combination  of  atomic  orbitals  expansion  to  study  the  ground-state  charge  density, 
i-tty  structure  factors,  directional  Compton  profile,  total  energy,  cohesive  energy,  equilibrium  lattice  constant, 
and  behavior  of  one-electron  properties  under  pressure  of  diamond  Good  agreement  is  obtained  with  available 
experiment  data.  The  results  are  compared  with  those  obtained  by  the  restricted  Hartree-Fock  model  the 
role  of  electron  exchange  and  correlation  on  the  binding  mechanism,  the  charge  density,  and  the  momentum 
density  is  discussed 


L  INTRODUCTION 

The  local  density  functional  (LDF)  formalism 
of  Hohenberg,  Kohn,  and  Sham.1’-  and  its  recent 
extension  as  a  local  spin-density  functional  formal¬ 
ism,1  form  the  basis  of  a  new  approach  to  the  studv  of 
electronic  structure  in  that  the  effects  of  exchange 
and  correlation  are  incorporated  directly  into  a 
charge-density-dependent  potential  term  that  is 
determined  self-consistentlv  f  mm  the  solution  of 
an  effective  one-particle  equ.it ton.  Applications  of 
the  LDF  formalism  to  atoms"  and  molecules’  have 
yielded  encouraging  results.  Similar  applications 
for  solids  are  complicated  by  ti)  the  need  to  con¬ 
sider  both  the  short-range  and  the  long-range 
multicenter  crystal  potential  having  nonspherical 
components,  (ii)  the  difficulties  in  obtaining  full 
self-consistency  in  aperiodic  system,  and  (iii) 
the  need  to  provide  a  basis  set  with  sufficient 
variational  flexibility.  Hence,  theoretical  studies 
of  ground-state  electronic  properties  of  solids 
in  the  LDF  formalism  have  been  mainly  limited  to 
muffin-tin  models  for  the  potential,7’3  non-self- 
consistent  schemes,9  treatments  of  simplified 
jellium  models10  or  spherical  cellular  schemes.11’13 

We  have  recently  proposed13,1'1  a  general  self- 
consistent  method  for  solving  the  LDF  formalism 
one-particle  equation  for  realistic  solids  using  a 
numerical- basis- set  LCAO  (linear  combination  of 
atomic  orbitals)  expansion  and  retaining  all  non¬ 
spherical  parts  of  the  crystal  potential.  We  have 
demonstrated  a  rapid  convergence  of  the  self- 
consistent  (SC)  cycle  when  the  treatment  of  the 
full  crystal  charge  density  is  suitably  apportioned 
between  real-space  and  Fourier- transformed 
reciprocal- space  parts  and  have  indicated  the 
large  degree  of  variational  flexibility  offered  by 
a  nonlinearly  optimized  (exact)  numerical  atomic- 
like  basis  set  We  have  shown  that  all  multi¬ 
center  interactions  as  well  as  the  nonconstant  parts 
of  the  crystal  potential  are  efficiently  treated  by 


a  three-dimensional  Diophantine  integration 
scheme. 

The  purpose  of  this  paper  is  to  illustrate  the  ap¬ 
plicability  of  our  method  to  real  systems  by 
studying  the  ground- state  electronic  properties  of 
diamond.  Diamond  has  been  long  considered  as 
a  prototype  for  covalently  bonded  insulators-5  and 
a  great  deal  of  experimental  work  has  been  done 
on  its  ground-state  properties,  including  co- 
hensive  energy,15  lattice-constant  studies.17  x-ray 
scattering  factors.11’1’  charge  density.13  and  di¬ 
rectional  Compton  profile.30*31  In  addition, 
theoretical  studies  on  its  ground- state  properties 
within  the  restricted  Hartree-Fock  iRHF)  model 
are  available*’3*35  so  comparison  with  the  predic¬ 
tions  of  the  LDF  formalism  is  possible.  Although 
the  eigenvalue  spectrum  (band  structure)  of  the 
local  exchange  Hamiltonian  for  diamond  has  been 
studied  previously  by  a  variety  of  first-principles 
techniques  [augmented  planes  waves  (APW),33’ 11 
orthogonalized  plane  waves  (OPW),33*30  pseudo- 
potential  OPW,11  LCAO,33"3'  and  cellular  meth¬ 
ods11],  the  ground-state  observables  related 
to  the  ground- state  crystal  charge  density  have 
received  much  less  attention.  While  our  method 
was  shown13  to  accurately  reproduce  the  band 
structure  of  diamond  as  obtained  by  other  tech¬ 
niques33’33  (when  correlation  and  self-consistency 
is  omitted  so  as  to  be  compatible  with  pre¬ 
viously  published  band- structure  models  and  the 
full  nonspherical  components  of  the  potential  are 
retained  in  both  calculations),  we  do  not  consider 
this  as  a  stringent  test  since  the  LDF  formalism 
in  its  "standard”  form  does  not  make  any  claim 
on  the  physical  significance  of  the  band  eigen¬ 
values  nor  are  these  eigenvalues  sensitive  enough 
to  the  details  of  the  basis  set  and  potential.36  In 
what  follows  we  present  our'  results  for  the  x-ray 
scattering  factors,  charge  density,  directional 
Compton  profile,  total  energy,  and  equilibrium 
lattice  constant  and  discuss  the  role  of  exchange 
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A  fully-seif-consistent  numencal-basis-set  linear-combmation-of-aiomic-orbitals  calculation  of  the  electronic 
structure  of  TiSj  is  reported  using  the  method  described  previously.  The  calculated  band  structure  differs 
considerably  from  those  previously  obtained  by  non-seif-consistent  muiTin-tin  models  Comparison  with 
experiment  shows  that  the  calculated  optical  properties  for  energies  helow  16  eV  and  the  various 
characteristics  of  the  valence  and  conduction  bands  agree  very  well  with  opticai-absorption  and  electron- 
energy-loss  data  as  well  as  with  photoemission,  x-ray  absorption,  and  appearance-potential  spectra.  A  small 
indirect  gap  (0  2-0  2  eV)  occurs  at  the  points  M  and  L  in  the  Bniloum  rone  with  a  larger  direct  gap  iO  8 
eV)  at  T.  We  suggest  that  the  characteristic  scmi-metallic  large  g  value  observed  experimentally  originates 
from  a  near  coincidence  of  the  band  gap  with  the  enhanced  spin-orbit  splitting  which  is  consistent  with  the 
soft-x-ray  data  and  our  band  model  The  bonding  mechanism  in  TiS;  is  discussed  in  detail;  it  is  snown  by  a 
direct  calculation  of  the  self-consistent  charge  density  and  the  transverse  effective  charge  that  the  system  is 
predominantly  covalent  with  small  static  ionic  character  and  large  dynamic  lomciiy  In  contrast  with  mutTin- 
tin  Xa  models,  the  bonding  is  found  to  be  largely  due  to  Ti  ds ip  to  S  3 p  bones  and  a  much  weaker  Ti  id  to 
S  ip  bond.  The  effects  of  mutTin-tin  approximation  and  self-consistency  are  discussed  in  detail.  Extrapolation 
of  these  results  to  the  case  of  TiSe,  is  made  and  the  possible  origin  of  its  charge-density  wave  is  discussed. 


L  INTRODUCTION 

The  renewed  experimental  and  theoretical  inter¬ 
est  in  the  layered  transition-metal  dichalcogenides 
has  brought  about  in  the  last  few  years  much  great¬ 
er  understanding  of  their  unusual  anisotropic  opti¬ 
cal  and  transport  properties/  and.  more  recently, 
those  properties  associated  with  the  observation  of 
charge -density  waves'2  in  some  of  these  sy?for"e. 
Perhaps  an  outstanding  exception  is  lT-TiS:  for 
which  the  available  information  is  still  in  a  state 
of  flux  and  conflict.  For  example,  while  optical 
experiments3  have  suggested  that  TiS_,  is  a  semi¬ 
conductor  with  a  gap  of  about  1-2  eV,  recent  tem¬ 
perature-dependent  resistivity  measurements  on 
high-purity  stoichiometric  samples 1,5  demonstrated 
a  semimetallic  behavior  and  a  remarkable  temper¬ 
ature-dependent  electrical  resistivity  proportional 
to  TJ  from  at  least  10-400  K.  Similarly,  the  theo¬ 
retical  description  has  had  its  difficulties.  .An 
early  non-self-consistent  Korringa-Kohn-Rostoker 
(KKR)  calculation  in  the  muffin-tin  (MT)  approxi¬ 
mation9  showed  a  fundamental  gap  of  2.0  eV  at  r 
and  a  smaller  indirect  gap  of  1.4  eV  between  r  and 
L.  Further,  while  some  of  the  calculated  inter¬ 
band  transition  energies  m  the  MT-KKR  model 
could  be  reconciled  with  the  observed  data  only 
after  a  rigid  shift  of  about  1.4  eV  was  introduced 
in  the  theoretical  transition  energies,  an  empirical 
non-self -consistent  tight -binding  model7  s  with  pa¬ 


rameters  chosen  to  fit  these  observed  optical 
transitions  yields  (a)  a  valence-band  width  about  a 
factor  of  2  too  small  (relative  to  the  x-ray  pnoto- 
emission  data3"11)  and  (b)  a  similar  underestima¬ 
tion  of  the  d-conduction-band  width  (compared  with 
the  observed  appearance-potential  spectra13).  Sim¬ 
ilar  contradictions  occur  in  the  various  one-elec¬ 
tron  models  derived  to  explain  absorption  mech¬ 
anisms  in  TiSa  while  the  semiempirical  (optical 
data)  model  of  Murray  etal.~  suggests  the  occur¬ 
rence  of  a  nonoverlapping  metal  d  band  in  the  gap 
formed  between  the  bonding  and  antibonding  metal- 
nonmetal  s-p  bands  lo-a*  gap),  the  crystal-field 
model  of  Husiman  efaf.,13  suggests  that  the  metal 
d  band  appears  in  the  valence-band  continuum.  By 
contrast,  x-ray  photoemission  of  Williams  and 
Shepherd10  places  the  lower  part  of  the  d  band  in 
partial  overlap  with  the  valence  band,  while  a  sec¬ 
ond  component  of  this  band  appears  in  the  antibond¬ 
ing  a*  region. 

In  this  paper,  we  present  the  results  of  an  ab 
initio  theoretical  study  of  the  electronic  properties 
of  TiS,  using  our  self-consistent  numerical-basis- 
set  linear-combination-of-atomic-orbitals  (LCAO) 
method.1'1'1'’  This  method  uses  an  accurate  non- 
linearly  optimized  (exact)  atomic  basis  set  for  all 
core,  valence,  and  some  virtual  states  of  the  atoms 
appearing  in  the  unit  cell.  All  interaction  and  over¬ 
lap  integrals  are  evaluated  accurately  using  a 
three-dimensional  numerical  Diophantine  mtegra- 
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We  review  the  present  state  of  understanding  of  the  electronic  structure  and  physical  properties  of  actinide  metals  and 
inter  metallic  compounds  as  derived  from  relativistic  APW  energy  hand  studies  of  some  of  the  light  (Th.  U.  Np.  and  Pul  and 
heavy  (Am.  Bk.  and  Cm)  metals,  the  intermetallics  L  Rh,  and  L  Ir,  and  i  SaCI  structure  compounds  such  as)  UC.  Emphasis  is 
placed  on  the  importance  of  Coulomb  correlation  and  the  roie  of  actinide-actinide  separation  in  determining  the  itinerancv 
(as  in  the  light  metals)  or  localization  i as  in  the  heavy  metals  i  of  the  5f  electrons  and  in  turn  their  resulting  magnetic  and  other 
properties.  The  UX,  systems  considered  are  significant  because  the  uranium  sites  are  si'fficiently  separated  to  fall  in  the  local 
moment  f  orbital  range  of  the  Hill  superconductivitv/magnetism  plot  but  show  no  magnetic  ordering.  Comparisons  to  recent 
optical,  de  Haas-van  Alphen  and  other  data  are  given  when  available  The  contrasting  cases  of  the  transition  and  rare-ea rtn 
systems  make  it  clear  that  the  5f  electrons  are  a  unique  species  which  offer  exciting  challenges  to  both  experimentalists  and 
theorists. 


1.  Introduction 

The  last  two  decades  have  witnessed  remarka¬ 
ble  advances  in  our  understanding  of  the  elec¬ 
tronic  structure  and  properties  of  a  wide  range  of 
materials.  Sparked  by  the  introduction  and  use  of 
highly  sophisticated  experimental  techniques,  a 
mass  of  experimental  knowledge  has  emerged 
which  challenged  previous  (and  oftentimes  simp¬ 
lified)  theoretical  models  and  brought  about  a 
more  unified  view  of  electronic  structure  and 
behavior  First  the  nohle  metals,  then  the  impor¬ 
tant  transition  metals  and  their  compounds  and. 
most  recently,  the  rare-earths  have  each  enjoyed 
a  period  of  research  which  resulted  in  a  high 
degree  of  sophisticated  understanding  previously 
thought  unlikely.  Such  a  period  has  only  recently 
arrived  for  actinide  research — called,  for  many 
reasons  and  only  somewhat  euphemistically,  the 
“last  frontier"  in  the  study  of  the  electronic 
structure  of  matenalv  Until  recently,  the  difficul¬ 
ties  of  obtaining  samples  of  sufficient  purity  and 
overcoming  problems  associated  with  their 
radioactive  nature  have  resulted  in  slow  progress 
and  limited  understanding  In  the  last  decade  the 
more  intensive  experimental  efforts  have  fos- 
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tered  intensive  theoretical  study  of  the  basic 
phenomena  and  have  given  rise  to  a  more  unified 
understanding  of  their  electronic  properties  jtj. 

That  the  actinides  enjoy  their  own  unique 
position  among  the  elements,  arises  from  an 
unusual  set  of  circumstances.  In  the  transition 
metals  the  s  and  d  valence  electrons  form  the 
conduction  bands  which  are  responsible  for  their 
electric,  magnetic,  and  optical  properties.  The 
rare-earths  are  unique  in  that  their  4f  electrons 
are  so  highly  localized  that  whereas  they  deter¬ 
mine  the  various  exotic  magnetic  structures  and 
properties  of  their  metals,  they  have  little  effect 
on  other  chemical  and  physical  behavior  which 
arises  from  the  transition  metal-like  structure  of 
their  5d  and  6s  valence  electrons.  The  5f  elec¬ 
trons  in  the  actinide  elements  are  not  as  well 
localized  as  the  4f's  in  the  rare-earths  and  have 
energies  which  are  close  to  those  of  the  6d  and  7s 
electrons.  This  produces  the  unusual  condition 
of  a  strong  "competition"  in  actinide  solids  be¬ 
tween  the  5f  electrons  and  the  6d  and  7s  elec¬ 
trons  in  determining  their  electronic  structure 
and  properties.  The  contrast  between  the  results 
on  the  actinides  and.  those  on  transition  and 
rare-earth  metals  makes  clear  that  the  5f  elec¬ 
trons  are  a  species  of  their  own  [1],  w  ith  interest¬ 
ing  and  complex  properties  that  present  great 
theoretical  challenges. 
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We  find  thjt  a  defect  state  treatment  of  localized  c.sciiations  m  L 1 1  witlnn  the  ImcjI  vlensitv  funciion.il  formalism 
accounts  remarkably  'veil  lor  the  observed  expermii. ni.d  (core  plus  .ipiu.il  c.ipi  c\«.it.;linnv  •  m  cmnra^;  to  the  failure 
of  the  one-eLvtron  band  model  We  draw  that  when  electron  relaxation.  NclMutcr.iction  juJ  ch.irec  polarization  ettects 
are  taken  into  account  by  treatine  the  excitation  as  a  localized  point  delect,  the  imnrow  d  lund  model  predicts  the  cor¬ 
rect  excitation  at’.J  interband  states. 


Despite  extensive  investigations,  the  nature  ot'the 
fundament ji  optical  gap  and  eore  excitations  in  ionic 
solids  -■  of  which  L.i E  is  considered  the  prototv  pc 
remains  as  a  challenge  to  expeumentalists  and  theorists 
alike.  Theoretically,  two  basic  models  have  been  ex¬ 
tensively  applied:  the  restricted  I lartree-Fock  (1117) 
model  [1]  and  the  local  density  t  LD)  model  (2  - 5 1 . 
With  the  recent  advances  in  linear  combination  of 
atomic  orbitals  I  LCAO)  techniques,  the  more  conven¬ 
tional  non-self-consistent  muffin-tin  schemes  were 
abandoned  in  favor  of  the  more  sophisticated  extended 
basis  set  sell-consistent  (non-muftle-tin )  LCAO  methods 
for  both  models.  Their  applications  have,  however, 
yielded  mixed  results:  (i)  local  exchange  calculations 
[3]  with  an  exchange  coel ilcient  a  close  to  1 .0  could 
reproduce  the  optical  gap  but  an  extended  Gaussian 
basis  (a  =  I )  studs  [3]  gave  calculated  one-electron 
energies  substantially  lower  than  experiment  m  the 
interband  region  [4]  and  yielded  the  suggestion  that 
the  observed  spectra  both  in  the  optical  gap  region 
(I  1-12  eV)  and  in  the  Li-K  excitation  region  (60- 
62  eV)  be  reinterpreted  as  Bloch-type  interband  transi¬ 
tions  instead  of  as  bound  excitons  -  in  marked  contra¬ 
diction  with  recent  experiments  [4-8].  (n)  restricted 
HF  calculations  [1]  revealed  a  pronounced  disagree¬ 
ment  of  the  one-electron  eigenvalue  diffeiences  with 
experimental  transition  energies  in  (lie  whole  spectral 
region  and  indicated  that  electron  correlation,  electron- 
hole  interactions  and  relaxation  corrections  (calculated 
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using  a  simplified  atomic  model)  are  necessary  to  brina 
the  lesultc  inti'  agieement  with  experiment.  It  was  fur¬ 
ther  stated  by  these  authors  |  1  |  that  the  local  exchange 
model  is  inadequate  lor  describing  excitations  in  these 
materials. 

We  report  results  which  show  that  a  proper  treat¬ 
ment  of  locaii/ed  excitations  m  the  LI)  model  accounts 
icmaikably  well  for  all  the  observed  experimental  data 
in  contrast  to  the  failures  of  the  one-electron  baud 
model.  Our  method  goes  beyond  the  conventional 
baud  model  by  considering  excitation  processes  as 
transitions  involving  point-defecl  stales  in  a  solid  and 
uses  t.tta!  (statistical)  energy  diffeiences  between  sena- 
laielv  calculated  ground  and  exeted  states  rather  than 
<  •ih'-dcar:  <ii  eneigy  differences  ut  a  ground  state  cal¬ 
culation  to  evaluate  the  relevant  excitation  energies. 
Speciticallv .  this  is  done  in  our  “small  periodic  cluster" 
(SIX')  model  [l,|  m  which  wc  perform  a  fully  self- 
consistent  band  structure  calculation  but  with  a  large 
crystallographic  unit  cell  (8-l(>)  atoms)  containing  a 
locally  excited  atom  at  its  center  instead  of  the  usual 
(perfect  crystal)  primitive  cell.  The  model  allows  for 
explicit  electron-relaxation  effects,  self-interaction  cor¬ 
rections  (for  the  non-existence  of  Koopmans'  theorem 
m  the  local  density  model  even  in  the  unrelaxed  limit), 
charge  polarization  coriectioii  and  correlation  (treated 
in  the  free-electron  approximation  I  eft  0  c  t  s . 

Our  starting  point  is  the  self-consistent  (SC)  band 
struciure  obtained  by  solving  the  one-particle  LD  func¬ 
tional  Hamiltonian  with  a  free-electron  exchange  po¬ 
tential  (using  a  =  2  ;3  )  and  the  electron  correlation  po¬ 
tential  ot  Smgwi  et  al.  1 1 0 1 .  Our  LCAO  basis  set  eon- 
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The  band  structure,  charge  density,  x-ray  scattering  factor  (and  their  behavior  under  pressure), 
equilibrium  lattice  constant,  and  cohesive  eneray  of  the  prototype  ionic  solid  LiF  were  determined  using  our 
recently  developed  self-consistent  numerical  basis  set  (non-muffin-tmi  linear-combmation-of-atomic-orbitals 
method, within  the  local-density  formalism  iLDF)  The  details  ot  the  bonding  and  the  etfects  of  exchange  and 
correlation  on  the  electronic  structure  are  discussed  with  reference  to  the  conventional  picture  of  ionic 
bonding.  Remarkably  good  agreement  is  found  with  the  observed  data  for  the  ground-state  properties  of  the 
system.  Contrary  to  the  results  of  previous  band  studies,  the  conventional  hand-structure  approach  to 
excitation  energies  (i.e..  identifying  them  with  the  band  eigenvalue  differences!  is  found  to  fail  completely  in 
accounting  for  the  observed  data  in  tne  entire  x-ray  unu  optical  spectral  region  wnen  fully  seif-consistent 
solutions  of  the  LDF  one-particle  equation  with  no  further  approximation  to  the  crystal  potential  are 
obtained.  It  is  found  that  in  the  presence  of  some  spatial  localization  of  the  initial  or  final  crystal  states,  the 
spurious  self-interaction  terms,  as  weil  as  the  polarization  and  orbital  relaxation  se;f-energy  effects  arc  of  a 
similar  order  of  magnitude  as  the  Koopmans'-like  interband  terms.  Irt  oruer  to  treat  these  effects  within  the 
LDF  self-consistcntly,  we  describe  the  excitation  processes  as  .-ansiuons  involving  poim-defect-ltke  stales  in 
the  solid  calculated  by  a  supercell  method  in  which  the  excitat.cn  energies  are  determined  as  total-energy 
differences  between  (separately  calculated)  excited-  and  ground-state  configurations.  The  excited  state  is 
represented  as  a  superiattice  of  locally  excited  sues  using  large  ■  v  and  i  tom  i  uut  ceils,  cacti 
containing  a  single  excited  sue.  We  find,  in  the  self-consistency  limit,  that  a  small  but  finite  degree  of 
spatial  localization  of  the  excited  states  exists  even  for  valence  excitations,  inducing  thereby  self-interaction 
as  well  as  self-energy  relaxation  and  polarization  effects  The  LDF  model  is  found  to  account  verv  well  for 
both  interband  and  exciton  transitions  over  the  enure  spectral  region  tlZ-b95  eV)  and  to  yield  definite 
predictions  regarding  the  exciton  bandwidtbs  and  senes  limits 


L  INTRODUCTION 

The  study  of  the  ground-state  electronic  prop¬ 
erties  of  molecules  and  solids  in  the  local-density 
formalism  (LDF)‘‘J  has  recently  attracted  con¬ 
siderable  attention.  Investigations  of  the  cohesive 
properties  of  small  molecules, metals, %i!  and 
rare-gas  solids1  have  elucidated  the  possibilities 
of  obtaining  a  reasonably  accurate  description  ot 
binding  energies  and  equilibrium  geometries  by 
incorporating  exchange  and  correlation  effects  di¬ 
rectly  into  a  one- body  potential.  Similar  studies 
on  compounds,  ionic  insulators,  and  covalent 
semiconductors  are  beset  with  the  difficulties  of 
having  to  consider  the  full  (non-muffin-tm)  crystal 
potential  and  to  explicitly  account  for  charge  redis¬ 
tribution  and  hybridization  processes  by  means  of 
a  fully  self-consistent  treatment.  We  have  recently 
developed10,11  the  fully  self-consistent  numertcal- 
basis-set  linear  combination  ot  atomic  orbitals 
(LCAO)  discrete  variational  method.  This  scheme 
permits  the  treatment  of  general  (analytic  or  nu¬ 
merical)  basis  functions  and  crystal  potentials  and 


the  determination  of  fully  self-consistent  solutions 
of  the  LDF  one-particle  equation  without  restrict¬ 
ing  the  iterative  path  to  mulfin-ttn  charge  densi¬ 
ties'"1  or  to  superposition  models.  Applications 
of  this  approach  include  the  stuav  of  the  ground- 
state  properties  of  diamond. ,J  boron  nitride. u 
titanium  disulfide, ;l  cadmium  suiprud.  and  LiAl. 

The  present  study  extends  our  treatment  to  tonic 
solids,  for  which  LiF  has  beer,  chosen  as  a  c-'-.o 
type.  We  first  consider  the  description  of  ground- 
state  properties  of  the  system,  such  as  the  band 
structure,  charge  density,  x-rav  scattering  fac¬ 
tors,  cohesive  energy,  equilibrium  lattice  con¬ 
stant,  and  behavior  under  pressure,  and  compare 
the  predictions  of  the  LDF  model  with  both  experi¬ 
mental  data  and  with  available  restricted  Hartree- 
Fock  (HF)  results.1'  We  then  consider  the  excited- 
state  properties  of  the  system  by  first  using  the 
standard  band  approach  to  excitation  energies 
(i.e..  viewing  them  as  differences  in  the  band 
eigenvalues  between  unoccupied  and  occupied 
bands).  Contrary  to  what  has  Deen  previously 
suggested, 16'18  we  find  that,  whereas  this  approach 
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Abstract 

Theoretical  determinations  of  spin  densities.  A.  J-  Freeman  (Physics  Depart¬ 
ment.  Northwestern  University.  Evanston.  III.  <.0201.  USA). 
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This  paper  describes  and  discusses  some  aspects  of  the  present  status  of 
theoretical  determinations  oi  spin  unu  charge;  densities  in  atoms,  rr.oleeuiar 
clusters, and  solids.  Seleced  results  of  recently  completed  calculations  are  used 
to  illustrate  "state  of  the  art"  methodology  and  accomplishments. 

1.  Introduction 

In  this  paper,  I  will  review  several  aspects  of  the  present  status 
of  theoretical  determinations  of  spin  (and  charge)  densities,  ob¬ 
viously  of  high  interest  to  the  Conference,  by  means  of  examples 
which  will  give  a  flavor  of  recent  developments  in  atoms,  molecu¬ 
lar  clusters  and  solids.  Because  of  space  limitations  and  for  sound 
pedagogical  reasons,  the  examples  given  will  be  limited  to  those 
investigations  which  have  been  carried  out  in  my  group.  A  more 
detailed  description  of  theoretical  methods  and  results  is  given 
elsewhere  [1],  Emphasis  will  be  placed  on  neutron  magnetic  form 
factor  studies  as  these  give  directly,  through  the  Fourier  trans¬ 
form,  the  magnetization  densities.  On  the  other  hand,  hypertine 
interaction  measurements,  which  provide  information  about 
br>th  charge  spin  distributions,  should  not  be  overlooked  by 
the  active  researcher. 

For  this  discussion,  it  is  convenient  to  write  the  charge  density 
ec(r)  and  spin  density  cs(r) 

ec (') 

8,(1-)  -e,t)~°,i(r)  (1) 

where  ec0(r)  and  p30 (r)  are  respectively  the  charge  and  spin 
density  for  an  electron  of  given  spin  orientation.  a(  -  1  or  I).  For 
generality,  I  have  written  o _.„(/•)  =  oi0(r)  although  if  a  single 
calculation  were  done,  then  the  equality  would  hold.  As  is  well 
known  unlike  the  case  of  neutrons  or  X-rays,  the  limitation  of 
the  hyperfine  measurements  as  a  probe  is  that  they  measure  the 
value  of  and  £>,  only  at  one  point,  namely,  at  the  position  of 
the  nucleus:  the  isomer  shift  measures  the  total  charge  density  at 
the  nucleus,  <?c(o);  for  a  magnetic  hyperfine  measurement  in  which 
core  polarization  effects  are  dominant,  the  relevant  quantity  is 
the  net  spin  density  at  the  nucleus. 

The  state  of  theory  will  be  the  subject  of  the  remainder  of  my 
talk.  Because  of  its  great  success  in  explaining  a  variety  of 
phenomena  in  solids,  there  arose  great  unrealistic  expectations 
that  theory  could  and  would  give  detailed  answers  to  the  ques¬ 
tions  raised  by  the  experimentalists.  Instead,  progress  has  been 
slow  because  of  the  difficulties  inherent  m  treating  an  inter¬ 
acting  many-body  system  which  cannot  be  solved  selectively 
for  the  particular  information  needed  by  the  experimentalist 
Fortunately,  a  number  of  developments  in  the  last  five  years  show 


i 

great  promise  for  more  rapid  progress  to  be  made  in  the  future,  and 
hence,  one  can  therefore  be  somewhat  more  optimistic  about 
the  degtee  of  theoretical  support  that  can  be  provided. 


As  has  been  the  case  for  some  time,  theory  is  in  good  shape  for 
treating  free  atoms  and  ions,  thanks  to  the  highly  successful 
seif-consistent-fieid  models.  The  newest  results,  dealing  with 
isomer  shifts,  isotope  shifts,  core  polarization  and  neutron 
magnetic  form  factors  have  been  obtained  from  fully  reiatiwstic 
Dirac-Fock  iDF)  calculations.  Here  we  can  rnemion  omy  some 
examples  of  this  work.  Free  atom  results  are  particularly  useful 
for  complex  materials  for  which  ah  initio  calculations,  of  the 
molecular  cluster  and  solid  state  types  to  be  described.  later  are 
clearly  impossible. 

2. 1 .  Relativistic  restricted  Dirac-Fock  results 

Recent  efforts  in  obtaining  accurate  free  ion  charge  densities 
have  centered  on  calculations  using  numerical  Hartree-Fock  or 
Dirac-Fock  theories  (with  exact  excnangei  or  their  approximate 
counterparts.  Hartree-Fock-Slater  and  Dirac-Fock-Siater  theo¬ 
ries  (when  the  Slater  free-eiectron  exchange  is  used).  A  systematic 
study  of  relativistic  electron  densities  and  isomer  shifts  in  3d. 
■id.  and  5 J  ions  has  been  made  recently  [2]  using  the  numerical 
self-consistent  field  program  of  Desciaux  [3].  Both  multicon¬ 
figuration  D-F  (for  both  point  and  finite  nuclei)  and  H-F 
calculations  of  the  changes  in  niO)  were  performed  (with  A.r-' 
assumed  constant)  for  each  set  of  ions,  thus  obtaining  relative 
values  of  the  isomer  shift  SAB  in  a  consistent  approximation. 
Thev  then  compared  5AB  values  with  experimental  isomer  shifts 
in  order  to  separate  the  nonrelativistic  and  relativistic  atomic 
contributions  to  the  isomer  shift  from  those  contributions  due 
to  the  chemical  environment.  In  each  of  the  cases  considered, 
these  workers  [2]  calculated  the  ionic  ground  states.  For  the 
Hartree-Fock  problem,  this  was  the  single-configuration 
Hund's  rule  ground  state  in  Russell-Saunders  coupling.  How¬ 
ever,  the  relativistic  single  particle  wavefunctions  arc  eigenfunc¬ 
tions  of  the  totjl  angular  momentum,  but  not  of  the  separate 
spin  and  orbital  angular  momenta;  thus,  the  Dirac-Fock  calcula¬ 
tions  are  performed  in  j-j  coupling.  However,  the  total  wave- 
function  is  not  a  single  j  -  j  configuration  but  rather  a  linear 
combination  of  j-j  configurations  of  total  angular  momentum 
J  -  L  -  S.  This  multiconiiguration  Dirac-Fock  wavefunction 
formally  approximates  as  closely  as  possible  the  I L-S)  Hartree- 
Fock  wavefunction. 

Dirac-Fock  caiculaiions  have  been  performed  by  Desciaux 
and  Freeman  [4]  for  the  divalent  and  trivalent  ions  of  the 
rare-earths.  The  importance  of  relativistic  effects  in  ihe  rate- 
earths  has  been  demonstrated  in  a  dramatic  way  by  means  of 
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A  first-principles  energy-band  study  of  the  metallic  rutile  phase  of  VO;  using  a  general  crystal  potential 
usd  an  expansion  of  the  Bloch  functions  in  a  linear  combination  of  atomic  oroiiais  is  reported.  The  results 
are  compared  with  previous  work  and  experimental  optical,  x-ray  absorption  and  emission,  arid  x-ray 
photoelectron  spectroscopy  data.  W{  obtain  a  large  density  of  stales  at  the  Fermi  energy;  the  Fermi  surface 
is  found  to  be  determined  by  the  two  lowest  d  bands,  at  the  bottom  of  the  "t  ”  manifold  which  is  split  by 
the  orthorhombic  field,  the  lowest-band  Fermi  surface  posses-.es  xome  nesting  features  corresponding  to  a 
nesting  vector  3  =  f  R  A  calculation  of  the  generalized  susceptibility  in  ihe  constant-matnx-elemem 
approximation  shows  the  existence  of  a  maximum  ai  the  zone  boundary  R  We  -uggest  that  the  formation  of 
a  charge -density  wave  wuh  wave  vector  q  =  T  R  accompanied  by  a  rer:<\::c  tatit.e  distortion  is  thus  possible, 
the  subsequent  condensation  of  phonons  at  the  point  R  eouid  then  explain  the  crystallographic  phase 
transition  observed  at  T  =  339  K 


L  INTRODUCTION 

Vanadium  dioxide  undergoes  a  crystallographic 
phase  transformation  at  T.  =  339 1  1  K  from  the 
monoclinic  distorted  rutile  structure  to  the  tetra¬ 
gonal  rutile  structure  which  is  accompanied  by  an 
abrupt  jump  in  the  metallic  conductivity1”  from 
10'3  to  10*  (0  cm)'1]  and  a  jump  in  the  magnetic 
susceptibility'1  [from  6.5  <  10'5  to  6.7  x  10" '  emu 
mole].  This  semiconductor  to  metallic  phase 
transformation  has  been  studied  very  extensively 
both  experimentally  and  theoretically.1  Unlike 
VjO-,  and  other  transition-metal  compounds  known 
as  "Mott-Hubbard  insulators"  which  show  antifer¬ 
romagnetic  ordering  in  the  low-temperature  phase. 
V03  remains  paramagnetic  below  T.  and  no  long- 
range  magnetic  order  is  observed  either  in  mag¬ 
netic  susceptibility.'  MossbauerU  or  nuclear-mag¬ 
netic-resonance  experiments.'  The  semiconducting 
monoclinic  phase  of  VO.  is  characterized  by  a 
pairing  of  the  vanadium  atoms  along  the  c  axis  re¬ 
sulting  in  a  doubling  of  the  unit-cell  size:  along  the 
c  axis  one  finds  alternatively  short  and  long  metal- 
metal  distances3  -1;  at  the  same  time,  each  vana¬ 
dium  "pair’’  is  tilted  with  respect  to  the  c  axis, 
the  vanadium  atoms  being  no  longer  at  the  center 
o i  the  distorted  ligand  octahedron.  It  Is  also  worth 
noting  the  large  further  distortion  of  the  octahedral 
geometry  in  the  low-temperature  phase,  with  an 
important  shortening  of  some  of  the  V-O  bonds.' 
This  pairing  of  the  vanadium  atoms  can  explain  the 
lack  of  localized  moments  and  account  for  the  mag¬ 
netic  properties  and  possibly  some  of  the  electri¬ 
cal  properties  of  the  semiconducting  phase,  as  the 


electrons  may  now  be  trapped  in  homopolar 
bonds.1'’  This  picture  has  been  also  useful  in  ac¬ 
counting  for  some  features  of  the  magnetic  be¬ 
havior11  of  the  alloy  V._rNb.O,  with  0.5. 

V  Electronic  structure  models 

These  structural  peculiarities  led  Goodenough1’ 
to  propose  a  model  of  the  electronic  structure 
which  can  explain  some  properties  of  the  metallic 
and  semiconducting  phases.  This  model  is  based 
on  molecular  orbital  theory  in  conjunction  with 
crystal  field  and  chemical  bonding  ideas.  In  the 
metallic  state,  the  interaction  of  cations  along  the 
c  axis  gives  rise  to  a  narrow  il  band  which  is 
overlapped  by  a  wider  z •  band  formed  by  antibond- 
ing  vanadium  d  and  ixvgen  p.  states.  1 1  the  low- 
temperature  phase,  the  crystal  distortion  is 
thought  to  be  responsible  for  three  effects:  a)  the 
t*  band  rises  aoove  rhe  Fermi  level  due  to  a 
stronger  vanadium  d-nxygen  p.  interaction,  re¬ 
sulting  from  the  shortening  of  the  distance  of  the 
vanadium  to  one  oxygen  atom;  tii)  the  J  band 
splits  due  to  the  pairing  of  V  atoms  along  the  c 
axis;  < iii )  the  lowest  part  of  the  d  band  is  filled 
and  separated  by  a  small  gap  from  the  upper  part, 
thereby  accounting  for  the  semiconducting  prop¬ 
erties  of  the  compound. 

The  essential  features  of  the  electronic  struc¬ 
ture  of  the  metallic  phase  of  VO.  proposed  by 
Goodenough1'1  have  been  confirmed  by  the  band- 
structure  calculation  of  Caruthers  etal.,1'  using 
the  combined  auymented-plane-wave-linear-com- 
birution-of-atomic-orbitals  tAPW-LCAO)  inter- 
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Results  of  ab  initio  self-consistent  energy-band  calculations  on  the  full  IS-atom/iunt- 
cell  ternary  alloys  .’/Rfc.B,  i\f  =Y,  Er  and  Ho'  are  reported  and  used  to  discuss  their 
superconducting  and  magnetic  properties — including  the  observation  of  reentrant  mag¬ 
netism  in  the  normal  state  at  T  <0.9  K  in  ErRh,B.. 

The  recent  discovery1  of  reentrant  magnetism  ing  from  -  2.5  to  -  12  K  and  noted  the  unexplained 

(at  T  =  0.9  K)  accompanying  a  transition  to  a  nor-  switchover  from  magnetism  to  superconductivity 

mal  metallic  state  in  ErRh4B,,  a  ternary  rare-  in  going  from  Ho  to  Er.  The  later  discovery  of  a 

earth  superconductor  with  7\5  =  8.7  K,  has  gener-  transition  from  superconductivity  to  magnetism" 

ated  considerable  interest  in  the  origin  of  this  in  the  same  system  upon  lowering  the  tempera- 

new  phenomenon.  It  has  also  raised  important  ture  of  the  sample  was  surprising.  For  our  stud- 

questions  regarding  our  understanding  of  magne-  ies  we  therefore  chose  ,l/  =  Ho  and  Er,  and  also 

tism  and  superconductivity  and  the  nature  of  their  Y  (which  has  no  4/  electrons),  and  a  T/  which  is 

possible  coexistence.  This  Letter  reports  results  almost  exactly  equal  to  that  of  the  Lu  alloy  'which 

of  ab  imtio  seLf-consistent  (SC)  energy-band  cal-  has  a  full  4/  shell). 

culations  carried  out  on  three  of  these  4ERh.B4  The  band  calculations  for  the  full  18-atom  unit  - 

systems,  with  ,V/  =  Er,  Y  (an  11.3-K  superconduc-  cell  structures  were  performed  seif-consistently 
tor),  and  Ho  (a  ferromagnet  below  6.6  K),  in  their  using  the  linear  muffin-tin  orbital  (LMTO)  meth- 
paramagnetic  states.  Very  similar  total  and  or-  od  and  other  related  methods.3  The  calculations! 

bital  (l)  angular  momentum  decomposed  atomic  scheme  is  essentially  the  same  as  the  one  used 

densities  of  states  i DOS )  are  found  for  the  three  earlier  for  a  number  of  .4-15  compounds4  cut  ex¬ 
systems  investigated.  The  peak  in  the  total  DOS  tended  to  cover  also  ternary  systems.  The 

*  at  title  Fermi  energy'  arises  from  the  4,1  states  at  LMTO  method3  is  a  computationally  rapid  and 

the  Rh  site  and  appears  to  be  responsible  for  the  efficient  method  for  obtaining  insight  into  the  ma- 
superconductivity  of  the  Y  and  Er  systems.  The  jor  physical  aspects  of  energy-band  structure, 

lower  5 d  DOS  at  the  rare-earth  site  could  give  By  using  the  logarithmic  derivative  as  a  parame- 

rise  to  magnetic  ordering  via  the  4/-5H  ion  the  ter,  the  usual  energy  dependence  of  the  elements 

rare-earth  site)  Ruderman-Kittel-Kasuva-Yosida  in  the  eigenvalue  matrices  can  be  neglected  to  a 
(RKKY)  interaction  in  HoRh  B,  (at  a  temperature  good  approximation.  Moreover,  by  replacing  the 
above  its  otherwise  possible  superconducting-  Wigner-Seitz-cell  integrations  by  integration 

transition  temperature  T, 5)  and  at  a  temperature  over  overlapping  Wigner-Seitz  spheres,  the  na- 
mll  below  the  observed  T, 1  in  ErRh,B4.  trix  elements  are  very  easily  set  up.  All  this 

In  this  new  group  of  tprnarv  borides,  MRh,B4,  makes  the  LMTO  method  much  more  rapid  than 

Matthias  et  alf  originally  reported  the  existence  the  augmented-plane-wave  and  Kohn-Korrinca- 

d  either  ferromagnetism  (M  =  Gd,  Tb,  Dy,  and  Rostoker  methods,  without  much  sacrifice  of  ac- 

Ho)  or  superconductivity  (M  =Y,  Nd,  Sm,  Th,  Er,  curacy. 

Tra,  and  Lu),  with  transition  temperatures  rang-  The  potential  contained  the  Hedin-Lundqvist 
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CONDUCTION  El  FCTRON  POLARIZATION.  MAGNETIZATION  DENSITIES  AND 
NEUTRON  MAGNETIC  FORM  FACTORS  IN  METALS* 
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The  problem  of  selling  up  ‘‘itinerant’'  and  "localized"  descriptions  of  3d.  -If  and  5f  electron  states  is  discussed.  Recent 
band-theory  predictions  of  conduction-electron  polarization,  maunetiz.ition  density  and  the  neutron  magnetic  form  tactor 
for  various  metals  are  described.  Comparisons  of  theory  with  experiment  show  good  agreement,  and  confirm  our  ideas 
on  the  range  of  validity  of  the  energy -band  description. 


1.  Introduction 

years  ago,  jf  we  had  been  asked^what 
arc  jhe  criteria  of  'itinerant-electron  mag¬ 
netism”  (the  title  of  this  conference),  most  of  us 
would  have  responded  with  n  discussion  ip.  terms 
of  the  spatial  extent  of  the  wave  function,  or  tin 
the  case  of  magnetic  carriers)  the  extend-*<ifK.-xs 
of  the  spin  density  distribution.  A  comparison 
between  theories  based  on  atomic  wave  func¬ 
tions  (as  representative  of  the  localized  limit) 
and  free-electron  wave  functions  ;as  represen¬ 
tative  of  the  other  extreme)  would  have  been 
included.  However,  with  the  accumulation  of  a 
mass  of  data-on  the  experimental  side  notably 
from  neutron  magnetic  scattering  investigations, 
and  on  the  theoretical  side  notably  from  energy 
bund  studies  of  transition  and  rare-earth  metals- 
il  became  gradually  recognised  that  this  simple 
picture  of  "itinerant”  and  "localized”  states  was 
unsatisfactory.  In  recent  years,  a  more  sophisti¬ 
cated  and  more  correct  description  has  been 
developed  which  instead  focuses  on  the  degree 
of  extension  or  localization  in  energy  of  the 
state  in  question.  Simply  pul,  whether  a  lo¬ 
calized  or  itinerant  description  ts  to  be  used 
depends  on  which  of  two  competing  terms  is 
dominant-Coulomb  correlation,  and  the  effec¬ 
tive  band  width  arising  from  the  overlap  of  d  or 
f  orbitals  on  neighbouring  atoms  and  their  hy¬ 
bridization  with  orbitals  of  differing  angular 
momentum  character.  In  the  intermediate  case 
where  neither  elfeci  can  be  neglected,  tn- 

*  Supported  by  the  National  Science  Foundation,  the  Air 
Force  Office  of  Scientific  Re'C.ireh  and  the  Energy 
Research  and  Development  Administration. 


terestir.g.  complex  and  unusual  properties 
result-as  in  the  actinides  [I], 

With  the  foregoing  as  the  stated  theme  of  this 
paper.  1  will  consider  several  aspects  of  the 
band  model  which  are  important  for  an  accurate 
one-particle  description  of  m..gnetiz;.: u.n  den¬ 
sities  in  metals.  Outstanding  among  thc-e  ,s  ;r.a 
consideration  that  eigenvalues  are  not  as  sen¬ 
sitive  as  ume  functions  to  the  approximate 
potential  used  in  the  calculation,  so  that  eacrgx 
bands  by  themselves  are  less  likely  to  provide 
detailed  explanations  of  electronic  tntei actions 
in  solids  than  are  studies  where  the  wave  func¬ 
tion  is  an  essential  element.  Thus  the  in¬ 
vestigation  of  phenomena  requiring  wave  func¬ 
tions  for  their  description  may  well  be  the  key 
to  assessing  the  limits  of  the  single-particle 
model  and  to  understanding  the  role  of  rrur.y- 
body  effects  in  metals.  It  is  for  this  reason  chat 
so  much  significance  attaches  to  the  recent  de¬ 
velopments  in  band-theory  technique  which 
nowadays  allow  us  to  make  meaningful  ario 
detailed  comparisons  with  experiment  involving 
the  wave  functions. 

However,  whereas  a  variety  of  experimental 
methods  yield  information  that  has  been  related 
to  the  electronic  bund  structure,  there  are  few 
kinds  of  experiment  which  yield  direct  in¬ 
formation  about  the  wait?  functions  in  mct„!s. 
Perhaps  the  one  outstanding  exception  i.  the 
determination  of  the  neutron  magnetic  form  ''ac¬ 
tor,  which  gives  the  Fourier  transform  ot  the 
magnetization  density.  Over  the  past  few  venrs. 
accurate  neutron-scattering  experiments  n.ivc 
show n^p1' the  inadequacy  of  the  communis -U'od 
free  atom  and  free  ton  Hartree-Fock  form  fac- 
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We  have  studied  the  possible  role  of  electronic  structure  on 
observed  phonon  anomalies  and  phase  transformations  by  means 
of  accurate  calculations  of  the  generalised  resoor.se  function, 
X(q),  for  non-interacting  conduction  e _e  oor  :r.s  . '  Sotn  Fermi 
surface  "nesting"  (parallel  sections  o:  Rer.ri  surface)  ar.i  so- 
called  "volume"  etfects  \,cara_lel  ele  :trsr.-.::le  bands  crossing 
the  Ep)  are  included  using  results  :f  esorsy  oar.u  calculations 
and  our  recently  developed  analytic  tetrar.edrcr.  linear  energy 
method.  Correlations  of  the  positions  of  phonon  anomalies 
and/or  the  occurence  cf  tnase  transformations  v:itn  large  maxima 
in  x(?)  are  found  ar.d  provide  strong  evidence  that  they  arise 
from  an  electronically  driver,  divergence  m  \  T)  :  ('*!'■  Lr.  1 1- 
TaS2  ar.d  1T-la3e2>  _arge  peaxs  in  x  are  found  at  she  obser¬ 
ved  c  values  along  7M  corresponding  sc  sr.e  ZZ'i .  )2)  In  bosh  Fi 
and  Pt^metals,  a  large  peak  is  found  as  precisely  Shat  c  value 
along  [110]  for  which  large  Rohr  anomalies  have  ceer.  found.  Ar. 
additional  peak  is  found  along  [lily  which  predicts  the  occur¬ 
rence  of  a  Kohr.  anomaly  is  this  direction  v  out  wnich  is  strong¬ 
er  in  Pd  than  in  pt ) .  (3)  Large  peaxs  in  \{\)  are  found  for 
both  NbC  and  TaC  at  precisely  those  c  values  at  which  soft  mo¬ 
des  were  observed  cy  Smith  and  Fiaser.  Maxima  in  x(?)  are  pre¬ 
dicted  for  other  directions.  The  locus  of  these  c-3x  values 
can  be  represented  by  a  warped  cube  of  dimension  ~/.2  (2U/a) 
in  momentum  space  -  in  striking  agreement  with  the  soft  mode 
surface  proposed  phenomenologically  by  ’.Veber.  No  such  peaks 
are  fount!  for  ZrC*  ar.d  HfC,  for  which  r.o  pncr.cn  anomalies 
have  been  observed.  (4)  For  the  tetragonal  rutile  phase  of 
V02>  we  find  a  large  maximum  in  x(l)  at  the  cone  boundary 
Thi3  suggests  the  possible  formation  of  a  CDW  with  q=FR.  The 
subsequent  condensation  of  phonons  at  R  could  then  explain 
■►he  crystallographic  (metal  to  semiconducting |  erase  transi¬ 
tion  at  1=340.  Recent  x-ray  diffuse  scattering  measurements 
by  Terauchi  ar.d  Cohen  of  Northwestern  confirm  these  predic¬ 
tions;  their  results  show  a  lattice  instability  near  the 
metal-insulator  transition  to  occur  at  the  R  point. 


1  .  INTRODUCTION 

Recent  theoretical  and  experimental 
developments  have  called  attention  to  the 
central  role  played  by  the  generalised  sus¬ 
ceptibility  function  x(?)  in  the  understand¬ 
ing  of  cany  physical  phenomena  in  solids. 
Since  x(l)  measures  the  response  of  the  sys¬ 
tem  to  an  external  (generally  spatially  in¬ 


homogeneous)  perturbation,  emphasis  has  been 
placed  on  possible  anomalous  cehaviour  which 
may  result  from  an  instability  in  the  conduc¬ 
tion  electron  gas  as  a  result  of  a  divergent 
in  this  response  function.  In  linear-re  spor.s« 
theory 

t(?)  =  x(m)/0-i(m)x(T))  (1) 

where  I(q)  is  the  electron-electron  ir.ter- 
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We  report  a  detailed  augmented-plane-wave  energy-band  study  and  wave-function  analysis  of  stoichiometnc 
PdH  which  shows  that,  even  though  the  Fermi  surface  of  PdH  is  qualitatively  similar  to  that  of  silver,  the 
ample  “proton  model"  is  not  valid.  Instead,  the  screening  of  the  proton  in  PdH  is  found  to  be  larger  than  in 
aa  isolated  H  atom  due.  in  pan.  to  the  formation  of  a  H-Pd  bonding  band  below  the  bottom  of  the  a -band 
complex.  This  result,  which  is  in  qualitative  agreement  with  Switendick's  earlier  calculation,  is  confirmed  by 
ultraviolet  photoemission  experiments.  A  partial  density-of-states  tDOSI  analvsis  in  the  energy  range 
spanned  by  the  six  valence  and  conduction  bands  reveais  the  quantitative  details  of  the  bonding  mechanism 
between  the  Pd  and  H  constituents.  At  the  Fermi  energy,  the  high  Pd  d  to  H  s  DOS  ratio  -  10.3  is  found 
to  be  far  higher  than  expected  in  silver,  despite  the  fact  that  the  Fermi-vurface  geometry  is  similar.  The 
Held-induced  conduction-electron  spin  density  at  the  proton  site  is  evaluated  from  the  wave  functions  at  the 
Fermi  energy  The  calculated  value  of  the  spin-lattice  relaxation  rate  arising  from  the  contact  term  in  the 
hypertine  interaction  is  found  to  be  in  good  agreement  with  the  experimental  value  of  Wiley  et  aL 


L  INTRODUCTION  - 

The  metal- hydrogen  systems  have  received  wide 
attention  due,  in  part,  to  the  large  number  oi  their 
technical  applications. i,‘  They  are  also  important 
for  the  tests  they  provide  ot  several  theories.  The 
Study  of  the  behavior  of  a  proton  in  a  simple  metal 
was  started  with  the  pioneering  work  of  Friedel,3 
and  later  extended  by  several  authors.1  It  is  now 
well  established  that  the  screening  of  a  proton  in 
an  electron  gas  is  not  well  described  by  linear 
screening  theories,  because  the  perturbation  intro¬ 
duced  by  a  proton  is  simply  too  strong  and  nonlin¬ 
ear  screening  effects  are  very  important.1  These 
theories  essentially  treat  the  proton  as  an  isolated 
impurity  in  an  electron  gas,  and  are  based  upon  a 
■self-consistent  solution  of  the  Schrodinger  equa¬ 
tion  using  the  Hohenberg-Kohn-Sharr/  local-den¬ 
sity-functional  formalism;  they  are  expected  to 
give  interesting  results  only  for  the  so-called 
“simple  metals”  which  are  known  to  absorb  an  ex¬ 
tremely  small  quantity  of  hydrogen.  While  elec¬ 
tron-gas  models  provide  useful  information  about 
variation  in  the  screening  with  variation  of  the  av¬ 
erage  inter  electronic  distance,  their  quantitative 
predictions  must  be  treated  with  caution,  as  the 
Jellium  model  is  recognized’’  to  be  a  quite  crude  re¬ 
presentation  because  even  the  so-called  simple 
metals  show  substantial  departures  in  their  Fermi- 
surface  properties,  etc.,  from  that  of  an  electron 
gas. 

Since  the  transition  metals  like  Pd  can  absorb 
a  large  amount  of  hydrogen,  the  proton  cannot  be 
treated  as  a  dilute  impurity.  Thus,  for  these  d- 
electron  systems  one  would  not  expect  to  obtain 
meaningful  resuits  from  electron-gas  theories.  To 


describe  properly  the  electronic  structure  of  tran¬ 
sition-metal  hydrides,  alternative  methods,  like 
the  augmented-plane-wave  !APW)  band- structure 
method,  have  been  used  bv  Switendick’  and  other 
authors/'5  Theoretical  studies  have  also  been  made 
of  nonstoichicmetrtc  compounds,  using  either  the 
supercell  method'  (with  an  assumed  rigid  order  of 
the  vacancies)  or  the  virtual-crystal  approxi¬ 
mation.5  Recently,  the  electronic  structure  of  dis¬ 
ordered  substoichiometric  PdHx  compounds  has 
been  obtained  using  the  average- /-matrix  approxi¬ 
mation/0  and  the  coherent-potential  approxi¬ 
mation.11 

Among  the  transition- metal  hydrides,  the  PdHx 
system  is  of  particular  interest  because  (i)  the 
superconducting  transition  temperature12  Tc  in¬ 
creases  with  increasing  values  of  x  (for  ,ts  0.8)  ar.d 
(ii)  the  reverse  isotope  effect13  on  T ,  (observed 
when  H  is  replaced  by  deuterium)  is  a  challenging 
problem  for  both  theorists  and  experimentalists. 
Miller  and  Satterthwaite11  proposed  a  qualitative 
explanation  of  (ii)  based  upon  changes  in  the  elec¬ 
tronic  structure  which  arise  from  differences  in 
the  zero-point  motion  amplitudes  of  the  two  iso¬ 
topes,  leading  to  a  stronger  bonding  of  the  hydro¬ 
gen  atom  (compared  to  D)  to  the  palladium  atom. 

An  alternative  explanation,  given  by  Ganmily,15 
proposes  that  the  onharmonicity  due  to  a  larsre 
zero-point  motion  amplitude  of  the  protons  leads 
to  a  different  electron-phonon  coupling  and  to  a 
reduction  of  the  electron-phonon  coupling  con¬ 
stant  \  when  D  is  replaced  by  H.  This  theory  as¬ 
sumes  that  the  electronic  properties  of  PdH  and 
PdD  are  identical  (e.g.,  same  density  of  states  at 
the  Fermi  energy. . .).  The  importance  of  inhar¬ 
monic  effects  in  the  electron-phonon  interaction  is 
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Superconductivity  and  magnetism  in  the  new  group  of  rare-earth  ternarv  borides.  MRh4B4,  has  been  reported  recently 
by  Matthias  et  al.  Particularly  striking  has  been  their  observation  of  re-entrant  magnetism  accompanying  the  return  to  a 
normal  state  at  low  temperatures  t T  =  0.9  K. )  in  the  superconductor  ErRhaBa  which  has  a  superconducting  transition 
temperature  of  8.5  K.  We  discuss  the  orcm  of  superconductivity  and  magnetism  in  these  systems  using  results  of  ab  imrio 
self-consistent  LMTO  energy  band  calculations  for  the  alloys  M  =  Y,  Er  and  Ho. 


1.  Introduction 

There  has  developed,  in  recent  years,  great  interest 
in  the  electronic  structure  and  properties  of  ternary 
compounds.  A  key  motivation  has  been  the  hope  that 
higher  superconducting  transition  temperatures.  7*c, 
and  higher  critical  field  materials  could  be  obtained. 
For  the  Chevrel  phases,  with  general  formula 
MiMo6S3  (where  M  is  a  metal  and  x  a  number  near 
1),  there  are  a  number  of  metallic  elements  which 
cause  the  compound  to  be  superconducting.  Supris- 
ingly,  Fischer  and  collaborators  [  1  ]  at  Geneva  and 
Rennes  found  that  the  paramagnetic  rare-earth  metals 
also  formed  superconductors  whereas  earlier  the  addi¬ 
tion  of  paramagnetic  rare-earth  impurities  was  always 
found  to  lower  Tc.  Recently.  Matthias  et  al.  [2J 
reported  that  in  a  new  group  of  ternary  borides, 
mrh4b4,  either  ferromagnetism  (M  =  Gd.  Tb.  Dy, 

*  Supported  by  Air  Force  Office  of  Scientific  Research,  the 
National  Science  Foundation  and  the  Energy  Research  and 
Development  Administration. 


Ho)  or  superconductivity  (M  =  Y,  Nd,  Sm,  Th.  Er. 

Tm  and  Lu)  existed.  In  table  1  we  list  the  transition 
temperatures  of  some  of  the  rare-earth  compounds. 
We  see  a  sharp  break  from  magnetism  (T^1)  and 
superconductivity  (7^)  in  going  from  Ho  to  Er.  Of 
considerable  interest  has  been  the  later  discovery  ( 3  J 
of  re-entrant  magnetism  (at  T  =  0.9  K.)  accompanying 
a  transition  to  a  normal  metallic  state  in  ErRh4B4.  a 
rare-earth  superconductor  with  7^.  =  8.7  K,  this  has 
also  generated  considerable  excitement  as  to  the  onam 
of  this  new  phenomenon.  It  has  also  raised  imporiar.i 
questions  regarding  our  understanding  of  magnetism 
and  superconductivity  and  the  nature  of  their  possi¬ 
ble  coexistence. 

In  this  paper,  we  discuss  the  origin  of  these  phe¬ 
nomena  using  results  of  ab  initio  self-consistent 
energy  band  calculations  carried  out  on  three  of  these 
MRH4B4  systems,  with  M  =  Er.  Yfan  I  I K  super¬ 
conductor)  or  Ho  (a  ferromagnet  below  t>.fc>  K)  in 
their  paramagnetic  states.  Very  similar  total  and 
orbital  (/)  angular  momentum  decomposed  atomic 
densities  of  states  (DOS)  were  found  for  the  three 
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1.  Introduction 

Ch.  2  has  given  a  review  of  the  theory  of  the  isomer  shift  and  some  basic  fea¬ 
tures  of  electron  densities  as  determined  theoretically  for  free  atoms  or  ions. 

In  this  chapter,  we  describe  in  some  detail  the  present  state  of  our  theoretical 
abilities  to  determine  electron  densities  in  condensed  matter.  It  is  apparent  at 
once  that  numerous  approximations  have  to  be  made  to  reduce  a  manv-body, 
many-center  problem  to  a  tractable  form.  The  great  variety  of  more  or  less 
severe  approximations  which  have  been  applied  in  reducing  the  many  electron 
problem  has  led  to  a  literature  with  a  rich  vocabulary  and  confusing  terminol¬ 
ogy.  One  purpose  of  this  review  is  to  try  to  clarify  the  differences  and  similar¬ 
ities  between  various  methods,  and  to  assess  the  state  of  theoretical  and  com¬ 
putational  development  in  predicting  the  wavefunction-dependent  properties. 

In  order  to  explore  the  range  of  validity  of  several  models,  we  discuss  a  selec¬ 
ted  number  of  experimental  results  where  comparisons  can  be  made. 

We  begin  by  discussing  the  highly  successful  self-consisient-field  models  for 
treating  free  atoms  and  ions.  The  effects  of  spin  and  symmetry  restrictions  on 
the  single  determinant  wavefunction  3nd  approximations  to  the  Hartree-Fock 
exchange  operator  are  surveyed.  The  importance  of  relativistic  and  finite- 
nuclear-size  phenomena  for  determining  quantitative  results  for  electron  den¬ 
sities  near  the  nucleus  is  stressed. 

Free  atom/ion  models  have  an  obviously  limited  utility  in  molecules  and 
solids,  due  to  chemical  bonding  and  electron  delocalization.  Perturbed  atom, 
and  finally,  fully  variational  molecular  orbital  theories  represent  one  response 
to  these  limitations.  It  is  apparent  that  self-consistent  molecular  cluster  models 
can  (within  a  one-electron  framework)  represent  a  satisfactory  generalization 
of  the  free  atom/ion  scheme  for  localized  properties.  We  discuss  the  encourag¬ 
ing  progress  which  has  been  made  in  recent  years,  and  find  indications  that  first- 
principles  MO  calculations  may  soon  prove  to  be  more  useful  than  extrapolat¬ 
ing  atomic  results. 

The  delocalized  Bloch  wavefunction  approach  for  periodic  solids  has  now 
reached  a  fairly  high  state  of  development,  particularly  for  metals.  We  display 
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We  proem  the  results  of  a  first-principles  fully  self-consistent  study  of  the  electiontc  properties  of  cubic 
boron  nitride  in  the  local-density  formalism  using  our  previously  published  numencal-basis-set  linear 
combination  of  atomic  orbitals  scheme.  The  resulting  band  structure  shows  considerable  disagreement  with 
previously  published  ofthogonalized-plane-wave,  augmented-plane-wave.  and  pseudopotential  studies.  A 
detailed  study  of  the  ground-state  properties  of  the  system,  such  as  x-ray  scattering  factors,  cohesive  energy, 
equilibrium  lattice  constant,  and  their  behavior  under  pressm-  yields  very  good  agreement  with  avi..ao:e 
experimental  data.  Reasonably  good  agreement  is  obtained  for  excited  state  properties  determined  by  optical 
and  x-ray  absorption  measurements.  The  bonding  characteristics  in  this  prototype  of  III- V  compounds  are 
rfi«g«w»t  in  detail  and  compared  with  results  of  our  previous  study  of  its  isoelectronic  homopolar  analog, 
diamond,  and  with  studies  on  the  hexagonal  graphitehke  modification  of  BN. 


L  INTRODUCTION 

The  extraordinary  properties  of  cubic  boron 
nitride  have  been  the  subject  of  intense  experimen¬ 
tal  investigations  since  its  discovery  in  1957  by 
Wentorf.1  The  simplest  III- V  compound,  it  is  iso- 
electronic  and  isostructural  witn  diamond.  Its 
hardness  matches  that  of  diamond3  while  its  effec¬ 
tive  ionic  charge  is  the  highest  of  all  known  III-  V' 
compounds.3*4  Like  carbon,  BN  crystallizes  in 
two  forms,  a  cubic  zinc-blende  and  an  hexagonal 
(graphite-like)  structure.  These  properties  make 
a  comparative  study  of  these  compounds  particu¬ 
larly  Interesting.  We  have  previously  studied  the 
electronic  properties  of  diamond5,8  in  the  local- 
density  formalism  (LDF)7'8  using  our  first-prin¬ 
ciples  self-consistent  numerical-basis-set  linear 
combination  of  atomic  orbitals  (LGAO)  method. ’•  10 
In  particular,  we  have  investigated  the  band  struc¬ 
ture,  charge  distribution  (as  measured  by  x-ray 
scattering  factors),  cohesive  energy,  equilibrium 
lattice  constant,  behavior  under  pressure,  and 
the  bonding  characteristics  induced  by  the  various 
electron  exchange  and  correlation  potentials  per¬ 
taining  to  the  local-density  formalism.  In  this 
study,  we  extend  our  investigation  to  the  isoelec¬ 
tronic  cubic  boron  nitride  in  an  attempt  both  to 
understand  the  electronic  structure  features  in  the 
heteropolar  zinc-blende  system  and  to  test  further 
the  applicability  of  the  LDF  to  quantitative  de¬ 
scriptions  of  covalent-ionic  systems. 

The  crystal  structure  of  cubic  boron  nitride  has 
been  established  by  Wentorf1  (lattice  constant  a 
*3.615*0.001  Aat25JC).  Its  optical  properties 
in  the  infrared  region  have  been  studied  by  Gielisse 
*t  al*  and  recently  by  Chrenko11  who  has  also  investi¬ 
gated  its  ultraviolet  absorption  spectra.  Some  pre¬ 
liminary  results  have  also  been  published  by  Phillip 


and  Taft12  and  HalperinandKatzir.13  The  spectra  of 
the  core  levels  have  been  studied  by  means  of  x-ray 
emission  spectroscopy. 14  * 18  The  distribution  of 
electrons  in  the  ground  state  recently  has  been  studied 
by  Weiss17  by  means  of  x-ray  crystallography. 

Theoretical  studies13'37  of  the  electronic  struc¬ 
ture  of  cubic  BN  include  the  work  of  Kleinman  and 
Phillips18  using  a  pseudopotential  technique  in 
which  the  heteropolar  antisymmetric  potential 
(which  vanishes  in  the  homopolar  diamond  case),15 
is  treated  as  a  small  perturbation  on  the  previously 
computed  diamond  potential.  Hamstreet  and  Fong25 
have  indicated  that  since  the  basic  pseudopotential 
cancellation  in  crystals  like  BN  made  up  of  first- 
row  atoms  is  incomplete  due  to  the  lack  of  p  char¬ 
acter  in  the  core  states,  a  nonlocal  correction  to 
the  pseudopotential  is  necessary.  Adjusting  em¬ 
pirically  their  local  and  nonlocal  pseudopotential 
parameters,  they  obtained  a  band  structure  that 
differs  substantially  from  that  of  Kleinman  and 
Phillips.18  In  a  subsequent  or  thogonalized-plane- 
wave  (OPW)  study,  Bassani  and  Yoshimine23  have 
similarly  faced  a  plane-wave  convergence  difficul¬ 
ty  due  to  the  lack  of  pseudopotential  cancellation 
for  crystal  states  having  p  symmetry.  (Their 
model  employs  atomic  Hartree-Fock  charge  ten¬ 
sities  to  represent  the  crystalline  density  and  the 
crystal  exchange  potential  is  linearized  with  re¬ 
spect  to  the  single-site  potentials).  Wiff  and  Ke- 
own33  and  Keown34  have  calculated  the  band  struc¬ 
ture  of  BN  using  the  augmented-plane-wave  (APW) 
method  with  a  muffin-tin  crystal  potential  in  which 
the  boron  and  nitrogen  potentials  were  obtained  by 
linearly  scaling  the  carbon  potential  and  adding  a 
constant  Madelung  term  (with  effective  charce 
chosen  empirically  to  yield  a  gap  of  8.8  eW  to 
account  for  the  partial  ionic  character  of  the  sys¬ 
tem. 
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Alex  Zunger*  and  A.  J.  Freeman 

Department  of  Physics  and  Astronomy,  and  Materials  Research  Center,  Northwestern  University,  Evanston,  Illinois  60201 

(Received  9  February  1978) 

Self-comment  ab  initio  studies  of  the  electronic-energy-band  structure  of  cubic  CdS  are  reported  within 
the  local-denaity-functiornl  (LDF)  formalism.  All  electrons  are  included  using  our  previously  reported  linear- 
comhinatirai-of-atomic-orbitals  method  in  a  numerical  basis  representation.  In  the  tint  set  of  calculations  we 
employ  the  same  lattice  constant,  exchange  (only)  potential,  and  computational  parameten  as  were  used  by 
Soiled  etaL  in  their  early  self-consistent  orthogonalized-plane-wave  (SCOPW)  investigation  so  that  a  direct 
comparison  of  results  can  be  made  and  the  validity  of  the  SCOPW  approach  for  covalently  bonded  4 d 
systems  can  be  assessed.  In  the  second  set  of  calculations,  the  Stukel  et  aL  computational  restrictions  are 
relaxed,  a  more  accurate  lattice  parameter  is  employed,  ana  the  fCohn-onam  exchange  and  the  Singwi  et  aL 
correlation  potential  are  used  to  obtain  the  local-density  formalism  solutions  to  the  band  problem,  including 
variation  of  the  band  structure  and  related  prope.-ties  with  pressure  (change  of  lattice  constant).  Comparison 
with  optical  and  x-ray  and  uv  photoemission  experiments  for  excitations  of  both  the  s-p  and  metal  d  bands 
m  the  5-19  eV  region  indicate  very  good  agreement.  The  direct  gap  at  T  is.  however,  found  to  be  0.5  eV 
(25%)  too  small,  a  discrepancy  similar  to  that  previously  found  in  nonempmcal  studies  for  other  heteropolar 
insulators  (e.g..  Ne  and  LiF).  This  is  traced  to  the  neglect  of  the  different  orbital  relaxation  at  the  T,,  and 
T|  band  edges  and  to  the  noncancellation  of  the  seif-interaction  terms  characteristic  of  the  local-density 
potential.  Simple  atomic  total-energy  models  for  tnese  effects  are  shown  to  bnng  this  gap  into  good 
agreement  with  experiment.  It  is  concluded  that  a  first-pnnciples  (parameter-free)  exchange  and  correlation 
LDF  model  describes  very  well  the  main  electronic-structure  features  of  the  system. 


L  INTRODUCTION 

There  now  exists  vast  literature  comprising  ex¬ 
perimental  and  theoretical  investigations  of  the 
electronic  structure  and  related  properties- of  the 
important  II-VI  compounds.  Of  obvious  importance 
for  understanding  and  interpreting  many  of  the  ob¬ 
served  properties  of  these  materials  is  a  detailed 
knowledge  of  their  underlying  electronic  band 
structures.  In  this  regard,  recent  optical  and 
photoemission  experiments  have  played  an  impor¬ 
tant  role  in  elucidating  key  features  of  the  band 
structure  over  an  increasingly  large  range  of  en¬ 
ergies  above  the  fundamental  band  gap. 

On  the  theoretical  side,  pioneering  ab  initio  en¬ 
ergy-band  studies1"*  were  carried  out  for  many  of 
these  materials  using  the  non-muffin-tin  seLf-con- 
sistent  orthogonalized-plane-wave  (SCOPW)  meth¬ 
od.  This  technique  has  been  applied  in  recent 
years  to  the  study  of  the  electronic  structure  of  a 
large  number  of  semiconductors  and  semimetals, 
and  has  resulted  in  a  considerable  gain  in  insight 
into  the  systematics  of  both  the  LII-V  and  III-VI 
compounds  as  well  as  in  a  better  understanding  of 
the  photoemission,  high-pressure,  and  optical 
data  for  these  materials.  At  that  time,  the  only 
popular  alternative  method  was  the  non-self-con¬ 
sistent  empirical  pseudopotential  approach0'7  which 
completely  avoided  the  specification  of  the  micro¬ 
scopic  potential  in  the  solid  (i.e.,  the  nature  of  the 
Coulomb  exchange  and  correlation  fields)  by  pa¬ 


rametrizing  empirically  the  first  few  Fourier 
components  of  an  unspecified  crystal  potential. 
Although  this  phenomenological  approach  has  pro¬ 
duced  a  great  deal  of  information  on  the  variation 
of  the  gaps  in  II-VI  compounds  and  the  origin  of 
the  important  band -structure  critical  points  in  the 
Brillouin  zone,  no  meaningful  comparison  of  the 
adequacy  of  the  SCOPW  approach,  as  well  as  that 
of  the  underlying  exchange  models,  could  be  made. 
More-recent  developments  in  semiempirical  non¬ 
local  (non-self-consistent)  pseudopotential  theory3 
have  improved  considerably  the  agreement  between 
the  calculated  band  structure  and  x-ray  and  uv 
photoemission  data  by  introducing  a  few  additional 
disposable  parameters  to  describe  the  geometry 
of  the  angular -momentum -dependent  model  pseudo- 
potential  terms.  Although  a  successful  paramet- 
rization  scheme,  little  could  be  learned  regarding 
the  microscopic  nature  of  the  elementary  inter¬ 
actions  which  lead  to  this  success.  Further,  due 
to  the  independence  of  the  parametrized  model  po¬ 
tential  on  its  associated  variational  wave  function, 
questions  such  as  self-consistency  and  the  nature 
of  the  crystal  dissociation  products  [which  are 
well-defined  local-density-functional  (LDF)  atoms 
in  both  linear-combination-of-atomic -orbitals 
(LCAO)  and  SCOPW  models  j  have  remained  large¬ 
ly  unresolved  and  have  inhibited  a  meaningful  com¬ 
parison  with  first-principles  models  based  on  well- 
understood  (but  less  successful)  crystal  potentials. 
Subsequent  studies  on  the  convergence  properties 
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The  energy  band  structure  of  TiSe,,  determined  in  the  local-denjuy  approach  yields  a  semimetaJ  (band 
overiao  0.20  +  0.05  eV)  with  holes  at  T  and  electron  pockets  only  at  L.  The  dimensions  of  the  electron 
pocket  indicate  the  presence  of  (7-8)  x  10"  camers/cm'  in  excellent  agieement  witn  coin  transport  and 
angular-resolved  photoemission  data.  The  observed  charge-density  wave  is  attributed  to  characteristic 
“volume"  effects,  t.e.,  nesting  of  parallel  electron-hole  bands  at  £f  separated  by  the  T-L  zone-boundary 
wave  vector. 


The  nature  and  unusual  properties  of  the  elec¬ 
tronic  state  of  TiS,  and  TiSe,  have  become  the  sub¬ 
ject  of  intense  interest  and  study.  The  question  as 
to  whether  TiS2  (and  TiSe,)  is  a  semiconductor  or  a 
semimetal  is  still  not  resolved  and  the  closely  cou¬ 
pled  question  of  the  origin  of  the  unusual  Ti  tem¬ 
perature  dependence  of  its  resistivity  is  still  not 
understood.  The  recent  observation1* :  of  a  2 
x  2 c0  lattice  instability  in  TiSe,  at  temperatures 
below  202 JK,  coupled  with  its  semimetallic  con¬ 
ductivity  above  this  temperature  in  the  most  nearly 
stoichiometric  crystals  (~  1 0-°—  1 0-1  carriers/cm3), 
has  served  to  intensify  interest  and  speculation  as 
to  its  electronic  structure.  Of  particular  concern 
has  been  the  mechanism  with  which  to  explain  the 
origin  o f  its  charge-density  wave  (CDW)  and  the 
failure  to  observe  such  an  instability  in  TiS,.  The¬ 
ory  has  not  added  enlightenment.  The  early  non¬ 
self-consistent  Korringa-Kohn-Rostoker  (KKR) 
muffin-tin  band  calculation3  for  TiS,  (TiSe,)  gave  a 
fundamental  gap  of  2.0  eV  (1.2  eV)  at  T  and  a 
smaller  indirect  gap  of  1.4  eV  (0.5  eV)  between  T 
and  L;  similar  semiconducting  characteristics 
have  been  predicted  by  empirical  tight- binding 
studies.4 

Recently,  our  fully-self-consistent  numerical 
basis  set  linear  combination  of  atomic  orbitals 
(ICAO)  general  (non-muffin-tin)  potential  calcula¬ 
tion  for  TiSj,5  showed  good  agreement  with  optical 
properties  for  energies  below  16  eV.  Contrary  to 
previous  expectations,  a  small  indirect  gap  (0.2- 
0.3  eV)  was  found  to  occur  at  the  points  L  and  M  in 
the  BriUouin  zone,  with  a  larger  direct  gap  (0.8 
eV)  at  r.  (A  recent  comprehensive  analysis  of 
transport,  magnetic  susceptibility,  Hall  effect, 
resistivity,  Seebeck,  etc.,6  data  yielded  a  gap  of 
0.2-0.5  eV,  in  agreement  with  our  predictions.) 
Since  extrapolation  of  our  results  for  TiS,  to  the 
case  of  TiSe2  is  crude  at  best,  and  since  muffin- 


tin3  or  semiempiricai  models4  are  not  able  to  re¬ 
solve  these  questions,  we  undertook  an  ab  initio 
study  of  the  electronic  band  structure  of  TiSe,  in 
order  to  understand  its  various  measured  proper¬ 
ties.  Of  direct  interest  was  to  compare  with  the 
first  angle- resolved  photoemission  studies  of 
Bachracrt  et  al.,~  which  showed  conduction- banu 
overlap  of  the  valence  band  along  the  f-.U  direc¬ 
tion  by  at  least  0.5  eV. 

The  band  structure  was  calculated  in  the  local- 
density  functional  (LDF)  approach,  using  our  pre¬ 
viously  published  numerical  basis  set  LCAO  dis¬ 
crete  variational  method.5-8  We  find  that  TiSe,  is  a 
semimetal  with  an  indirect  negative  gap  (0.20  eV) 
between  r  and  L ;  contrary  to  some  previous  ex¬ 
pectations,2  no  electron  pocket  was  found  at  M. 

The  dimensions  of  the  electron  pocket  at  L 
(0.20 LH,  0.25 LA,  and  0.5L.V)  indicate  the  pre¬ 
sence  of  (7-3)  x  IQ20  conduction- band  carriers  ''cm3 
in  the  perfect  crystal.  This  is  in  very  good  agree¬ 
ment  with  the  independently  reported  transport1 
and  angle-resolved  photoemission  experiments.1 

In  our  calculations,  the  linear  basis  set  con¬ 
sisted  of  Ti  Is  to  4 p  orbitals  and  Se  Is  to  4 d  orbi¬ 
tals,  obtained  as  accurate  self-consistent  (numer¬ 
ical)  solutions  to  the  corresponding  sindie-site 
LDF  one-particle  equations.  About  3500-4000 
Diophantine  integration  points  per  unit  cell,  to¬ 
gether  with  a  cutoff  distance  of  27  a.u.  for  the  di¬ 
rect  space  lattice  sums,  were  needed  .as  was  the 
inclusion  of  the  Se  4 d  and  Ti  virtual  orbitals)  in 
order  to  obtain  a  numerical  accuracy  of  2-3  mRy 
in  the  band  eigenvalues  for  energies  below  (  ,+  1 
eV.  The  initial  crystal  potential  was  constructed 
from  an  overlapping  superposition  charge  density 
made  up  from  atomic  Ti  3d2 4s 2  and  Se  4s;4/>4  solu¬ 
tions  of  the  LDF  equations.  The  potential  included 
the  local  exchange  p‘/3  (with  a  2/3  coefficient)  and 
the  nongradient  correlation  contribution.10  All 
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We  show  that  observed  changes  in  the  nature  of  the  conducting  state  of  lT-TiSe-,  In  go¬ 
ing  from  the  normal  semimetallic  state  to  the  charge-density— wave  semiconducting  state 
can  be  successfully  modeled  by  variation  of  a  single  structural  parameter,  2  ,  which  mod¬ 
ulates  the  Ti-Se  bond  length.  These  ab  initio  band-structure  results  lead  to  a  number  of 
interesting  experimental  consequences. 


Current  interest  in  charge-density -wave  (CDW) 
phenomena  in  metals,  especially  the  layered  di- 
chalcogenides  of  the  transition  series,  remains 
high,  particular!,  as  regards  their  origin  and  the 
physical  properties  which  accompany  the  transi¬ 
tion  to  the  CDW  state.1'4  Recent  advances  in  ab 
initio  band-structure  theory  avoid  the  approxima¬ 
tions  underlying  the  more  conventional  approach¬ 
es  to  the  one-electron  problem  (e.g. .  muffin-tin 
approximations,  lack  of  self-consistencv,  short¬ 
comings  in  the  basis-set  representations,  etc.), 
and  have  been  found  to  describe  accurately  the 
properties  of  the  normal  state  of  these  compounds. 
In  particular,  the  nature  the  conducting  state 
o£  TiSj5’6  and  TiSe,7  (semiconductor  and  semi- 
metal,  respectively),  which  had  been  highly  con¬ 
troversial  for  some  time.8  has  now  been  correct¬ 
ly  predicted;  even  more  strikingly,  these  calcu¬ 
lations  were  able  to  predict  the  dimensions  of  the 
electron  and  hole  pockets  in  TiSe2,7  in  excellent 
agreement  both  with  transport  measurements7 
and  with  recent  angle-resolved  photoemission 
data.9  Similar  studies  on  TiS23  and  VSe3 10  were 
used  to  analyze  successfully  a  wealth  of  experi¬ 
mental  data,  including  x-ray  and  uv  photoemis¬ 
sion,  optical  reflectivity,  infrared  effective 
charge,  appearance-potential  spectra,  and  trans¬ 
port  properties,  and  have  provided  a  coherent 
microscopic  picture  of  the  chemical  bonding  in 
these  systems. 

Questions  pertaining  to  the  electronic  proper¬ 
ties  of  the  CDW  state,  however,  remain  largely 
unresolved.  Although  the  basic  phenomena  in¬ 
volving  phonon- mediated  interactions  between 
the  electron  and  hole  pieces  of  the  Fermi  surface 
(FS)  have  been  phenomenologically  understood 
for  some  time  on  the  basis  of  simple  momentum- 
conserving  geometrical  FS  models,  a  theoretical 
description  of  the  underlying  changes  in  the  elec¬ 
tronic  structure  of  the  CDW  state  is  still  lacking. 
This  has  a  simple  practical  reason.  Unlike  the 
unit  cell  of  the  normal  phase,  which  contains  a 


tractable  number  of  atoms,  the  CDW  state  is  ex¬ 
tremely  complex  (e.g.,  in  TiSe.  it  contains  24 
atoms  and  128  of  its  720  electrons  are  m  valence 
states'  and  so  cannot  be  treated  by  the  present 
generation  one-electron  methods  with  an  accura¬ 
cy  comparable  to  the  relevant  band  overlaps  and 
bond  gaps  (0. 1-0.2  eV).11  In  particular,  the 
lar^e  number  of  internal  structural  degrees  of 
freedom  in  the  CDW  phase,  makes  it  practically 
impossible  to  establish  a  direct  correlation  De¬ 
tween  the  structural  deformations  accompanying 
the  CDW  formation  and  the  pertinent  electronic 
properties. 

We  have  isolated  a  simple  structural  variable 
in  lT-TiSe;  that  allows  the  electronic  response 
to  the  observed  complex  pattern  of  l2rtice  distor¬ 
tions  in  the  CDW  state  to  be  modeled  in  a  simple 
and  direct  way.  This  simple  structural  degree 
of  freedom  is  the  internal  crystal  parameter33  z 
that  modulates  the  Ti-Se  bond  length  and  deter¬ 
mines  the  height  of  the  Se  atoms  above  the  metal- 
layer  planes.  The  simplicity  of  the  model  per¬ 
mits  a  detailed  study  of  the  electronic  response 
to  such  structural  deformations  and  provides  a 
direct  characterization  of  the  conducting  states 
of  the  material  in  terms  of  its  electronic  struc¬ 
ture.  .Among  other  results,  variations  of  this  pa¬ 
rameter  in  a  first-principles  energy-band  calcu¬ 
lation  are  found  to  transform  the  system  from  a 
semimetallic  state  to  a  semiconductinglike  state 
characteristic  of  the  CDW  phase.  The  results  of 
these  energy-band  studies  lead  to  a  number  of  in¬ 
teresting,  and  possibly  important,  experimental 
consequences,  including  the  effect  of  pressure  on 
the  conduction  properties  and  the  prediction  of 
the  absence  of  a  surface  superlattice  due  to  atom¬ 
ic  relaxation  at  the  surface. 

The  band  calculations  were  performed  with  use 
of  the  numerical-basis-set  linear  combination  of 
atomic  orbitals  (LCAO)  method  described  previ- 
ously5'12;  briefly,  it  consists  of  a  first-principles 
LCAO  approach  to  the  local-density -functional 
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The  calculated  surface  generalized  susceptibility,  determined  from  aft  initio  self-con- 
sl stent  thin-fllm  eaergy  bands  for  the  unreconstructed  phase  of  W(001),  displays  a  promi¬ 
nent  peak  at  A7  when  matrix  elements  and  local-field  corrections  are  included  and  sup- 
porta  proposed  charge-density-wave  interpretations  of  the  observed  reconstruction  to  the 
c(2x2)  phase.  These  results  are  consistent  with  the  Debe  and  King  parallel- shift  model 
and  the  suggestion  that  W(001)c<2*  2)-H  at  room  temperature  represents  an  impurity- 
stabilized  substrate  reconstruction. 


Surface-structural  phase  transitions  are  be¬ 
coming  a  subject  of  intense  experimental  and  theo¬ 
retical  investigation.  A  prototypical  example  is 
given  in  the  low-energy-electron-diffraction 
(LEED)  studies  of  the  tungsten1' 2  (001)  surface 
which  show  a  temperature-dependent  phase  tran¬ 
sition  to  a  (/2xv  2)R45°  or  c(2x2)  structure  when 
the  temperature  is  lowered  below  about  300  K. 
This  transformation  seems  to  be  of  second  order 
and  is  reversible  on  varying  only  the  tempera¬ 
ture.  [  A  similar  transition  has  been  observed  on 
Mo(OOi).*]  Recent  investigations1' 2  conclude  that 


chemisorbed  impurities  (hydrogen  in  particular) 
need  not  be  present  on  the  surface  when  the  c(2 
x2)  structure  is  observed,  thus  implying  that  the 
transition  is  characteristic  of  the  clean  surface. 

Felter,  Barker,  and  Estrup1  and  Debe  and  Kinar 
have  suggested  that  relatively  small  periodic  dis¬ 
tortions  could  account  for  the  phase  transition, 
and  they  have  noted  that  this  is  compatible  with  a 
charge-density-wave  (CDW)  mechanism  as  in  the 
layered  transition-metal  dichalcogemdes.  A 
LEED  intensity  analysis3  for  W(001)  at  100-140  K 
supports  a  model  proposed  by  Debe  and  King4  in- 
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A  fint-pnnciples  self-consistent  (non-mufTm-tin)  band-structure  calculations  on  IT-VSe,  is  used  to  anaivze 
some  of  the  transport  and  optical  properties.  The  material  appears  to  be  metallic  with  characteristic  overlaps 
between  the  metal  d -based  and  nonmetal  p -based  bands.  These  overlaps  are  found  to  be  sensitively 
modulated  both  bv  the  change  of  the  (anomalously  high t  c/a  ratio  and  by  the  sandwich  height  parameter  z 
(determining  the  V-Se  bond  length).  This  leads  to  a  number  of  interesting  predictions  regarding  the 
electronic  structural  changes  associated  with  the  charge-density-wave  instability. 


L  INTRODUCTION 

The  origin  of  charge-density  waves  (CDW)  and 
their  reLation  to  the  observed  anomalous  behavior 
of  the  layered  transition-metai  dichalcogenides 
continues  to  be  of  great  theoretical  and  experimen¬ 
tal  interest.  Of  the  various  dichalcogenides,  the 
IT  family  of  the  group-V2  compounds  may  be 
among  the  most  interesting  for  study  because  of 
the  wide  variety  of  phase  transitions  and  associated 
anomalous  properties  which  they  exhibit.1'2  For 
the  ir-TaSz  and  lT-TaSe,  systems  the  Korringa- 
Kohn-Rostoker  (KKR)  energy-band  calculations 
of  Myron  and  Freeman3  gave  Fermi  surfaces  which 
are  very  similar  for  both  systems,  having  cross 
sections  which  are  approximately  constant  In 
planes  that  are  perpendicular  to  the  z  axis  and 
which  can  be  nested  by  approximately  the  same 
wave  vector  parallel  to  the  I'M  direction.  Using 
the  energy -band  structure  and  three-dimensional 
Fermi  surface,  Myron  et  al .*  calculated  accurately 
the  generalized  susceptibility  \(q) ,  which  arises 
from  both  Ferml-surface  nesting  and  band-struc¬ 
ture  effects  for  bands  just  above  and  below  the 
Fermi  energy  E r  Maxima  were  found  at  those  q 
values  which  correspond  to  the  CDW  nesting  vec¬ 
tors  found  by  Wilson  et  al.1  Recently,  Woolley 
and  Wexler1  determined  the  band  structure  and 
Fermi  surface  of  tF-TaS,,  -TaSe2,  and  -VSe2 
using  the  Layer  method  for  band  calculations  in  a 
so-called  double  muffin-tin  approximation  to  the 
potent  iaL 

In  our  recent  local-density  self-consistent  (SC) 
non-muffin-tin  investigations  of  the  electronic 
structure  and  properties  of  the  group-IVB  ma¬ 
terials  lT-TiSj®  and  -TiSe,,7  we  found  profound 
changes  to  occur  relative  to  the  results  of  the 
early  non-self -consistent  muffin-tin  KKR  results 


obtained  previously.3  Whereas  the  muffin-tin  KKR 
results  gave  a  semiconductor  with  a  fundamental 
gap  of  2.0  and  1.2  eV  at  r  ana  a  smaller  indirect 
gap  of  1.4  and  ~0.5  eV  between  r  and  L  for  TiS, 
and  TiSe,,  respectively,  our  self-consistent  nu¬ 
merical  basis  linear  combination  of  atomic  or¬ 
bitals  (LCAO)  method  reducea  the  calculated  di¬ 
rect  gap  in  TiS,  to  0.3  eV  and  the  indirect  gap  be¬ 
tween  0.2  and  0.3  eV.  For  TiSe,,  our  SC-LCAQ 
results  showed  it  to  be  a  semimetal  with  an  in¬ 
direct  overlap  of  0.2  eV  between  F  and  L — in 
agreement  with  transport2  and  angle-resolved 
photoemission12  experiments.  These  results  ap¬ 
pear  to  confirm  the  validity  of  the  local  density 
functional  energy-band  approach  for  describing 
the  electronic  structure  of  the  layered  dichal¬ 
cogenides. 

In  addition  to  the  standard  energy -band  study  of 
TiSe,.  ’  e  also  investigated11  the  effect  of  internal 
structural  changes— specif icallv,  the  eifect  of 
varying  (i)  the  so-called  sandwich  height  param¬ 
eter  z,  which  determines  the  vertical  position^  of 
the  Se  atoms  above  and  below  the  metal  plane,  and 
(ii)  the  c/  a  ratio  on  the  electronic  structure.  We 
found  that  the  observed  changes  in  the  nature  of 
the  conducting  state  in  going  from  the  normal 
semimetallic  to  the  charge-density-wave  semi¬ 
conducting  state  can  be  successfully  modeled  by 
variation  of  a  single  structural  parameter  z. 
Among  other  results,  variations  of  this  parameter 
in  the  energy -band  calculations  are  found  to  trans¬ 
form  the  system  from  a  semtmetal  to  a  semicon¬ 
ducting-like  state  characteristic  of  the  CDW  phase 
and  to  lead  to  a  number  of  possibly  important  ex¬ 
perimental  consequences,  including  the  effect  of 
pressure  and  intercalation  on  the  conducting  prop¬ 
erties. 

As  a  consequence  of  this  wo»k  on  the  group-IVB 
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The  electronic  structure  of  the  high-temperature  rutile  phase  of  NbO,  is  studied  by  the  lineanzed- 
*agmented-plane-wave  method.  Potentials  constructed  by  superposition  of  neuiral-atom  and  ionic-charge 
densities  »re  used  to  explore  variability  of  the  electronic  band  structure.  A  rigid-band  scheme  is  shown  to 
accurately  describe  optical  absorption  of  the  rutile  phase  of  NbO.  stabilized  by  the  addition  of  20  at  ao  Ti  as 
measured  by  Raccah  et  al  Differences  between  the  band  results  for  rutile  NbO;  and  the  optical  absorption 
measurements  on  the  low-temperature  body-centered  tetragonal  phase  of  NbO,  are  attributed  to  band 
splittings  induced  by  lattice  distortion  which  occur  at  the  phase  transition.  The  static-electron  response 
function  x(<V  is  calculated  in  the  constant-matm-elements  approximation.  In  contrast  lo  ihe  case  of 
isoelectromc  VO,,  no  Fermi-surt'ace  nesting  features  are  observed,  and  y(q)  is  found  to  be  structureless  in 
the  vicinity  of  the  point  P  =  (1/4.  1/4,  1/2)  which  has  been  associated  with  a  possible  soft-mode  phonon 
instability  responsible  for  the  lattice  transformation. 


L  INTRODUCTION 

Niobium  dioxide  shows  a  structural  pnase  trans¬ 
formation  (rutile  to  body-centered  tairagonal) 
around  810°C  accompanied  by  rapid  changes  in 
resistivity  and  magnetic  susceptibility.  The  high- 
temperature  rutile  phase  is  metallic,  while  the 
distorted  phase  is  found  to  be  semiconducting.1"' 
The  transition  is  less  abrupt  than  in  the  isoelec - 
tronic  compound  VO,  and  has  been  characterized 
as  second  order  by  neutron-diffraction  studies.' 

It  is  possible  to  prepare  good  single  crystals  of 
NbOj  in  contrast  to  the  first  transition  series 
compounds,  so  this  system  provides  a  good  test¬ 
ing  ground  for  comparing  theories  of  electron- 
phonon  interactions  ana  electronically  driven 
transitions  with  experiment.  In  addition,  pseudo¬ 
binary  compounds  such  as  TijNty.,0,  are  found 
useful  in  studying  effects  of  d-electron  occupancy 
on  electrical,  optical,  magnetic,  and  structural 
properties.®'9 

Neutron  and  x-ray  studies  have  provided  accurate 
values  for  atomic  positions  in  both  phases  of  NbO,. 
The  low-temperature  distorted  phase  is  understood 
to  be  essentially  equivalent  to  a  dimerization  of 
the  metal-atom  chains  found  along  the  rutile  c 
axis,  analogous  to  the  famous  Peierls  distortion. 

In  the  case  of  VOJt  x-ray  diffuse  scattering 
measurements  by  Terauchi  and  Cohen10  confirm 
the  existence  of  soft  phonon  modes  centered  at  the 
point  R  -  of  the  Briliouin  zone  as  predicted 


by  a  band- structure  study.11  However,  the  mech¬ 
anism  of  the  phase  transition  in  NbO,  is  less 
clear5'12-13  and  neutron-scattering  experiments 
have  failed  to  detect  low-frequency  soft-mode 
phonon  branches15  particularly  at  the  wave  vector 
P  -  (i,  ?,  i)  thought  to  be  associated  with  the  static 
lattice  distortion. 

In  this  paper  we  study  the  electronic  band  struc¬ 
ture  and  the  generalized  susceptibility  x(q)  in  order 
to  assess  the  possible  role  of  electronic  instabil¬ 
ities  in  the  observed  structural  phase  transforma¬ 
tion.  Since,  as  is  well  known,  first-principles 
calculations  on  such  systems  of  rather  low  sym¬ 
metry  are  difficult,  some  calibration  against  ex¬ 
periment  must  be  made  prior  to  estimating  the 
validity  of  the  theoretically  determined  electronic 
response  functions.  Thus  it  is  especially  irr.cort- 
ant  for  our  work  that  Lu  et  al.1*  have  made  exten¬ 
sive  optical  reflectance  measurements  on  single- 
crystal  NbO,  and  TiJNb1.,02  and  have  derived  the 
associated  absorption  function  «,(w).  These  data 
provide  the  basis  for  a  detailed  comparison  be¬ 
tween  theoretical  band  structures  and  experiment. 

For  the  band  studies  we  use  the  recently  devel¬ 
oped  linearized-augmented-plane-wave  (LAPW) 
method  of  Koelling  and  Arbman.15  As  discussed 
in  Sec.  n,  this  method  has  proven  to  provide  an 
efficient  and  accurate  approach  for  complex 
structures.  Using  a  rigid-band  scheme,  the  re¬ 
sults  of  our  calculations  are  found  to  describe  ac¬ 
curately  the  optical  absorption14  of  the  rutile  phase 
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LOCAL  DENSITY  APPROACH  TO  BULK  AND  SURFACE  CHARGE  DENSITIES* 


A.  J.  Freeman 
Physics  Department 
Northvescam  University 
Evanston,  IL  60201 


I.  INTRODUCTION 

As  is  clear  from  some  of  the  ocf.e r  lectures  giver.  at  this  Ins¬ 
titute,  Che  current  popularity  of  energy  band  theory  stems  from  it3 
successful  application  to  the  study  of  increasingly  diverse  problems 
in  saLid  scare  physics.  Recent  new  sophis ticacad  exoerimencs  on 
both  traditional  materials  and  those  having  complex  crystallographic 
structures  have  demanded , '  however .  r.cc  only  theoretical  descriptions 
of  elgenva lue  phenomena  but  also  detailed  and  precise  wave  functions 
with  which  Co  determine  noc  only  cnarge  and  spin  densities  but  also 
Che  expectation  values  of  different  observable  operators.  Such  a 
demanding  test  of  the  predictions  of  or.e-e lectron  cneory  has  Che 
additional  virtue  of  permitting,  by  their  comparison  with  experiment, 
accurace  determinations  of  the  relative  magnitude  and  importance  of 
many-body  effects  in  real  solids.  Thus,  as  emphasised  by  Norman 
March,  there  has  develooed  considerable  incarasc  in  applying  the 
Hohenberg-Kohn-Sham  .1,2)  local  density  functions-  (LDF)  formalism, 
and  it3  recent  extension  as  a  local  soin  dans  i^lk  functions  1  i'LSDF) 
formalism  (3),  to  the  investigation  of  various  ground  state  proper¬ 
ties  of  3oIids  despite  the  usual  difficulties  of  solving  the  asso¬ 
ciated  one-particle  equation  characterised  by  a  multi-center  non- 
•pherical  potential. 

Applications  of  the  LDF  formalism  to  atoms  74,5)  and  molecules 
(6)  have  yielded  encouraging  resulcs.  Similar  3oplicacions  for 
solids  are  complicated  by  (i)  the  need  to  consider  both  the  short 
range  and  the  long  range  muLticenter  crystal  potential  having  non- 
spherical  components,  'll)  the  difficulties  in  obtaining  full  self- 
consistency  in  a  periodic  system  3nd(iii)  the  need  to  provide  a 
basis  sec  with  sufficient  variational  flexibility.  Hence,  cheorecical 
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Tile  observation  of  magnetism  and  superconductivity  in  both  the  ternary  rare-earth  rhodium  bondes 
(MRh.B.)  and  the  Chevrel  phase  temanes  M.Mo.S,  (with  x  close  to  1)  and  most  recently  (i)  the  discovery 
of  re-entrant  magnetism  accompanying  a  transition  to  a  normal  metallic  state  at  T<TJ  (the 
superconducting  transition  temperature)  and  (n)  the  coexistence  of  magnetism  and  superconductivity  have  ' 

generated  considerable  excitement  as  to  the  origin  of  these  phenomena.  We  discuss  the  origin  of 
,  magnetism  and  superconductivity  in  the  MRh.B,  compounds  and  re-entrant  magnetism  in  ErRh,B,  using 
the  results  of  ab  initio  self-consistent  LMTO  energy  band  calculations.  The  total  and  separate  I- 
decomposed  contributions  to  the  density  of  states  tDOS)  arising  from  the  two  M.  eight  Rh  and  eight  B 
atoms  per  unit  cell  are  used  to  estimate  their  various  contributions  to  magnetic  ordering  (via  the  4f-5d 
RJCKY  interaction)  and/or  superconductivity  (using  a  Gaspan-GyoriTy  modei  to  obtain  the  electronic 
connbutiotis  'o  A.  the  electron-phonon  coupling  paramer-r.  and  McMiilan  st  ong  coupling  theory  to 
obtain  TJ).  Comparisons  are  made  with  recent  experimental  results  and  several  experimental  speculations 
about  the  re-entrant  magnetism  state  are  presented  including  the  possibility  of  a  mixed  (ferromagnetic  and 
superconducting)  state. 

PACS  numbers:  74.70.Lp.  74.70.Ps.  75.30.Hx 


nmtoDUCTioN 

The  discovery  (1,2]  of  magnet  Is™  and  supercon¬ 
ductivity  In  rare-earch  ternary  compounds  has  brought 
together  the  study  of  these  seemingly  opposing  col¬ 
lective  phenomena.  In  particular,  the  observations 
of  re-entrant  magnetism  ln^  ErRh^B^  and  -  Ho^.jMojSg 
and  rase  racently  the  coexistence  of  antlferro- 
aagnatlsn  and  superconductivity  In  several  of  Che 
Chevrel  ternaries  have  raised  questions  regarding 
Chet*  phenomena  and  our  fundamental  understanding 
of  the  interactions  leading  to  magnetism  and 
aup«  t  co  aduc  t  iv  1  cy . 

In  the  ternary  boride  MRh^B^  compounds  (with 
H  “  e  rare-eerth  metal)  Macchias  et_  el.  [1]  found  chat 
either  f arromagneclsm  (M  »  Gd,  tb ,  Dv,  and  Eo)  or 
superconductivity  (M  -  Y,  Nd,  So,  Er,  Tm,  and  Lu) 
existed  et  temperatures  T  <  12  K.  Table  I  gives 
both  the  aagnetlc  transition  temperatures  (T?)  and 
the  auperconducclag  transition  temperatures  (T?). 

The  behavior  shown  In  Table  I  Is  in  3harp  contrast 
to  that  observed  In  Che  ternary  Chevrel  phase  compounds 
“rf06s8  and  KxMojSeg  where  each  rare-earth  compound 
ma  found  to  be  superconducting  except  for  Ce  and 
Ea,  The  sharp  breae  in  properties  between  the  H  ■  Ho 
and  Er  boride  compounds  Is  of  especially  great 
interest— with  magnetism  for  Ho  (Tj?  ■  6.56  K)  and 
anpareonductlvlcy  for  Er  (T*  *  8.7  K) ,  even  though 
the  effective  magnetic  moments  of  both  RE  Ions  (10.6  Ug 
for  Bo  and  9-6  Ug  for  Er)  differ  by  only  a  small 
taount  (-101).  Surprisingly,  ErRh/^B;,,  which 
becomes  superconducting  ac  8.7,  was  found  to  become 
aagnetlc  it  I  *  0.9  K  with  the  return  of  the  system 
to  e  normal  conducting  state  ("re-entrant  magnetism" 
in  e  superconductor). 

Be  have  been  studying  the  origin  of  these 
phenomena  using  results  of  ab  Initio  self-consistent 
band  structure  calculations  carried  out  on  five  of 
these  ternary  borides,  notably  M  -  Y  (which  has  no 
*f  electrons),  Gd  (with  the  largest  total  spin). 


Ho,  Er,  and  Lu  (which  has  a  closed  4f  shell).  We 
are  particularly  concerned  with  Che  competing 
mechanisms  for  producing  magnetism  and  supercon¬ 
ductivity  In  the  rare-earth  systems.  However,  only  a 
qualitative  understanding  Is  ac  present  possible 
because  our  band  calculations  are  for  the  paramagnetic 
states. 

The  magnetic  and  superconducting  interactions 
are  Included  In  the  sense  of  perturbation  theory. 

A  more  complete — and  much  core  dlf f Iculc — cheorv 
would  solve  che  basic  Hamiltonian  wlch  boch  magnecic 
and  superconducting  Interactions  Included  sU-ul- 
taneously.  Such  a  treatment  would  provide  a 
unified  view  of  the  phenomena  of  magnetism  and 
superconductivity.  Here  we  consider  noc  only 
phenomena  in  which  me  nagneclc  Interactions 
inhibit  superconduccvilty ,  but  also  chose  In  uhleh 
the  superconducting  ordering  prevencs  che  onset  of 
magnetism.  This  discussion  leads  co  the  suggestion 
of  a  possible  new  experiment  and  the  occurrenc  of  a 
mixed  (magnecism  and  superconductivity)  stace  for 
ErRh^Bi . 

METHOD 

Since  some  details  of  our  energy  band  cal¬ 
culations  [4]  have  been  given  elsewhere,  and  due  to 
space  limitations,  we  note  here  only  chac  che 
calculations  were  done  self-cons  Is cently  for  the 
full  18  atoms  per  unit  cell  using  the  Linear  Muffin 
Tin  Orbital  (LMTO)  and  related  methods  with  a 
potential  thac  contained  Che  Kedin-Lundqulsc 
treatment  for  exchange  and  correlation  In  the 
relaclvlsdc  Dirac  equadon  but  nmout  spla-orblc 
coupling . 

EESU12S 

Since  che  basic  energy  : <  . -ucturea  are 
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The  generalized  susceptibility.  x(<7),  in  Pd  and  Pt  for  q  along  the  { 1 00  j .  [  1 10 1.  [Ill],  and  [120]  directions  was  deter¬ 
mined  from  their  APW  and  R  APW  energy  Dand  structures,  respectively,  using  the  analytic  tetrahedron  linear  energy  scheme 
of  Rath  and  Freeman.  The  band  structures  were  previously  found  to  yield  Fermi  surface  radii,  temperature  dependencies  of 
the  static  magnetic  susceptibility.  \t  T).  resistivity,  and  a  nun  lattice  relaxation.  T)T.  in  very  rood  agreement  with  cm  er;- 
ment.  In  the  x(<7>  calculations,  we  used  2048  tetrahedra  in  1  48th  irreducible  8Z  and  the  energy  eigenvalues  tor  bands  4.  5. 
and  6  which  cross  the  Fermi  energy  as  titled  to  a  Fourier  series  representation.  The  mtraband  parts  of  \tq)  at  q  -  'j  for  "otn 
metals  are  found  to  agree  with  the  density  of  states  at  the  Fermi  energy  to  without  0.5%  Our  results  show  that  tne  domi¬ 
nant  contribution  to  xintra  wises  from  the  dominant  band  5  whose  "jungle-gym"  FS  has  strong  nesting  features;  the  mam 
peak  for  Pd  occurs  at  the  same  q  value  t=  0.65  it /a)  for  q  along  the  ( 0q0| .  |tj.  q,  0] .  and  [<?,  q.  q  |  directions.  The  locus  of 
this  main  peak  is  a  square  in  the  (0,  0,  1)  plane.  The  maximum  of  Xmtra  for  q  along  the  [  1 10|  and  [111)  directions  are 
23%  and  13%,  respectively,  higher  than  the  value  of  \lq)  at  q  =  0.  For  q  along  the  [ 0 1 0 1  and  1 1  201  directions,  the  peak  is. 
however,  lower  than  the  value  of  x,ntra  4t  q  =  0.  Hence,  while  phonon  anomalies  are  predicted  for  the  [110]  and  j  1 1 1 1 
directions,  no  anomaly  is  predicted  for  either  the  ( 1001  or  [120]  direction.  The  predicted^  value  for  the  [  1 1 0 ]  anomaly, 
q  *  0.65  ir/a  is  close  to  the  experimental  value  of  ~0.7  ma.  Although  there  may  be  a  hint  of  an  anomaly  at  0.56  [111]  in 
the  measurements,  a  more  detailed  investigation  of  this  region  is  called  for.  For  platinum,  xintra  ior  <7  along  the  [010], 
[110|  and  [1U|  directions  has  mam  peaks  which  occur  at  q  =  0.68  ma.  0.75  ma.  and  0.85  nia,  respectively.  Here  too.  this 
main  peak  comes  from  the  nesting  of  the  jungie-gy  m  Firmi  surface  which  is  not,  however,  as  tlat  as  that  of  palladium. 
Anomalies  are  predicted  (although  weaker  in  Pt  than  in  Pd)  along  [110]  and  [111]  but  not  along  [  100]  and  1120].  The 
[110]  anomaly  is  lose  to  the  measured  q  value  ( -0.7-0.8  n,a).  Also  in  agreement  with  experiment,  we  predict  a  weaker 
[110]  anomaly  for  Pt  than  for  Pd.  In  both  Pd  and  Pt.  weaker  anomalies  are  predicted  for  the  ( 1 1 1 1  direction  than  for  the 
[110]  direction. 


1.  Introduction 

The  role  of  electronic  structure  in  determining 
phonon  anomalies  in  metals  has  become  a  subject  of 
intense  interest.  Strong  impetus  for  these  investiga¬ 
tions  has  been  given  by  the  correlation  of  spin  den¬ 

*  Based  on  work  performed  under  the  auspices  of  the  US 
Air  Force  Office  of  Scientific  Research,  and  National 
Science  Foundation  and  the  Department  of  Energy. 

**  Present  address.  Physics  Department.  Case  Western 
Reserve  University,  Cleveland.  OH  44106,  USA. 

•**  Science  Applications  Inc.,  Schaumberg,  11  60195,  USA. 


sity  waves,  charge  density  waves  (CDW)  and  phonon 
anomalies  with  maxima  in  the  generalized  suscepti¬ 
bility,  x(fl).  f°r  the  non-interacting  conduction  elec¬ 
trons.  Some  time  ago,  Kohn  [l]  showed  that  for  a 
free-electron  gas,  the  existence  of  a  sharp  Fermi  sur¬ 
face  leads  to  non-analyticities  in  the  phonon  disper¬ 
sion  relations  (now  known  as  “Kohn  anomalies  ’) 
which  occur  at  phonon  momentum  values  equal  to 
twice  the  momentum  of  the  electrons  at  the  Fermi 
surface.  Subsequently,  Taylor  f 2 J  treated  tiie  effect 
of  various  radii  of  curvature  of  a  single  band  Fermi 
surface,  Afanasev  and  Kagan  [3|  gave  a  more  detailed 
discussion  for  the  case  of  a  Hat  Fermi  surface,  and 
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Resume.  —  Utilisant  une  version  relativiste  de  la  methods  LAPW.  nous  presentons  des  resultats  d'une  etude 
detaillee  de  la  structure  de  bande  eiectromque  ainsi  que  de  quelques  proprietes  physiques  de  l'uranium  Jans  la 
phase  orthorhombique.  Les  densites  d'etats  totales  et  partielles  sont  utilisees  pour  discuter  les  resultats  de 
photo-emission  de  Veal  et  Lam.  Les  surfaces  de  Fermi  des  bandes  predominantes  6  et  7  presentent  une 
symetrie  elevee  ainsi  que  quelques  particularites  topoiojnques  t  nesting )  qui  peuvent  etre  comparees  au\ 
mesures  actuelles  d'effet  dHvA  d’Arko  et  Schirber.  Le  calcui  du  facteur  de  forme  neutromque  en  presence 
d’un  champ  magneuque  est  en  bon  accord  avec  les  mesures  de  .Ylagiic  et  al. 

Abstract.  —  We  present  some  results  of  a  detailed  study  of  the  electronic  band  structure  and  some  physical 
properties  of  orthorhombic  uranium  using  a  relativistic  version  of  the  LAPW'  method.  Total  and  projected  iby 
orbital  angular  momentum)  densities  of  states  are  used  to  discuss  the  photo-emission  data  of  Veal  and  Lam. 
The  Fermi  surfaces  of  the  dominant  6th  and  7th  bands  are  found  to  have  high  s>mmeiry  (and  some  nesting 
features)  of  interest  to  dHvA  measurements  being  pursued  by  Arko  and  Schirber.  The  theoretical  magnetic 
field  induced  neutron  form  factor  is  found  to  be  in  good  agreement  with  the  measurements  of  Maclic  et  al. 


The  anomalous  physical  properties  of  a-U  have 
made  it  a  subject  of  great  interest  for  experimental 
and  theoretical  study.  Unfortunately,  its  complex 
structure  (orthorhombic  with  4  atoms/unit  cell)  and 
the  resulting  lack  of  symmetry  have  made  its  theore¬ 
tical  study  by  conventional  energy  band  methods 
exceptionally  difficult  and  costly.  Thus,  while  we 
were  able  to  study  the  light  and  heavy  actinide 
metals  [1]  in  their  cubic  (high  temperature)  structu¬ 
res,  the  study  of  a-U  awaited  development  of  a 
computational  scheme  which  was  not  only  rapid  and 
efficient  but  which  avoided  the  so-called  asymptote 
problem  which  plagued  our  earlier  efforts. 

In  this  paper  we  report  some  results  of  our  relati¬ 
vistic  energy  band  studies  on  a-U  using  a  relativistic 
version  of  the  linearized  augmented  plane  wave 
method  [2-3]  (LAPW).  This  scheme  has  been  suc¬ 
cessfully  applied  to  the  study  of  some  actinide 
compounds,  notably  UGe,  [4],  Ulr,  [2],  and 
URh,  [5].  We  have  determined  the  energy  band 
structure,  density  of  states  (DOS),  orbital  angular 
momentum  projected  DOS  (for  use  in  analysing  the 
XPS  experiments  of  Veal  and  Lam  [6],  detailed 
Fermi  surface  cross-sections  in  close  collaboration 
with  the  work  of  Arko  and  Schirber  [7]  reported 
elsewhere  at  this  conference),  wavef unctions  and 
magnetic  field  induced  spin  densities,  neutron  ma¬ 
gnetic  form  factors  (in  close  collaboration  with  and 
to  understand  the  measurements  of  Maglic  et 

(*)  Supported  by  the  U  S.  NSF,  the  AFOSR.  and  the  DOE. 


al.  [8]).  and  generalized  susceptibilities,  ,*(q)  (for 
investigating  possible  electronically  driven  phonon 
anomalies  and  charge  density  waves).  Because  of 
space  limitations  we  are  able  only  to  give  a  brief 
indication  of  this  extensive  work  here-  A  full  report 
is  being  prepared  for  publication  elsewhere. 

The  three  lattice  parameters  and  the  atomic  posi¬ 
tional  parameter,  u,  were  obtained  by  extrapolating 
the  temperature  dependent  X-ray  measurements  of 
Barrett  et  cJ  [9]  to  the  temperatures  assuming  that 
the  system  was  cooled  under  pressure  in  order  to 
avoid  the  anomalies  associated  with  transitions  at 
T  <  43  K.  Using  these  parameters  (a  =  5.360  A. 
b  =  1 1.086  A.  c  =  9.328  A.  and  /x  =  0. 102)  our  mo¬ 
del  warped  muffin  tin  potential  was  constructed  in 
the  standard  overlapping  charge  density  model  using 
the  Kohn-Sham-Gaspar  (a  =  2/3)  exchange.  In  our 
calculation  we  use  the  f’  d:  s'  configuration.  Diffe¬ 
rent  estimates  including  a  pseudo-self-consistent 
model  (as  was  used  earlier  for  UGe,)  to  determine  an 
improved  choice  of  configuration,  indicate  that  a 
configuration  in  which  roughly  0.25  to  0.5  electron 
are  transferred  from  the  f  to  the  d  levels  is  more 
appropriate.  The  calculated  band  structure  was  fit¬ 
ted  with  a  Fourier  series  and  used  in  our  determina¬ 
tions  of  the  Fermi  surface.  DOS  and  y(<i>. 

Figure  1  shows  the  l  =  3  angular  momentum  pro¬ 
jected  DOS.  We  see  a  good  deal  of  structure  in  the 
total  DOS  (not  shown)  arising  from  the  lower  sym¬ 
metry  and  the  hybridized  -et  of  s.  p.  d.  and  f  bands. 
There  is  a  rapidly  increa'  fig  total  DOS  just  above  EF 
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The  electronic  and  magnetic  properties  of  a  coherent  modulated  structure  consisting  of 
three  atomic  layers  each  of  Cu  and  Ni  perpendicular  to  lllll  were  determined  from 
ab  initio  self-consistent  spin-polarized  energy-band  calculations.  We  predict  a  substan¬ 
tial  reduction  in  spin  magnetization  and  contact  hyperfine  fields  it,  at  the  Ni  layers  (in 
contrast  with  enhancements  reported  earlier  from  ferromagnetic  resonance  measure¬ 
ments),  charge  <0.1*  ")  and  spin  (~0.01mb)  transfer  to  the  Cu  sites,  and  large  at  the 
Cu  nuclei. 


PACS  numbers:  71.25.Pi,  73.40  Jn,  75.10. Lp 

Considerable  Interest  has  been  aroused  by  the 
observation  of  greatly  enhanced  elastic  proper¬ 
ties  of  short-wavelength  composition- modulated 
bimetallic  films1  and  especially  by  the  significant¬ 
ly  larger  magnetization  density  of  a  short-wave¬ 
length  (30  A)  CuNi  thin  composition- modulated 
film2  compared  to  that  of  pure  M  as  deduced  from 
ferromagnetic  resonance  experiments.  White 
and  Herring*  have  proposed  that  the  anomalously 
large  magnetization  value  is  indicative  of  a  com¬ 
plex  magnetization  distribution,  including  nega¬ 
tive  regions,  which  can  substantially  increase 
ferromagnetic  resonance  frequencies. 

In  this  Letter,  we  present  the  first  theoretical 
determination  of  the  electronic  structure  and  mag¬ 
netic  properties  of  a  realistic  layered  coherent 
modulated  structure  (CMS).  Ab  initio  self-con¬ 
sistent,  semirelaiivistic,  spin-polarized  linear 
muffin-tin  orbital  (LMTO)  energy-band  calcula¬ 
tions  have  been  carried  out  on  a  50-50-composi- 
tion  CMS  consisting  of  three  atomic  layers  each 
of  Cu  and  NI  modulated  along  the  l  111]  direction.* 
Ws  find  no  evidence  for  a  complex  magnetization 
distribution.  Instead,  the  spin  magnetization  of 
the  NI  layers  is  reduced  relative  to  that  of  pure 
NL  The  central  M  layer  moment  is  0.5Qmb 
whereas  that  of  the  interface  Ni  layers  is  0.37 
The  large  reduction  of  the  "interface"  Ni  layer 
moments  is  consistent  with  the  magnetically  dead 
layers  reported  earlier  for  Ni  layers  deposited 
on  a  Cu  substrate.4  The  calculated  core-polariza¬ 
tion  hyperfine  fields,  H„,  at  the  Ni  sites  also 
show  a  reduction  proportional  to  that  found  for  the 
spin  magnetization.  The  small  Cu  magnetic  mo¬ 
ments  (-0.01 /i8)  result  in  large  Hn  values  which 


are  substantially  larger  than  any  Knight  shifts  ob¬ 
served  for  Cu. 

Our  CMS,  shown  in  Fig.  1,  has  six  atoms  per 
unit  cell  and  four  inequivalent  layers  denoted  as 
"bulk"-like  for  the  central  Cu  and  Ni  and  “inter¬ 
face”  for  the  remaining  Cu  and  Ni.  Since  deposi¬ 
tions  of  a  metal  A  on  a  metal  B  result  in  a  lattice 
constant  for  the  initial  layer  which  is  between 
that  of  A  and  B ,  we  ferm  our  CuM  structure  with 


FIG.  1.  CuNi  CMS  modulated  along  l 111)  (without 
local  strain).  At,  A;,  and  A3  define  the  six-atom  unit 
cell;  >  and  b  denote  “interface”  and  “bulk.” 
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Abstract.-  we  describe  and  discuss  several  recent 
this  rapidly  developing  area  of  research  which  i3 
tpecezoscopis  e . 

The  study  of  surface  phenomena  has  added  a 
new  and  important  area  of  research  which  has  yet 
to  be  exploited  fully  by  Che  Mossbauer  spectro¬ 
scopic  t.  As  discussed  at  the  Conference,  this  is 
particularly  true  as  regards  transition  metals 
where  recenc  advances  in  experimental  and  theore¬ 
tical  methods  for  studying  bulk  and  surface 
properties  have  Increased  our  knowledge  about 
these  phenomena  buc  have  also  raised  important 
questions  about  our  fund,  -stal  unders tend ir~  in 
s  dramatic  way.  Here  1  will  mention  tvo  areas, 
surface  magnetism  and  surface  phase  transitions, 
sdiich  are  particularly  Important  examples  of 
rsceat  developments  in  this  fast  growing  field. 

A.  Surface  States.  Surface  Magnetization  and 


A  striking  example  of  some  of  the  conflicting 
results  Is  the  contradictory  conclusions  about 
the  validity  of  the  Stoner-Wi  V.farth-S later  band 
theory  of  fermaagnetlsm  drawn  from  electron  spin 
polarization  experiments  Cl],  The  possibly 
important  role  of  surface  phenomena  In  reconcil¬ 
ing  these  difficulties  has  been  recognized  in 
recent  years  and  theoretical  efforta  at  describing 
the  electronic  structure  of  surfaces  have  intensi¬ 
fied,  These  investigations  share  the  common 
problem  with  bulk  studies  of  treating  localized 
d  electrons  along  with  the  Itinerant  s-p  electrons. 
For  bulk  ays terns ,  considerable  progress  hss  been 
mads  in  the  last  few  years  in  this  direction  [2]. 


developments  in  the  study  of  surface  phenomena  in 
becoming  of  great  interest  to  the  Mds3bauer 

Similarly,  sophisticated  theoretical  methods 
are  being  developed  for  accurately  determining  the 
electronic  structure  of  transition  metal  surfaces. 
These  use  the  bulk  energy  band  methods  In  the  local 
density  formalism  as  extended  to  the  case  of  a 
thin  film  [3-15],  As  in  bulk  systems,  a  number  of 
thin  film  studies  for  transition  metals  has  shown 
that  self-consistency  plays  an  essential  role 
[6-12,13]. 

I  d«*cuss  in  detail  the  work  which  we  [9-12] 
have  carried  out  on  ferromagnetic  transition  metal 
films  using  our  recent  extension  of  the  lo.eal  spin 
density  functional  formalism  to  thin  films.  These 
are  the  first  spin  polarized  ab  initio  self-consis- 
tent  energy  band  studies  of  ferromagnetic  transi¬ 
tion  metal  films-  that  are  thick  enough  [9- layers 
for  Hl(QOl)  and  7  layers  for  Fe(001)]  to  accurately 
determine  the  energy  dispersion  and  spatial  charac¬ 
ter  of  the  surface  states  and  their  effects  on  the 
surface  spin  polarization  layer  by  layer  spin  and 
charge  distributions  and  layer  projected  density 
of  states. 

One  striking  feature  of  the  results  la  the 
Friedel  type  oscillation  Induced  by  the  surface 
discontinuity  in  the  layer  spin  magnetic  density  - 
but  not  in  the  charge  density.  A  result  of 
particular  interest  here  is  the  prediction  of  a 
decrease  (by  2(71)  in  the  surface  layer  magnetic 
moment  in  Ml  and  a  larger  Increase  in  the  case  of 
Fa.  There  is  no  evidence  for  magnetically  'dead' 
layers . 
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W«  report  detailed  results  of  ab-imuo  self-consistent  LMTO  energy  band  studies  of  several  Chevrel  phase  temanes.  EuMo^Sj 
and  SnMogS,.  In  these  unique  calculations,  all  electrons  of  the  full  15  atoms  /  unit  cell  are  treated  explicitly  with  a  self-consistent 
Hedin-Luadqvist  exchange  and  correlation  potential  in  the  relativistic  Dirac  equation  but  without  spin-orbit  coupling.  As  in 
oar  earlier  work  on  the  rare-earth  ternary  borides,  the  total  and  separate  tbv  atom)  /-decomposed  contributions  to  the  density  of 
slates  (which  show  the  Eu  atoms  to  be  magnetically  isolated)  are  used  to  discuss  their  magnetic  and  superconducting  properties. 


The  discovery  that  the  ternary  compounds  of 
the  rhodium  borides.  MRh4B4,  and  the  Chevrel 
phase  molybodenum  chalcogemdes.  MMo6Ss  (or 
Ses),  exhibit  (i)  superconductivity  and/or  mag¬ 
netism  including  (ti)  the  co-existence  of  the  two 
and  (iii)  re-entrant  magnetism  in  several  important 
systems  has  generated  great  interest  in  the  origin  of 
these  phenomena  both  experimentally  and  theoreti¬ 
cally.  Experiments  on  these  systems  have  far  out¬ 
stripped  theoretical  understanding  and  have  gener¬ 
ated  numerous  unsolved  questions. 

In  the  case  of  the  rare-earth  ternary  rhodium 
borides  our  self-consistent  LMTO  energy  band 
studies  have  provided  a  qua|it?»'v“  understanding 
of  the  underlying  electronic  structure  and  some  of 
the  basic  phenomena  observed  [1],  These  results 
were  also  used  to  predict  [1,  2]  a  possible  mixed 
state  in  ErRh4B4  at  temperatures  above  the  re-en¬ 
trant  magnetism  temperature — a  state  which  may 
have  been  seen  in  the  magnetization  vs.  magnetic 
field  measurements  of  Ott  et  al.  [3].  The  complexity 
of  the  crystal  structure  of  the  Chevrel  phase 
ternaries  have  made  difficult  even  non-self-con¬ 
sistent  energy  band  studies.  Instead,  early  studies 
were  carried  out  with  approximate  molecular 
cluster-tight  binding  [4]  and  (limited  basis  set) 
localized  orbital  methods  [5],  We  here  report  re¬ 
sults  ab  initio  self-consistent  LMTO  energy  band 
studies  of  several  of  the  MMosX,  compounds  (M 
«■  Eu  and  Sn)  which  include  all  electrons  and  all 
15  atoms  in  the  unit  ceil. 

The  band  calculations  have  been  performed  by 
use  of  the  Linear  Muffin  Tin  Orbital  (LMTO) 
method  [6].  The  scheme  of  calculation  is  essentially 
the  same  as  has  been  used  for  calculations  on  A 15 
compounds  (7]  and  ternary  borides  (1). 

The  crystal  structure  data  of  Marezio  et  al.  [8] 
on  PbMo4S,  has  been  used  in  all  our  Chevrel  phase 
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calculations,  but  without  distortion  of  the  Mo- 
octahedron  and  with  a  rhombohedral  angle  of  40:. 
Further,  if  one  considers  spherical  symmetric 
potentials  around  each  site,  this  leaves  only  4  types 
of  atoms  in  the  structure,  whereas  including  the 
distortion  would  make  most  of  the  15  atoms  per 
cell  inequivalnet.  In  this  simplified  structure  there 
is  one  Sn  site,  six  equivalent  Mo  sites,  but  two 
types  of  S-sites. 

The  potential  used  the  Hedin-Lundqvist  treat¬ 
ment  of  exchange  and  correlation.  The  charge  den¬ 
sities  from  the  core  states  as  well  as  from  partly 
occupied  f-states  were  recalculated  in  each  itera¬ 
tion  by  using  the  actual  MT-potential.  The  basis 
set  included  s.  p  and  d  states  for  all  atoms,  while  in 
the  three-centre  terms  in  the  LMTO  matrices  /w 
<.  3.  This  results  in  an  eigenvalue  problem  of  di¬ 
mension  135  x  135.  The  sharp  structure  in  the 
total  DOS  shown  in  fig.  1  for  EuMo6Ss  arises  from 
the  very  flat  nature  of  the  energy  band  structure. 
Very  similar  results  were  also  found  for  SnMo6S3 
which  is  a  superconductor.  These  flat  bands  yield 
high  electron  velocities  at  £F  and  appear  to  explain 
the  unusually  high  upper  critical  fields.  H,:.  found 
in  these  superconductors.  We  see  that  there  is 
considerable  structure  in  the  total  DOS.  particu¬ 
larly  around  £F.  which  arises  from  the  Mo-4d 
electrons.  From  a  partial  DOS  calculation  [cf.  table 
1],  we  find  that  there  is  a  high  4d  DOS  at  £F  which 
is  favorable  for  superconductivity.  There  is  a  dis¬ 
tinct  gap  in  the  DOS  just  above  £F  which  falls  in 
the  middle  of  the  Mo-d  bands  and  a  smaller  band 
gap  between  the  Mo-d  and  S-p  states  as  was  also 
inferred  or  seen  in  the  earlier  calculations  [4.  5j.  As 
for  most  other  high  £c  compounds,  the  Fermi  level 
falls  just  2  electrons  below  the  middle  of  the  bond¬ 
ing-antibonding  “gap'"  of  the  metal  d-states.  (This 
explains  why  the  temanes  Mo;RejSs,  which  have  2 
more  electrons  per  cluster,  are  semiconductors.) 
Usually  for  Mo  compounds,  the  d  bands  are  oc- 
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self*consistent  band  calculations 

T.  Jarlborg 

Physics  Department.  Northwestern  University.  Evanston.  Illinois  60201 

% 

A.  J.  Freeman  , 

Physics  Department.  Northwestern  Umversitv.  Evanston.  Illinois  60201 
and  Argonne  National  Laboratory.  Ary  mine.  Illinois  604 JO 
(Received  20  July  1979.  revised  manuscript  received  27  February  1980) 

AH-atom-all-electron  self-consistent  semirelativishi;  linear  muffin-tin  orbital  energy-band 
studies  are  reported  for  TiBe;.  ZrZn,  (at  ambient  and  high  pressures),  the  high- 7"^  superconduc¬ 
tor  ZrVj,  and  the  ultralow -Tc  superconductor  YAK  Total  and  partial  t by  atom  type  and  / 
values)  density  of  states.  Stoner-like  parameters,  electron-phonon  coupling  parameter  h  and  su¬ 
perconducting  transition  temperatures  Tc  (determined  within  the  rigid  muffin-tin  approximation 
and  McMillan's  strong  coupling  theory  with  and  without  paramagnon  contributions),  and  their 
behavior  under  applied  pressure  are  used  to  discuss  the  origin  of  their  observed  magnetism 
and/or  superconductivity.  From  an  analysis  of  the  results,  we  suggest  possible  high-fr  super¬ 
conductivity  at  high  pressure  for  ZrZn;.  and  its  unlikely  occurrence  for  TiBe-  unless  dominant 
soft  phonon  modes  exist. 


I.  INTRODUCTION 

The  cubic  Laves  or  Cl  5  compounds  of  type  MgCu2 
exhibit  a  number  of  interesting  phenomena  including 
(i)  the  much  studied  v  eak  ferromagnetism1-9  in 
ZrZn2  and  (ii)  correlations  between  high  supercon¬ 
ducting  transition  temperatures  Tc  and  lattice  instabil¬ 
ities1*^11  ZrV2  and  HfV,,  which  are  well  known  in  the 
more  familiar  A 1 5  and  rocksalt  transition-metal  com¬ 
pounds.  Interest  in  superconductivity  and  itinerant 
magnetism  in  the  C15's  has  recently  been  revitalized 
foiiowing  the  report  of  itinerant  antiferromagnetism 
in  TiBe2  by  Matthias  et  at.  13  This  discovery  followed 
the  expectation  raised  by  the  theoretical  proposal  of 
Enz  and  Matthias19  that  the  itinerant  ferromagnetism 
ofZrZn2  arises  from  inhibited  "p-state"  pairing.  In 
this  view,  the  Cooper  pair  interaction  is  thought  to  be 
repulsive  at  short  distances  (j  state)  due  to  a  dom¬ 
inant  soft  phonon,  leaving  only  p  and  high  pair  states 
to  be  attractive.  Thus,  the  s-state  repulsion  would 
contribute  to  the  Stoner  factor  sufficiently  to  tip  the 
balance  toward  magnetism  and  to  prevent  the  oc¬ 
currence  of  p-state  superconductivity.  They  further 
predict  that  above  a  critical  pressure  at  which  the 
magnetism  is  destroyed  in  ZrZn2,  p-state  pairing  su¬ 
perconductivity  would  appear  and  in  TiBe2  would 
coexist  with  any  antiferromagnetic  ordering. 

Although  extensively  studied  experimentally,  only 
approximate  model  calculations3  9  have  been  per¬ 
formed  for  the  electronic  band  structure  of  ZrZn2 
treated  as  a  diamond  lattice  of  zirconium  atoms  with 
the  zinc  atoms  entirely  neglected.  No  energy-band 
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studies  have  been  made  of  TiBe2  (or  the  high-Tr  su¬ 
perconductors).  Thus,  it  has  not  been  possible  to  as¬ 
sess  the  possible  electronic  contributions  to  either  a 
more  conventional  approach  to  superconductivity 
and/or  magnetism  or  to  the  inhibited  p-state  pairing 
idea.  This  paper  reports  the  first  ail-atom-all¬ 
electron.  self-consistent  semirelativistic  energy-band 
calculations  for  TiBe2,  ZrZn2  (at  ambient  and  high 
pressures),  the  high-T,.  superconductor  ZrV3.  and  the 
ultralow- rc  superconductor  YAI2.  Total  and  partial 
(by  atom  type  and  /  values)  density  of  states  (DOS), 
calculated  Stoner-like  parameters  to  study  the  condi¬ 
tions  for  magnetism,  electron-phonon  coupling 
parameters  A,  and  Tc  values  and  their  behavior  under 
pressure  are  presented  and  used  to  provide  a  qualita¬ 
tive  understanding  of  these  phenomena  including  the 
possibility  of  superconductivity  in  ZrZn2  and  TiBe2 
with  the  application  of  pressure  without  invoking  p- 
state  pairing. 

II.  METHOD 

The  energy-band  structures  were  determined  self- 
consistently  using  the  linear  muffin-tin-orbital 
(LMTO)  method. :o  The  CIS  (dfl2)  structure  is  a 
very  closely  packed  structure  with  71%  of  the  volume 
filled  by  touching  spheres  and  the  site  symmetry  is 
high  for  the  two  constituents.  This  makes  the  Cl  5 
structure  especially  suitable  to  study  with  the  LMTO 
band  method  using  spherically  symmetric  potentials. 
The  basis  set  included  s,  p.  and  </  orbitals  ( -  2  ) 
for  each  site,  while  the  three-center  terms  included  / 
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Results  of  ab  initio  self-consistent  linear  muffin-tin  orbital  energy-band  studies  for 
EuMOjS,,  GdMogS,,,  SnMosS,,  and  SnMocSe,  including  ail  electrons  in  all  15  atoms/unit 
cell  are  reported.  The  large  charge  transfer  from  both  the  Eu,  Gd,  and  Sn  sites  and 
from  the  Mo  sites  to  the  chalcogens  is  shown  to  be  the  driving  mechanism  with  which  to 
explain  their  unusual  magnetic  and  superconducting  properties. 

The  unusual  magnetic  and  superconducting  tal  and  theoretical  interest.1  These  Chevrel- 

properties  of  the  ternary  molydenum  chalcogen-  phase  compounds,  with  general  formula  A/tMo«Ss 
ides  have  attracted  a  great  amount  of  experimen-  (or  Se9)  and  jr  between  1  and  4,  occur  for  a  large 
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As  described  elsewhere  in  these  proceedings,  the  Chevrel 
phase  molybdenum  chalcogenides,  MMogSg  (or  Seg) ,  exhibit  a 
number  of  unusual  properties  including  high  Tc  superconductiv¬ 
ity  and  the  highest  upper  critical  fields,  Hc2,  known.  They 
also  exhibit  magnetism  and  its  coexistence  with  superconduc¬ 
tivity  and  the  phenomenon  of  re-entrant  magnetism  at  low  tem¬ 
peratures  (see,  for  example,  Fischer,  1978,  and  references 
cited  therein) . 

From  the  theoretical  side,  the  complexity  of  the  crystal 
structures  has  made  difficult  even  non-self-consistent  energy 
band  studies  of  the  Chevrel  phase  ternaries.  Instead,  early 
studies  were  carried  out  with  approximate  molecular  cluster- 
tight  binding  (Mattheiss  and  Fong,  1977;  Andersen  et  al.,1978) 


2 Supported,  by  the  AFOSR  (grant  No.  76-2948),  the  N$F 
(grant  Nos.  DMR  77-22776  and  DMR  77-22646  and  through  the 
Northuestem  University  Materials  Research  Center ,  NSF  grant 
No.  76-80847 ) ,  and  the  COE. 
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I.  INTRODUCTION 


The  observation  in  lanthanum  of  a  high  superconducting 
transition  temperature,  T  ,  and  its  large  pressure  derivative, 
dTc/dp,  has  led  to  wide  speculation  that  2  pairing  mechanism 
involving  f  electrons  is  responsible.  Under  pressure,  Tc 
rises  sharply  from  6K  at  ambient  pressure  and  saturates  at  a 
value  of  nearly  13K  around  200  kbar.  In  addition,  at  53  kbar 
and  low  temperature,  a  phase  transition  occurs  which  is 


^ Supported  in  part  by  the  Air  Force  Office  of  Scientific 
Research  (Grant  76-2948 ),  the  Department  of  Energy  and  the 
Rational  Science  Foundation  (through  the  Northwestern  Univer¬ 
sity  Materials  Research  Center,  Grant  76-80847 ) . 
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Self-consistent  linearized  augmented-plane-wave  study  of  the  electronic  structure  and 
superconductivity  of  fee  lanthanum  under  pressure 

Warren  E.  Pickett* 

Department  of  Physics  and  Astronomy,  Northwestern  University,  Evanston.  Illinois  60201 

A.  J.  Freeman  * 

Department  of  Physics  and  Astronomy.  Northwestern  University,  Evanston.  Illinois  60201 
and  Argonne  National  Laboratory,  Argonne,  Illinois  60439 

D.  D.  Koelling 

Argonne  National  Laboratory,  Argonne,  Illinois  60439 
(Received  16  October  1979) 

We  report  the  resul  '  of  a  Unearned  augmented-plane-wave  calculation  of  the  electronic  structure  of  fee  La  at 
throe  lattice  constants  corresponding  to  ambient  pressure,  50,  and  120  Itbars.  The  Kohn-Sham-Gaspar 
approximation  for  exchange  and  correlation  is  used  and  the  potential  is  allowed  a  fully  non-muffin-tin  form.  The  / 
bands  lie  — 2-2.5  eV  above  the  Fermi  level  and  are  ~  I  eV  wide,  resulting  in  a  very  small  (0  05  electrons)  localized  / 
occupation.  Under  pressure  the  /  bands  nse  and  broaden  appreciably,  resulting  in  only  a  slight  increase  in  / 
occupation.  The  rigid-muflin-tin  approximation  for  the  eiectron-phonon  interaction  a  overestimates  the 
superconducting  transition  temperature  T,  by  40^-.  but  we  find  that  tne  drastic  increase  in  Tr  under  pressute  can 
'  be  attributed  primarily  to  changes  in  the  electronic  stiffness  7.  Structural  transitions  which  occur  at  25  and  53 
khan  may  be  related  to  changes  in  Fermi-surface  topology  which  we  find  to  occur  approximatelv  at  these  pressures. 


L  INTRODUCTION 

The  observation  in  lanthanum1  of  a  high  super¬ 
conducting  transition  temperature  Te  and  its  large 
pressure  derivative  dTjdP  has  led  to  wide  specu¬ 
lation  that  a  pairing  mechanism  involving  /  elec¬ 
trons  is  responsible.  Under  pressure,  T,  rises 
sharply  from  6  K  at  ambient  pressure  and  satu¬ 
rates  at  a  value  of  nearly  13  K  arouna  200  kbars, 
as  shown  in  Fig.  1.  Two  separate  questions  arise: 
(l)  Wiy  is  La  a  high-Te  superconductor  at  all, 
since  the  isoelectronic  elements  Sc  and  Y  are  not 
good  superconductors?  (ii)  What  is  the  cause  of 
the  large  pressure  derivative,  dTe  dP^l  K  10 
kbar,  which  is  among  the  largest  known  for  super¬ 
conductors? 

From  Fig,  1,  the  slope  aT^'dP  is  seen  to  show 
structure  at  25  and  53  kbars.  In  addition.  Balster 
and  Wittig1  have  measured  a  concomitant  anomaly 
in  the  low- temperature  electrical  resistivity  at  53 
kbar.  As  La  is  known3  to  be  fee  at  room  tempera¬ 
ture  in  a  wide  pressure  range  around  53  kbar,  they 
have  proposed  that  these  anomalies  arise  from  an 
Lsostructural  (fee  — fee)  phase  transition  which 
terminates  at  a  critical  point  below  room  tempera¬ 
ture.  Especially  considering  the  decrease  in  dTc 
dP  by  a  factor  of  4  (Ref.  1)  at  this  transition,  this 
revives  questions  concerning  the  complex  inter¬ 
relation  of  lattice  instability,  phonon  softening, 
and  high- temperature  superconductivity. 

The  superconductivity  of  La  was  attributed  by 
Rondo5  and  many  others**5  to  /-electron  mecha- 
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nisms  of  various  kinds,  with  the  large  pressure 
derivative  of  Tc  resulting  from  variation  with 
pressure  of  the  energy  of  the  /  state.  The  explana¬ 
tion  in  terms  of  /-electron  states  lying  above,  but 
very  close  to,  the  Fermi  level  £rwis  also  sup¬ 
ported  by  other  peculiarities  observed  in  La.  The 
temperature  dependence  of  the  magnetic  sus¬ 
ceptibility5  was  interpreted  to  be  due  tc  an  ef¬ 
fective  localized  moment  of  one- half  Bohr  mag¬ 
neton,  and  the  Knight  shift,  also  normally  tem¬ 
perature  independent,  was  found  to  increase  by 
20no  in  cooling  from  room  temperature  to  1.5  K. 
The  linear  specific-heat  coefficient*  was  found  to 
be  very  large  and  the  thermal- expansion  coeffi¬ 
cient3  was  found  to  be  negative  i.i  certain  tempera¬ 
ture  regions.  These  phenomena  all  seemed  to 
suggest  the  occurrence  of  significant  4/r-electron 
character  in  La,  and  it  has  been  argued1’10  that 
La  must  be  a  "4/- band"  metal.  In  addition,  the 
low  melting  point  of  La  had  led  Matthias  et  al.n  to 
postulate  an  /-electron  concentration  of  0.7  elec¬ 
trons/  atom. 

However,  the  development  of  reliable  tunnel 
junctions  led  to  the  unambiguous  assignment  of  La 
as  a  relatively  strong-coupling  d-electron  super¬ 
conductor,13*15  with  electron-phonon  coupling  con¬ 
stant  \=  0.8-0. 9.  While  this  partially  accounts  for 
the  large  electronic  specific  heat,  which  is  en¬ 
hanced  by  the  factor  (1  +\),  it  still  leaves  the 
possibility  that  the  pressure  enhancement  of  Te 
is  related  to  /character  in  the  electronic  struc¬ 
ture.  Ratto,  Coqblin,  and  CalleanP  have  attributed 
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APPENDIX  C 


Lectures  and  Seminars  presented  by  Arthur  J.  Freeman  ,  1979 

• 

Theory  Colloquium,  Oak  Ridge  National  Laboratory,  Oak  Ridge,  Tennessee, 

Feb.  7,  1979.  Talk:  "Local  Density  Approach  to  Bulk  Solids  and  Surfaces". 

Physics  Department  Colloquium,  Louisiana  State  University,  Baton  Rouge, 

Louisiana,  Feb.  8,  1979,  Talk:  "Electronic  Structure  of  Bulk  Solids,  Surfaces 
and  Chemisorption". 

Seminar,  Los  Alamos  Scientific  Laboratory,  Los  Alamos,  New  Mexico,  Feb.  2, 

1979,  "Re-entrant  Magnetism  in  Rare-earth  Superconducting  Compounds". 

Physics  Department  Colloquium,  University  of  Illinois  at  Chicago  Circel, 

Chicago,  Illinois,  Feb.  21,  1979,  Talk:  "Electronic  Structure  of  Bulk  Solids, 
Surfaces  and  Chemisorption". 

Symposium  on  Muon  Spin  Rotation  and  Associated  ProbLems,  April  12-14,  1979, 
Dogashima,  Japan,  Invited  Paper:  "Muon  Hyperfine  Anomalies". 

Talks  given  at  the  Institute  of  Physics,  Peking,  People's  Republic  of  China, 

April  17,  1979,  "Local  Density  Theory  of  Electronic  Structure  of  Solios"; 

April  18,  1979,  "Reentrant  Magnetism  in  Rare-earth  Superconducting  Ternary  Compounds 
April  19,  1979,  "Surface  States  and  Surface  Magnetism  in  Ferromagnetic  Transition 
Metals ". 

Talks  given  at  the  University  of  Nanking,  Nanking,  People's  Republic  of  China, 

April  23,  1979,  "Local  Density  Theory  of  Electronic  Structure  of  Bulk  and  Surfaces"; 
April  24,  1979,  "Reentrant  Magnetism  in  Rare-earth  Superconducting  Ternary  Compounds 

Talks  given  at  Fudan  University,  Shanghai,  People's  Republic  of  China, 

April  25,  1979  ,  "Local  Density  Theory  of  Electronic  Structure  of  Bulk  ana  .'urtaces"; 
April  26,  1979,  "Reentrant  Magnetism  in  Rare-earth  Superconducting  Ternary  Compounds 

Talks  given  at  Hangchow  University,  Hangchow,  People's  Republic  of  China, 

April  27,  1979,  "Local  Density  Theory  of  Electronic  Structure  of  Bulk  and  Surfaces". 

American  Physical  Society  Meeting,  Chicago,  Illinois,  March  19-23,  1979, 

Leader  of  group  session  "Band  Theory  of  Metals". 

Tnt .  tnnf  on  Layered  Materials  and  Intercalates,  University  of  Nijmegen, 

The  Netherlands,  Aug.  28-31,  1979,  Member  of  the  Organizing  Committee  and 
Chairman  of  the  2nd  plenary  session. 

ICM'79,  Munich,  Germany,  Sept.  3-7,  1979,  Talks:  "Magnetic  Interactions  in 
Chevrel  Phase  Ternary  Compounds"  and  "Surface  States,  Surface  Magnetization  and 
Electron  Spin  Polarize' ion  in  Ferromagnetic  3d  Metals";  Member  Int.  Program 
Committee;  Organizer,  Chairman  and  Invited  Speaker  at  Panel  Symposium  on 
Surfaces  and  Magnetism;  Speaker  at  Closing  Plenary  Panel  "Trends  In  Magnetism". 

Int.  Conf.  on  Mossbauer  Spectroscopy",  Portoroz,  Yugoslavia,  Sept.  10-14,  1979, 
Invited  paper:  "Electronic  Structure  of  Surfaces,  Surface  Magnetism  and  Surface 
Phase  Transitions"  (withC.  S.  Wang,  H.  Krakauer  and  M.  Posternak);  Chairman 
and  Organizer  of  Panel  on  Surfaces. 

Talk  at  ETH,  Zurich,  Switzerland,  Sept.  17,  1979,  Physics  Department  Colloquium, 
"Surface  Magnetism  in  Ferromagnetic  3d  Metals". 


Arthur  J.  Freeman 


American  Physical  Society  Talks  Presented,  March  1979,  Chicago,  Illinois 

"Near  Infra-red  Luminescence  from  CuCl"  (with  G.K.  Wong,  J.  Bajaj  and  C.W.  Chu). 

"Infra-red  Absorption  in  CuCl"  (with  F.  Wagner,  Jr.,  M.  Vjeinert  and  C.W.  Chu). 

"Self-Consistent  Electronic  Structure  of  Ternary  Chevrel  Phase  Compounds", 

(with  T.  Jarlborg). 

"Electronic  Structure  of  CuCl"  (with  C.S.  Wang,  T.  Jarlborg  and  M.  Weinert). 

"Electronic  Structure  of  Cu-  and  Ag-Halides"  (with  T.  Jarlborg). 

"isostructura 1  Phase  Transitions  and  Superconductivity  in  fee  La"  (with 
W .  g.  Pickett  and  D.D.  Koelling). 

"Local  Density  Functional  Study  of  the  O'  -  y  Isostructural  Phase  Transition 
in  Cerium"  (with  W.E.  Pickett  and  D.  D.  Koelling). 
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Lectures  and  Seminars  Presented  -  1980 

Physics  Department  Colloquium,  University  of  Indiana,  <B loomington ,  Indiana, 

February  12,  1980:  "Magnetism  and  Superconductivity  in  Ternary  Rare-earth  Compounds" 

Sanibel  Symposia,  Palm  Coast,  Florida,  March  4-20,  1930:  "Electronic  Structure  of 
High  Temperature  Superconductors;  C-15's  and  Chevrel  Phase  Compounds" 

International  Symposium  on  the  Physics  of  Actinides  and  Related  4f  Materials, 

Zurich,  Switzerland ,  April  9-11,  1980. 

Invited  papers : 

"f-Electrons  in  Solids:  Electronic  Structure  and  Properties  of  Actinides 
and  Rare-Earths"  ("Opening  talk) 

"f-electron  Bonding,  Electronic  Structure  and  Phase  Transitions" 

Chairman  of  Panel  Session:  "Localized  vs  Itinerant  f-electrons" 

International  Conferent  on  Physics  of  Magnetic  Materials,  Jaszoviec,  Poland, 

April  12-20,  1980.  ''Magnetism  of  Surfaces,  Interfaces  and  Modulated  Structures" 
(invited  paper);  Organized  and  Chaired  Symposium  on  Trends  in  Magnetism. 

Colloquium  Physics  Department,  ETH,  Zurich,  Switzerland,  April  23,  1980,  "Magnetism 
and  Superconductivity  of  the  Chevrel  Phase  Compounds" 

Colloquium  Physics  Department,  University  of  Basel,  Basel,  Switzerland,  April  28, 
1980,  "Magnetism  Surfaces,  Interfaces  and  Modulated  Structures" 

Colloquium  Physics  Department,  ETH,  Zurich,  Switzerland,  April  29,  1980, 

"Surface  States  and  Electronically  Driven  Surface  Phase  Transitions" 

Colloquium  Physics  Department,  University  of  Geneva,  Geneva,  Switzerland, 

May  2,  1980,  "Magnetism  and  Superconductivity  of  Chevrel  Phase  Compounds" 

International  Conference  on  Ordering  in  Two  Dimensions,  Lake  Geneva,  Wisconsin, 

May  28-30,  1980  "Generalized  Susceptibilities  and  Electronically  Driven  Surface 
Phase  Transitions  (invited  paper  given  by  H.  Krakauer) 

V  International  Conference  on  Hyperfine  Interactions,  Berlin,  West  Germany, 

July  21-25,  1980,  "Magnetism  of  Surfaces,  Interfaces  and  Modulated  Structures" 
(Invited  talk)  (unable  to  attend  due  to  illness  of  son) 

Yamada  Conference  IV  on  Physics  and  Chemistry  of  Layered  Materials,  Sendai, 

Japan,  September  8-10,  1980,  "Electronic  Band  Structure  of  Transition  Metal 
Dichalcogenides  and  Artificial  Layered  Coherent  Modulated  Structures" 

(Invited  paper)  (unvable  to  attend  due  to  illness  of  son) 

International  Conference  on  Ternary  Superconductors,  Lake  Geneva  Wisconsin, 

September  24-26,  1980  "Spin  Polarization  and  High  Critical  Fields  in  Chevrel 
Compounds"  < Invited  talk) 
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'folks  presented  at  the  American  Physical  Society  Meeting  in  New  York,  24-28  March  1980 
Invited  talks 

"Magnetic  Isolation  and  Superconductivity  in  Rare-Earth  Chevrel  Compounds" 

T.  Jarlborg 

i 

"Surface  Resonances,  Generalized  Susceptibilities  and  Electronically  Driven 
Surface  Phase  Transitions" 

H.  Krakauer 

Other  talks 

"Self-Consistent  Electronic  Structure  of  Coherent  Metal  Superlattices :  Cu-Ni" 
(with  T.  Jarlborg) 

'’Magnetism/Superconductivity  in  C15  Compounds  from  Self-consistent  LMTO 
Band  Calculations" 

(with  T.  Jarlborg) 

"Hyperfine  Anomalies  in  Muonic  Atoms  From  Breit  Interactions  and  Unrestricted 
Dirac-Fock  Solutions" 

(with  J.V.  Mallow,  J.P.  Desclaux  and  M.  We inert ) 

"Self-consistent  Electronic  Structure  of  the  Contracted  Tungseten  ''001) 

Surface" 

(with  M.  Postemak  and  K.  Krakauer) 

"Self-consistent  Electronic  Structure  of  A  ■’(111)  Films  and  Chemisorption 
Bonding  of  Oxygen  on  A ’(111)" 

(with  D.-S.  Wang,  H.  Krakauer  and  M.  Postemak) 


